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ABSTRACT 


The monograph is devoted to the study of dynamic processes in technical systems for various 
technological purposes. A new approach and methodology is proposed for a systemic, synergistic 
approach, taking into account the influence of energy fields of physical and mechanical effects, 
transformation and inversion of types of energy action. Models and equations of motion of dis- 
crete and continuous dynamic systems, dispersed media in the spectrum of their processing are 
considered. Changes in the parameters of subsystems are revealed: working mediums, mechanical 
systems, the processes of their interaction are investigated on the basis of taking into account 
their stress-strain state. The processes of grinding, sorting, mixing, compaction of materials and 
media are considered. The intensification of physical and mechanical processes, methods and means 
of their creation was achieved by the systematization and complexity of approaches due to the joint 
consideration of the mutual influence of the internal properties of subsystems. The processes of 
material processing by superresonant, subresonant and multi-mode parameters are investigated. 
Methods for determining effective parameters and modes of their operation are proposed. The 
processes of grinding, sorting, mixing, compaction of the processing medium in the field of vibra- 
tion load are described. As a result of the research carried out, new properties of the behavior 
of discrete-continuous systems under power load conditions are revealed. For the first time, the 
stresses and deformations of both working bodies and media were taken into account to create 
energy-saving vibration systems for various technological processes. The carried out scientific 
research makes it possible to obtain the laws of change in the state of dispersed media under the 
action of power loads by technical systems, new technological and design solutions were proposed. 


KEYWORDS 


Technical systems, technological media, grinding, sorting, mixing, compaction, cavitation, mo- 
dels, continuous and discrete parameters, force loads, stresses, deformations, energy, synergetics, 
vibration amplitude, vibration frequency, laws of motion. 
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CIRCLE OF READERS AND SCOPE OF APPLICATION 


It is intended for students, masters, post-graduate students, doctoral students of higher edu- 
cational institutions of a technical profile, and can also be useful for designers, research workers 
of industrial mechanical engineering. 

The processes of material processing by superresonant, subresonant and multimode parame- 
ters are investigated. Methods for determining effective parameters and modes of their operation 
are proposed. The processes of grinding, sorting, mixing, compaction of the processing medium in 
the field of vibration load are described. As a result of the research carried out, new properties of 
the behavior of discrete-continuous systems under power load conditions are revealed. For the first 
time, the stresses and deformations of both working bodies and media were taken into account to 
create energy-saving vibration systems for various technological processes. The carried out scienti- 
fic research makes it possible to obtain the laws of change in the state of dispersed media under the 
action of power loads by technical systems, new technological and design solutions were proposed. 

The monograph will be of interest to the countries of Europe, America, Asia. 
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INTRODUCTION 


|. Nazarenko 


Technological technical systems occupy a leading place in many technological processes in 
various sectors of the national economy. The classification features and the proposed criteria have 
been determined, according to which the assessment of the design and technological parameters 
of the existing technical systems of industrial engineering has been carried out. A method for 
studying technical systems for technological purposes has been developed, which provides for 
the implementation of analytical studies and the rational use of the regularities of changes in the 
internal properties of the machine-medium system based on the joint use of models with discrete 
and distributed parameters. 

The study of the dynamic parameters of the vibration crusher has been carried out and the 
dependences of the influence of the ratio of vibration masses on the grinding process have been 
established, and algorithms for calculating the parameters of a vibration three-mass crusher have 
been formulated. 

The physical and mathematical model of the screen is substantiated, on the basis of which the 
equations of motion of the material over the sieve are compiled and the distribution of the ampli- 
tude of oscillations of the material particles from the oscillation frequency of the sieve is obtained. 
The influence of the vibration frequency and the angle of inclination on the sorting efficiency have 
been established. An algorithm and method for calculating the main parameters of a vibration 
screen have been developed. 

It has been determined that the use of a vibration concrete mixer makes it possible to reduce 
the duration of mixing by 2.5—3 times in comparison with conventional free mixing concrete mixers. 
It has been proven that the process of preparing hard concrete mixtures and improving the quality 
of concrete mixtures occurs due to the destruction of defective aggregates and a more uniform 
distribution of the binder throughout the mixed mixture. 

An algorithm and method for calculating the main parameters of a vibration concrete mixer 
have been developed. 

A model of the hybrid system «vibration plant — medium» is proposed, in which active and reac- 
tive resistances to wave coefficients are taken into account. The results of calculations of changes in 
the parameters of vibration machines for the processes of compaction of concrete mixture confirm 
the implementation of accounting for the reactive and active resistance of the concrete mixture. 

Predictive mathematical models have been created to determine the stress state of the soil 
during compaction at characteristic points of the massif under the pipelines, which make it possible 
to obtain rational values of the parameters of the working bodies for soil compaction. An algorithm 
and method for calculating the main parameters of vibration systems for compaction of working 
media have been developed. 
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The assessment and substantiation of methods for studying the parameters of acoustic treat- 
ment of technological media have been done. The functional relationships between the acoustic pa- 
rameters of the cavitation apparatus and the rheological properties of the processing technological 
media have been revealed. The process of staged acoustic treatment of technological media is de- 
scribed. Analytical dependences are obtained for determining the main parameters of the cavitation 
apparatus and the regularities of the processes of staged acoustic treatment of technological media 
under conditions of energy minimization are established. An algorithm for substantiating the rational 
choice of the structural-parametric system «acoustic apparatus — technological medium» is proposed. 

The main states of reliability of a technical system and its elements are investigated by the 
example of vibration platforms for compacting concrete mixtures. A model for assessing reliabi- 
lity based on fuzzy logic, implemented in the MATLAB b2020 medium, has been developed. When 
processing statistical data on the MIBF of the propeller shafts of vibration platforms, the Weibull 
distribution law was established and the possibility of failure-free operation of such structural ele- 
ments was determined. 

In the course of the research, finite element models of vibration machines were developed 
to implement the transmission of energy in the vertical direction with a fixed frequency, with 
a polyphase direction of the action of external forces, with the spatial direction of movement in the 
vertical and horizontal planes. 

Investigations of the nature of distribution and numerical values of stresses and strains in the 
forming structure of the unit are carried out, depending on the angle of instantaneous action of 
the external force of vibrators and variable frequencies. 

otress concentrators in structural elements of machines for technological purposes are deter- 
mined, taking into account the workloads. 

The static and dynamic loads of the support contour of a truck crane have been investigated. 

The obtained research results are taken into account for making further design decisions when 
creating such machines. 
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ASSESSMENT OF THE CURRENT STATE OF PARAMETERS 
AND OPERATING MODES OF TECHNOLOGICAL 
TECHNICAL SYSTEMS 


|. Nazarenko, O. Dedov, |. Bernyk, A. Bondarenko, 
A. Zapryvoda, M. Nazarenko, |. Pereginets 
ABSTRACT 


The classification features of technological technical systems, the further development of the 
theory of the working process are given. The classification features are determined. Only the 
classification features that determine their purpose and the nature of the work they perform are 
considered. An assessment of the design and technological parameters of the existing techni- 
cal systems of industrial engineering has been carried out. New and improved evaluation criteria 
have been used. Reliable assessment and analysis made it possible to formulate more reasonably 
methods for solving the problem. Continuously improving, the technology leaves a reflection on the 
directions of development of technical system designs. On the basis of this, working hypotheses 
were formulated, the essence of which lies in modeling the medium and technical system as a single 
system with their own dynamic individuality. The system «technical system — medium» is represen- 
ted by a complex hybrid (mixed) dynamic system, in which the technical system is a system with 
discrete parameters, and the medium is with distributed parameters, and this system is reduced 
to the calculated one in the form of a system with discrete parameters. 


KEYWORDS 


Technical system, medium, classification, criteria, assessment, analysis, modeling, discrete, 
distributed parameters, hypotheses. 


11 DETERMINATION OF THE CLASSIFICATION CHARACTERISTICS OF TECHNICAL SYSTEMS 
AND PROCESSING MEDIUMS IN THESE STUDIES 


Technical systems can be processes and objects, depending on the nature of the connections 
between the elements of the system. The main goal of the development of technical systems is to elimi- 
nate the contradiction between the growing needs of society and the capabilities of today’s technical 
systems. Any technical system is primarily a physical object. And the correct choice of fundamental, 
that is, physical, foundations of functioning will lead to its viability and efficiency. The dialectic of the 
development of technical systems lies in the fact that new and well-known technical solutions are 
effectively combined in a newly created or improved system. 

The technical systems of industrial engineering include machines (equipment), which include a num- 
ber of elements, devices designed to perform certain technological operations and functions (Fig. 1.1). 
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O Fig. 1.1 Structural diagram of the machine 


For the development of general principles for the creation of machines, for the further deve- 
lopment of the theory of their working process, the classification features of these machines have 
been determined. The classification signs take into account the fundamental and design features 
of the structure of vibration machines, their purpose and the nature of the work they perform. 

In this work, the most significant two classifications of machines are used: constructive-tech- 
nological and structural. Structural and technological classification implies the determination of the 
use of a machine for performing a particular work process (crushing, sorting, compaction, etc.) 
with a note of the design features of their working bodies. 

Structural classification characterizes the general fundamental features of the structure of 
machines of dynamic action and the nature of the interaction of their working bodies with the 
processed medium. This classification creates the necessary prerequisites for the development of 
a unified theory of machines, guides and facilitates their further development. 

The parameters of work (the magnitude of the amplitude and the type of trajectory of oscilla- 
tions of the working body) of dynamic machines are determined by the nature of the working for- 
ces, the frequency of oscillations, a feature of elastic links and the ratio of the masses of individual 
parts of the system. Therefore, the structural classification in the monograph is determined by the 
signs of the degree of freedom of the dynamic scheme of the machine, the nature of its drive, the 
characteristics of elastic links. 

The type of work performed by machines and the peculiarities of the interaction of their work- 
ing bodies with the processed medium determine, on the one hand, the fundamental structure of 
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CHNOLOGICAL TECHNICAL SYSTEMS 
the machine, and on the other hand, put them in an equivalent relationship in assessing the working 
conditions. In this regard, one of the main classification features of a machine is the nature of the 
work it performs (features of the interaction of its working body with the medium being processed 
and the type of overcoming different types of resistance). 

All machines are divided into several groups that are fundamentally similar in structure and 
nature of interaction with the processed medium. 

The following operations, which are homogeneous with respect to the nature of external resistance 
to the working bodies, can be noted: the transfer of energy to different media during their processing. 

In the classical theory of oscillations, all dynamical systems are usually divided into systems 
with one, two, three, n and an infinite number of degrees of freedom (systems with distributed 
parameters). Such a distribution of dynamical systems presents certain conveniences with regard 
to their theoretical study. Dynamical systems with an equal number of degrees of freedom are 
described by the same number of differential equations solved by similar methods. 

By the type of drive, which is the third classification feature, the monograph examines vibration 
machines with a power drive. Currently, machines with a kinematic drive have little application and 
are not used in research. 

The fourth classification feature is based on the characteristics of elastic links, which deter- 
mine the nature of the restoring forces and internal resistances of the vibration machine. According 
to the characteristics of elastic connections, all investigated vibration machines are divided into 
machines with linear and nonlinear characteristics of stiffness [1]. 

A feature of these studies is taking into account the conditions of interaction with the process- 
ing medium [2], as a result of which there is a need to determine the classification features of such 
mediums. Such information is necessary for the justified direction of the formation of the algorithm 
and methods of technical systems of dynamic action. In such conditions, the key features of the 
classification are the state of the medium, the nature of the qualitative behavior of the medium, 
rheological properties and dynamic parameters. Depending on the ratio of the yield point t at pure 
shear to atmospheric pressure p,, the process medium in its state can be [3]: 


—hard plastic t/P, > 1; 
— liquid-plastic t/P, = 1; (1.1) 
— liquid at: t/P, <1. 


By the nature of the qualitative behavior under dynamic load in the medium (1.1), wave phe- 
nomena may appear or these phenomena are absent. Accounting for wave phenomena in the work- 
ing medium can be estimated from the ratio of the wave propagation time ¢.,, and the oscillation 
period 7 [4] under the following conditions: 


on (1.2) 
co] 
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Under the first condition (1.2), the process of oscillation of the medium can be considered slow 
and the elastic wave can be neglected. In this case, accelerations and deformations are determined 
exclusively by forces that are constant and unchanged in the direction of the action of external 
forces on such a medium. 

When the second condition (1.2) is satisfied, the motion of the medium is determined not only 
by forces, but also by elastic waves. Movement, speed and acceleration along the majestic change 
and depend not only on time, but also on the coordinate in the direction of which external forces 
act on such a medium. 

So, in the first case, inertial forces dominate, and in the second, elastic ones. In most ca- 
ses, real mediums to be investigated fall between these conditions (1.2). That is, as a rule, it is 
necessary to take into account the elastic, inertial and dissipative properties, which is accepted 
in this work. 


12 ASSESSMENT OF DESIGN AND TECHNOLOGICAL PARAMETERS OF EXISTING 
TECHNICAL SYSTEMS OF INDUSTRIAL ENGINEERING 


The design and technological parameters are assessed using functional criteria according to 
the appropriate method [5]. The methodology is based on the principles of parameter assess- 
ment (Fig. 1.2) of the same type of machines using the technical characteristics of existing 
technological technical systems. Subsequently, depending on the purpose of the technical system, 
the numerical values of the criteria are determined (Table 1.1) and, on their basis, conclusions are 
drawn on the formulation of research tasks for certain processes. 

The application of the above criteria, fulfilled for machines for crushing, sorting, mixing and 
compaction, are taken into account in the corresponding sections of the monograph. Their analy- 
tical form is shown in Table 1.1. 

Table 1.2 shows the numerical values of the main parameters of vibration crushers. 

Calculations according to this Table 1.2 obtained numerical criteria for evaluating the para- 
meters of vibration crushers, confirming significant differences (Table 1.3). 

The main parameters of resonant technical systems for sorting materials are presented 
in Table 1.4. 

As follows from Table 1.4, resonant screens, having a one-two-mass structure according 
to the screening scheme, are three and four-mass. The screen parameters are also different: 
the vibration amplitude (Table 1.4) has a wide range from 0.6 to 10 mm, foreign production — 
0-14 mm. They are also excellent in vibration frequency, which is reflected in one of the main 
indicators — vibration acceleration. In terms of the power of a foreign-made rotor for such parame- 
ters (amplitude and frequency of oscillations), they have lower values, which indicates a possible 
greater oscillation amplification factor, and therefore greater efficiency. 

Technical characteristics of gravity concrete mixers are given in Tables 1.5, 1.6. 
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@ Table 1.1 Criteria for evaluating technological and design parameters of technical systems for dynamic 
processes of processing technological media 


on Ff wo Fro 


Criterion 


Effect of mass on performance 

Impact of energy consumption on performance 

Effect of power on ground 

System synergy factor (efficiency) 

The rate of change in the stress state of the medium in time 
Speed 

Acceleration 

Dynamism of the system 

Mixing drum volume ratio 


The ratio of the mass of the medium to the mass of the technical system 


@ Table 1.2 Numerical values of the main parameters of vibration crushers 


The largest 


Tyne Productivi- Mass, Power, _ size of the 

yP ty, m*/h kg kWt input mate- 
rial, mm 

VSHCHD 440x800 22 1500 60 390 

(REC «Mechanort-tekhnika) 

VSHCHD 65x400 2 chambers 2.5 1400 14 50 

(REC «Mechanort-tekhnika) 

VSHCHD 130x300 1 1500 22 110 


(REC «Mechanort-tekhnika) 


VSHCHD 600x150 (aaa 1082 )7.5 120 
(Mining Academy in Krakiw) 

VSHCHD 600x800 47 2900 60 9 00 
(REC «Mechanort-tekhnika) 

VSHCHD 1200x1500 188 9400 110 1 000 
(REC «Mechanort-tekhnika) 

SVDA-1000 938 8116 480 1 000 


(Skochinsky Institute 
of Mining) 


Analytical 
dependence 


K ,=Plm 
k,=P/R 
k,=P/m 
K,=E/E 
v,, =do/dt 
k, =X,-@ 
k, =X,-o° 
k= aig 
K,=V,/V, 


K,=m,/M,, 


The smallest 
size of the 
source ma- 
terial, mm 
40 

10 

20 

15 

70 

100 


300 


Machine 
volume, 
mm 


10.92 


3.6 


2.89 
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@ Table 1.3 The value of the criteria for evaluating the parameters of vibration crushers 


Type 


VSHCHD 440x800 (REC «Mechanort-tekhnika) 
VSHCHD 65400 2 chambers (REC «Mechanort-tekhnika) 
VSHCHD 130x300 (REC «Mechanort-tekhnika) 
VSHCHD 600x150 (Mining Academy in Krakiw) 
VSHCHD 600x800 (REC «Mechanort-tekhnika) 
VSHCHD 1200x1500 (REC «Mechanort-tekhnika) 
SVDA-1000 (Skochinsky Institute of Mining) 


@ Table 1.4 Parameters of resonant technical systems for sorting materials 


Type 


VR-50 
VUR-80 
VKVS 
PEV-1-0,5xa 
PEV-2-0,5x5 
PEV-3-9x2,5 
79-15 

PVG-2 
GRL62 
PEV-2-4«12 
EV 


RMS-2,25 
RMN-4,5 
RS 

DV 

GC 
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Scheme 


One-mass unbalanced 


Two-mass balanced 
Two-mass balanced 
Two-mass 
Two-mass 
Two-mass 
Two-mass 
Two-mass balanced 
Two-mass balanced 
Two-mass 
Two-mass 
Two-mass 
Two-mass 
Two-mass 
Two-mass 
Two-mass 
Three-mass 
Four-mass 
Four-mass 


@ Table 1.5 Technical characteristics of Johnson-Ross gravity concrete mixers 
Standard size and loading volume 


yd? 
4.5 


m 

3.44 
4.59 
6.12 
7.65 
o.16 
11.4 


| 

3 440.7 
4587.6 
6 116.8 
7 646.0 
OAS 3c 


i 11 469.0 
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K1 K2 K3 K4 K5 
0.0015 0.38667 250.00 8.75 2.0147 
OC0TS OW ses 100008 a0 0.6940 
0.0007 0.0455 68.18 5.50 0.3460 
OOO Or 030s a4arey sree del aals: 
0.0019 0.7833 416.67 7.14 2.6184 
COGS a9 OST ASO Se OO Rolo 
0.0116 1.9542 169.08 0.01 — 
Amplitude of Vibration fre- Power, Maximum acce- 
vibrations, mm = quencyrad/i kWt leration, m/s? 
4.0 Joo 10 24 
4.0 97.4 10 38 
0.65 294 1.6 do 
0.60 314 OL: 60 
0.70 314 0.5 70 
0.80 314 2.0 80 
0.80 314 2.0 80 
4.0 78.9 30 24 
10 be es 2/ 
0.90 314 4.0 90 
OD lO 92.3...105 0...0 Bo ic0 | 
Oral Non 209 6 39...91 
0...14 NUS:<e09 6 36...01 
2) Me: 105...209 6 36:20 
0...14 105...209 3.9...0 | 36...91 
ee 105...209 Spec. e620 
Dinh 105...209 6 36...51 
Dele Go 20S de8 age 5 
10.12 105...209 ome: Ole e 
Drum rotation drive power, kW Number of Mass. t 
engine power _ total mixing blades ; 
1x29.8 29.8 3 6.804 
2x22.4 44.8 3 V2 701 
2x29.8 Od.6 4 14.515 
exon 74.6 4 16.33 
2x44.8 89.6 4 17.917 
2x06.0 We 4 19359 
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@ Table 1.6 Specifications of BMH Systems of The Roll Master Gravity Reversible Mixers 
Loading volume, m? Drum diameter, m Drive power, kW Productivity, m?/hour Mass, kg 


3 2.9 2x2e 85 12 000 
4.5 2.8 4x18 We 15 000 
6 3.19 4x22 165 18 000 
g 3.39 4x37 245 21 000 


The variety of currently used concrete mixers can also be explained by an unreliable compa- 
rative assessment of their technical, technological and operational parameters. Based on a review 
of the designs of prototypes of concrete mixing machines (Tables 1.5, 1.6), the key parameters 
were identified that affect the mixing process. According to the analysis of technical and economic 
indicators and parameters, the tasks of researching vibration concrete mixers are determined. 

Thus, studies using the above technique have established that for the processes of grinding [6, 7], 
sorting [8], mixing [9], compaction [10, 11], cavitation treatment [12-14], research of para- 
meters, modes and reliability machines [15-20] there is a significant difference in the numerical 
values of the criteria of the studied class of technical systems. The main reason lies in the use of 
simplified models, the absence of modeling and consideration of the joint movement of the working 
bodies of the technical system and the processed medium or material. 


13 DISCUSSION OF RESULTS AND IDENTIFICATION OF MAIN DIRECTIONS OF RESEARCH 


For the development of general principles for the creation of technological technical systems, for 
the further development of the theory of the working process, classification signs are determined. 
The monograph considers only the classification features that determine their purpose and the 
nature of the work they perform. This has some limitations of the existing classification signs. 

An important point in the analysis of existing works is the assessment of the design and tech- 
nological parameters of the existing technical systems of industrial engineering. New and improved 
evaluation criteria have been used. Reliable assessment and analysis made it possible to formulate 
more reasonably methods for solving the problem. Continuously improving, the technology leaves 
a reflection on the directions of development of technical system designs. On the basis of this, 
working hypotheses and research methods were formulated, the results of which are presented in 
the corresponding sections of the monograph. 

Working hypotheses for researching technical systems for technological purposes: 

— consideration (modeling) of the medium and the technical system as a single system with its 
own dynamic personality. The embodiment of this principle guarantees the movement of the studied 
systems in a given or established mode of operation; 

— «technical system — medium» system is a complex hybrid (mixed) dynamic system, in which 
the technical system is a system with discrete parameters, and the medium is with distributed 
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parameters, and this system is reduced to the calculated one in the form of a system with discrete 
parameters, in which wave phenomena of the medium and are represented by contact force; 

— an increase in the efficiency of technical systems for various technological purposes is 
achieved through the efficient use of energy supplied to the medium both on the main and on the 
Superharmonic resonance mode of oscillations of the system under study, as well as with a variable, 
time-controlled mode of operation; 

— reliability of the machine is ensured by a rational combination of the effect of shock and 
vibration at low frequencies, as well as the use of reliable vibration exciters. 

Research methodology for technical systems for technological purposes: 

— analytical research and rational use of the regularities of changes in the internal (elastic- 
inertial and dissipative) properties of the «machine — medium» system; 

— conducting experimental studies based on the determination of rational operating modes of 
technical systems for various technological purposes: 

— reliability assessment based on a small number or absence of failures: 

— study of the stress-strain state of metal structures of machines of static and dynamic action; 

— establishment of regularities in the distribution of normal and contacting stresses for diffe- 
rent shapes of the cross-section of elements of metal structures. 

The proposed method of compiling a mathematical model [1, 2] is based on the idea of con- 
sidering at the first stage a separate medium, the stress-strain state of which is under the load of 
the force of any law of its change. 

The next stage of the analytical solution of the motion of the medium is based on the premise that 
the structure of the dependence reflecting its state at any moment of motion consists of two compo- 
nents — one is capable of accumulating energy that passes from one form to another (reactance), and 
the other is energy dissipation (active resistance). The resulting formula has a discrete form in structure, 
however, its constituent parameters and characteristics take into account wave processes in the medium. 

The final stage is the drawing up of the equation of the joint motion of the vibration sys- 
tem «machine — medium» by the appropriate combination of reactive and active components of 
the resistance of the machine and the medium, for which the method of vector diagrams is used. 


CONCLUSIONS TO SECTION 1 


1. For the development of general principles for the creation of technical systems for technolog- 
ical purposes, the further development of the theory of their working process, classification features 
and proposed criteria have been determined, according to which the assessment of the design and 
technological parameters of the existing technical systems of industrial engineering is carried out. 

2. Working hypotheses of research of technological technical systems are formulated, the 
implementation of which is a guarantee of the reliability of the obtained results of the movement 
of the systems under study in a given or established operating mode. 
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3. A method for researching technological technical systems has been developed, which pro- 


vides for the implementation of analytical research and the rational use of the regularities of 
changes in the internal properties of the machine-medium system based on the joint use of models 
with discrete and distributed parameters. 
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RESEARCH OF PROCESSES OF PRODUCING MATERIALS 
BY TECHNICAL POWER LOADING SYSTEMS 
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S. Oryshchenko, O. Fedorenko, S. Tsepelev 
ABSTRACT 


An increase in the effective working process of the crusher is achieved when it operates in the 
vicinity of resonance, i.e. at a high value of the coefficient of resonant amplification of oscillations. 
The optimal values of the amplitude of the oscillations of the masses have the greatest values, 
provided that the coefficient of regulation of elasticity in frequency n is in the range of 0.6<n<1. 
Within the frequency range of 14.5 Hz<f<26 Hz, material is effectively crushed. The numerical 
values of the coefficient are established, which determine the rational ratios of vibration masses. 
The studies carried out and their results indicate that the main stresses causing the destruction of 
the material are shear stress. Analytical equations are obtained that reflect the picture of energy 
consumption for the destruction of the material. Such studies are planned as a continuation of the 
topic under consideration through experimental studies and consideration of several members of 
the series and optimization of the parameters of crushing machines. The proposed approach to 
studying the energy characteristics of material destruction in the crusher chamber with a guar- 
anteed zone of stability of parameters in the resonance zone can be used for other processes. 
such processes include the destruction of materials during dynamic cutting of soils, processing of 
materials in mills and sorting of materials, which are widely used in Ukraine abroad. An algorithm 
and method for calculating the main parameters of a vibration crusher have been developed. 


KEYWORDS 


Vibration crusher, material, grinding, physical, mathematical model, equation, resonance, para- 
meters, destruction, stability zone, stress, deformation, energy consumption. 


2.1 SELECTION AND JUSTIFICATION OF PHYSICAL AND MATHEMATICAL MODELS OF TECHNICAL 
SYSTEMS FOR MATERIALS GRINDING PROCESSES 


Vibration jaw crushers are a relatively new direction in the development of crushing techno- 
logy, which are highly energy efficient and are used in various industries, including construction. 
The most promising of the known vibration jaw crushers are crushers with two movable jaws and 
crushers with a heavy body and a pendulum. On the other hand, the creation of such machines 
is hindered by the severity of providing a resonant mode of operation in which the efficiency of 
crushers increases significantly. 
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Let's consider a model of a working fluid in terms of interaction with a machine. Let's consider 
the layout of the elements of the vibration system under study in the following order (Fig. 2.1): 
vibration exciter — machine body — working medium — striker. Let's take: m — total mass of the ma- 
chine; c — coefficient of elasticity of elastic connections; m,,— mass of the material; F,(¢) — driving 
force; x, — movement of individual elements of the crusher structure; D — average diameter of 
the input material. 





O Fig. 2.1 Layout of vibration system elements 


The equation of motion of our system will be an equation of the form [1]: 


ox 1 0°x 
cad een (2.1) 
ou’ c* ot® 
where c — speed of the wave propagating in the medium. 
The solution to equation (2.1) will be the dependence [1]: 
x(ut)=[ Aeas r+ Bsin2r snot (2.2) 
C lf 


Coefficients A and B are determined by the boundary conditions and, by their physical content, 
reflect the vibration amplitudes in the corresponding contact zones. 
Forces acting on the contact between the crushing plate and the material: 


—F., =F, sinaat — (m,X, + 2X,C,— X.0, + MyXp + Xp0,— X40; + 2CX—205X,+X4C,), (2.3) 


cr.m 


where F.._ — reaction force of the material from the action of the crushing jaw (equal to the 


cr.m. 


crushing force with the opposite sign); 
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F, sin oot — (m,X, + 2X,0, — XpC, + MyXp + XpC, — X40, + 26jXq— 2CXq + XqC; J} 
— crusher motion equation. 


On the other hand, the expression for the reaction of the material to the crushing plate can 
be written as follows: 


Fm = OS =E5 (**) (2.4) 
: Ou 


Equating expressions (2.3) and (2.4) let's obtain: 


x . m,X, + 2X,C, — XC, + MX, + 
F5( =Fsinnt-| i oT! et ee , (2.5) 
ou +Xo6,— XC, + CC :Xq — 20jXq + Xqlon 


Let’s find the unknown derivatives dx/du and x = ox*/0°t by differentiating equation (2.2): 


2 ={-ASsin@v + 82c05% sno, (2.6) 
ou Cc .)hOU C C 

2 
Br =| Acosy + Banyo st (2.7) 
Ot C C 


To find the displacements of the first and third masses included in equation (2.5), let’s use the 
solution in the form of the following functions: 


‘ [t] =A, sin ot; on 


X,[t] = A, sino. 


Thus, an equation can be derived to determine the vibration amplitudes of the contact zone: 
r 

B, = Ata( “| (2.9) 
C 


Ac, -2A,c, + A,c, —F, —w°Am, 


2c, — w°m, + E5( © }an{ 
C C 


As can be seen from equation (2.10), the effect of the material on the dynamics of the mouse 
is determined by a parameter ES(w/c)tan(rw/c) that can be turned into the following form [1]: 


A =- (2.10) 
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m.0° (‘) (2.11) 
W 


where m — the mass of the material. 

The ratio (toy /w) = y is called the coefficient of added mass. It takes into account the dif- 
ference between the mass m,, and its static value m,,, [1]. 

When designing any machine, great attention should be paid to the rational selection of kine- 
matic parameters. The kinematic parameters mainly include the geometric characteristics of the 
machine as a whole and its individual parts. 

One of the features in the design of vibration crushers is the presence of vertical movements 
of the crusher relative to the falling material. 

Expanded, the dependency for defining the material path will be as follows: 


| - g(sinlo.] — 1, cosa) eee 


a (2.12) 
Af m,.@ 


where ll, — coefficient of friction; m,. — the mass of the vibration parts of the machine (mass of 
the crusher). 

Knowing the total length of the crushing plate, let's find the number of cycles of application of 
loading to the material during its stay in the crushing chamber (h/1, ) = n,. 

On the other hand, the required number of load cycles for material destruction is calculated 
by the formula [2]: 


n’ = (5,/8,), (2.13) 
where 6, — depth of entry of the crushing plate into the body of the crushed material, at which the 
material is destroyed, m; 5, — depth of entry of the crushing plate into the body of the material at 


a given energy of a single impact. 
The value 5, depends on the dependence [2]: 


5, =(2 V)/(¢,k,5,); (2.14) 
where E. = mv*/2 — impact energy, J; v — speed of application of the load to the material, m/s; 
c, — wave propagation speed, m/s; k, — dynamic factor; S, — contact area, m*. 

2.2 INVESTIGATION OF THE STRESS-STRAIN STATE OF THE MATERIAL 


Study of the stress-strain state of a solid material. 
The stresses arising on the platforms of the selected element are shown in Fig. 2.2. 
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O Fig. 2.2 Working body model 


A complete description of the stress-strain state in the linear theory of elasticity (spherical 


coordinates) includes the following dependences [3]: 
— static equations: 


00,,, : 1 OT, 19%,, 
dp psin(o) 00 eas ey) 
OT 1 00, 19% 1 
Ee ep i ie eek N: 
dp psintg) 00 p ag (toe + 2cta(o))) + 
OT 1 OT 1 00 1 
P i) OP _ f): 
op psinl) 00 *D dp +5 (50 ~ Gu) ota(9)—8t9p) +L =O 


— deformation equations: 
































aa gretweate)} 


itm, a(v) _ a(w), tau 
Toa 5 sin(p) 00 ne 0 ae 0) pop 
1 1 07. ow 
(2 ( sin(@)) + a 


p sin() 
— stress equations: 





Yeo 


O,, =2GE,, +E; T5 = GY ,9: 
Og, = 2GEg, + AB; T,, = GY: 


Ou, = 26E,,, + he: To = OY 693 
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(2.15) 


(2.16) 


(2.17) 
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where p, 8, @ — suitable coordinates for the placement of the stressed body element; u, v, w — move- 
ment in the direction of the corresponding coordinate axes; 6,,,09,0,, — normal stresses at the 
respective sites; T,5,T,,:T 9, — Shear stresses at the respective sites; M, N, L — volumetric for- 
ces; G = E/2(1+ v) — shear modulus; v = 4/e(A + G)—Poisson’s ratio; A = Ev/(1+ v)1—-ev)) — 
Lame constant; e = 1/p* olp*u)/dp + (1/psin(@)) (Colw sin(@)) /d@) + ov/d@) — volumetric relative 
deformation; €,,, E99, €4 — linear deformations; y,, Y,;, Yer — Shear deformation. 

To solve the problem, let's use the displacement method. 

Combining equations (2.15)—(2.17) by alternate substitution, and performing some mathema- 
tical transformations, let’s obtain the equation of elastic balance of the body in displacements [4]: 


(1+ 26) psin(@)de/dp — 26| dé, /00 — a(E, sin(@)) /00| = 0; 
(A. + 26)(1/sin(p))de/00-26| a, /ag-a(pE, )/dp |= 0: (2.18) 
(4. +2G)sin(p)de/dq — 26] d(pé, sin(p)) /ap — a, /26 |. 


where &,,59.5; — components of rotation, determined by the following formulas: 


= 1/2psin(p)| aw/00 _ a(vsin(@))/d9 |; 
E, = 1/2p| ou/dp — 0(pw)/ap |: (2.19) 
Eo = /2psin(p)| a(pv sin(@)) /op ~ au/00 |. 
Therefore, to solve this problem, it is necessary to find the displacement functions. These 
functions must satisfy three main conditions [4]: 
1) be harmonious or biharmonic functions: 
2) transform equation (2.18) into an identity; 


3) stresses, expressed in terms of displacement functions, must satisfy the compatibility equa- 
tions on the surface, which are written as follows [5]: 


O,, =O,,/ + T 9M + T 9M 
Goy = Tog! + GyglM + Topi (2.20) 
Ooy = Tol + TygM + Oyo, 


where /, m, n — direction cosines. 
In expanded form, the displacement functions for a spherical coordinate system are written 
as follows [6]: 


2Gu, = 4(1- v)I (x, cos8 + x, sin®)sin + x, cosy | —0M/dp; 
2Gw, = 4(1- v)I (x, cos0 + x, sin®)cos@ — x, sin @ | —(/p)(oM/d@); (2.21) 
2Gv, = 4(1- v)| (=x, sin@ + x, cos®) | —(1/psing)(oM/28), 

where M=¥, + p| (x, cos6 + x, sind)sing + x, cos @ |. 
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In this problem, the body forces are insignificant, so they can be neglected. 

Therefore, from the above, it is possible to conclude that the solution of equations (2.21) is 
reduced to finding three harmonic functions ¥,,, Xo, X3- 

such functions are volumetric spherical functions, which, in general, can be written as follows: 


4 = VA, p’ +B, p "NC, cos(me) + D,, sin(mO} CL, Pl LH) + Fn (11), (2.22) 
where tt =cosl@] — module; P_ (i) — Legendre function of the first kind, degree n and order m; 
Q__{u) — Legendre function of the second kind, degree n and order m; @ (u) — Legendre function of 
the second kind; P,(u) — Legendre function of the first kind; A,, B,C, Das by ms Fm — constants. 

Now, knowing the displacement by formulas (2.16), (2.17), it is possible to find deformations 
and stresses. Let's consider the most atypical condition under which n1—n3, m1—mée are diffe- 
rent integers, for this —1< <1. For subsequent calculations, let’s take: n1=4, m1=1; n2=5, 
m2=2; n3=6, m3=3. Then equations (2.22) can be written as follows: 


Xa(mint) = -((p%, +B,)(sin[ 8 JH, + cos[ 0 |M,) /1 2/1-p° 0°] x 
-32 + 230u? —210u* + 


XI} 4 5u(3 10? + 7g) hig? 30 (3 —10u? + tae op ; 
1-y 
(p'"'A, +B, )(sin[ 20 ]H, + cos[ 20 |M, ) 
X2(m2,n2) — ate) x 
64-686" +1260" —630u° +105.(—1+w2) (-1+ 31a E+ 
- (fea a 
+ 210 (14?) (-14-3y?)L,, 
Laima) = -{(p4, +B,)(sin['36]H, +cos[30]M,) /4(1- 2)” "| x 


512+ 7326? — 20286‘ + 
+20370p° - 6930u° + 315u(—1+?) x 


(-3+1 1p) a 


Study of the stress-strain state taking into account the rupture zone in the body 
of destruction. 

The realistic crack shape is an elliptical shape. In this form, by changing the ellipticity, it is 
possible to realize different degrees of curvature of the crack boundary. 

The solution to this problem is complex and very voluminous, due to the use of biharmonic 
functions. Therefore, for further research, let’s use the technique proposed for solving the stress- 
strain state of a body near the rupture zone, which has the shape of an ellipse [3, 4]. 


Fag +630 (—1+p2) (-3+1 11?) by, |. (2.23) 
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In general, the problem is divided into two separate ones: 

1) an elliptical crack under internal pressure: 

2) an elliptical crack under shear conditions. 

For a better understanding of the problem statement, let’s present a diagram of the location 
of a plane crack in a three-dimensional body (Fig. 2.3). 

The limiting conditions for solving the problem of an elliptical crack under the action of loading 
are as follows: 


o=-F, €=0, w=0, 0=0. (2.24) 


To solve the problem, the gravitational potential outside a homogeneous elliptical plate of the 
following form was used [4]: 


el x? i 2 dJ 
=—||—— — —1|—., (2.25) 

* sll eater J lin 
where V(I) = |S (a + J)(b? + yr — function of J; J — continuous integration; a, b — respectively 


large and having semi-axes of an ellipse. 
Constant A is found by the following formula [4]: 





A= —ab°F/4uE(m), (2.26) 


n/2 
where E(m) = J a (sin[ @ ]) do — complete elliptic integral of the second kind; @ — ampli- 
0 


tude of the integral; m= sinlo] — modulus of the integral; o — modular angle. 





O Fig. 2.3 Scheme for determining the stress-strain state 
in the presence of a rupture in the crushing material 
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After some transformations and calculations, the displacements from the action of normal 
forces through the elliptical coordinates will have the following form [4]: 


Bao) eee | 


|| poeta yemey,aesmecomn]| 1 lle ss)oo]" 
a el en’, |e 7 


: abe ie o)(E- 0)(b +8)" 


(1-2v)E[u, +] Pee —2(1—y) 
9 Az dnu, cnu, 
w= a ge [é (ab — ot) —a°b? (0+ 3) (a +b) 00 | ; (2.27) 
+ ——— -— 


We 1/2 
a(&—0)(&-8)(a° +8)" (b° +6) 
where €, o, 0 — elliptical coordinates; € = const — family of ellipsoids; o — hyperboloid with one 
cavity; © — hyperboloids with two cavities; F [ u, | = J dn*BdB — function; u, — variable; m’ — argu- 
0 


ment of elliptic functions (modulus of the integral), (m’ +(m)° =1; snu,, cnu,, dnu, — elliptic 
functions of C. Jacobi. 

Now let’s return to solving the general problem of a body with a crack. Using the principle 
of superposition of the stress states of two problems, namely, the stress-strain state of a body 
without a crack and in the vicinity of a crack, it is possible to write the following equations: 


Ue =U, +U,,; CO, =O, +9, (2.28) 


where u., o,— displacement and stress, respectively, calculated for a body without a crack; u_,, 0, — cor- 
responding displacements and stresses calculated for bodies with a rupture zone. 

since the displacements indicated in formulas (2.27) were determined for elliptical coordi- 
nates, and preliminary calculations were performed for a spherical system, therefore, when per- 
forming further calculations, it is necessary to switch to a single coordinate system. Let's accept 
the spherical coordinate system as more convenient for further operation. 

After simple transformations of the system of equations (2.27), let’s obtain the value of the 
displacements in spherical coordinates: 


7 A (2%r° —70r“p + 105r°p* 
”  8r’ | — B4r°p? + 35rp* — Bp° 


"|[-8+(-5+4vjou[20)) + O[>-rP 7 


2A \aettanead 
= +cVv 


ee ~ Isin2ol+O0fo—rl’. (2,29) 
ap —B4r°p’ + 35rp* — 6p” | el oa 
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It should be noted that these equations were obtained on the basis of the expansion of formu- 
las (2.27) in power series. The calculations took into account the condition — & =~ p*; 0 = 0 =O. 


To obtain the stress components, it is necessary to substitute formulas (2.29) into the equa- 
tion for determining the stresses. 


Finally, the stresses in the body in the presence of a rupture zone will be equal: 
4 3 2.2 
Ela eres 
42r° (14 vy) —85r* (7+ 4v)p + 105r? (5 + 2v) p* — 
i= 372A —4er* (13+ 4v)p° + 70r(4+v)p' -3(19+ 4v)p° + 
—1/ pr 14r° (14 v) — 105r“p — 35r? (—7 + 2v)p* + 
+3] +14r?(-19+ 8v)p* —70r(-24v)p'+ — |cosl2o! 
+(—29 + 16v)p° 
_ ~ 70r'p + 105r'p" a + 4y + 3cos[2@]) — 
—84r°p* + 35rp* — 6p” 
42r° (14 v) —85r* (7+ 4v)p + 105r° : + 2v)p* — 
aa — 42r* (13+ 4v)p* + 70r(44 v)p* -3(19 + 4v)p” + 
—X 14r° (1+ v) —105r“p — 35r’ ( —7+2v)p* + 
+3] +14r°(-19 + 8v)p* — 70r( aed )p'+ — |coslee] 
+(-29 + 16v)p" 
3G (21r° — 70r‘p +105r%p° — 84r’p? + 35rp* —6p°) + 
42r° (1+ v)—35r* (7 + 4v)p + 105r* (5 + 2v) p* — 
bee’ (13+ 4v)p° + 70r(4+v)p* hin on 
a) 4 32 
Pee = a5" els ween aula J. 
+ 14r° (14 v) —105r“p — 35r’ (-7 + 2v)p* + ) |cosl2q] 
3A| +14r* (19+ Bv)p* —70r(-2 + v)p* + 
+(-29 + 16v)p° 
478 21r° (2 + vy) —35r* (5 + 4y)p + 


Tyg =O yy = az] + 31509 (14 v)p? —42r* (7 + Bv)p*+ [sinl2l;t,, =0. (2.30 
+ 35r(-4+5v)p" +9(3—4v)p° 
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Thus, knowing all the components of stresses and displacements for two stress-strain states, 
it is possible to determine the total energy of destruction. 

Based on the superposition of equations for determining displacements without a crack and 
equations for determining the displacements of a body with a crack, the corresponding graphs were 
constructed [4]. The graphs of the dependence of the principal stresses on the size of the major 
semiaxis of the crack ellipse are shown in Fig. 2.4. The graph of stress versus zenith angle @ is 
shown in Fig. 2.9. The change in stress op depending on the radius ris shown in Fig. 2.6. 
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O Fig. 2.4 Graph of stress dependence on the size of the major semiaxis of the ellipse: 
a — op-—a; b — 08-4; ¢ — Tp@—a 





O Fig. 2.5 Graph of stress changes depending on the angle of application of the force @: 
a— Op-; b- Tpg- 
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O Fig. 2.6 Graph of changes in normal stresses o, depending on the radius r 


As follows from the graphs (Fig. 2.4), the dependences of stresses on the size of the major 
semiaxis of the ellipse have a linear characteristic within the framework of the calculations, the 
numerical values of which were used to determine the parameters of the crusher. 

The significant influence of the angle of application of the force on the change in stresses is 
noted in the range of 0.3—0.5 (Fig. 2.5). 

A more complex law of variation of normal stresses depending on the radius r is noted in the 
range of 9-12 mm (Fig. 2.6). 


2.8 RESEARCH AND ANALYSIS OF MATERIAL GRINDING MACHINE PARAMETERS 


A mathematical model of an experimental jaw crusher. The physical model of a vibration jaw cru- 
sher is shown in Fig. 2.7. When considering the physical model, the structural elements of the machine 
take an elementary form, the gravity forces of individual links are applied at certain points, the forces 
of resistance to destruction are presented in the form of a total force applied to the jaw center [7, 8]. 





O Fig. 2.7 Physical model of a three-mass vibration crusher 
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When compiling the equations of motion of the vibration system of a jaw crusher, let’s use 
Hamilton's variational principle. It establishes that the variation in kinetic and potential energy 
plus variation in the work of non-conservative forces during any time interval should be equal 
to zero [9-11]. 

To compose the equations of motion, let's consider the structural diagram of a vibration 
jaw crusher. The investigated jaw crusher is made in the form of a three-mass resonance sys- 
tem (Fig. 2.7), in which rectilinear vibrations are realized and the dynamics of which occurs ac- 
cording to a three-mass scheme. 

The vibration exciter 1 with the base plate is the first or active mass, the movable body 3 
is the second mass, the middle plate with inner armor 0 is the third mass. Each mass performs 
horizontal oscillations along the x-axis along the generalized coordinates, respectively, and has its 
own inertial parameter. 

The active mass is set in motion by the centrifugal force arising from the rotation of the unba- 
lanced masses of the vibration exciter. The crusher is adjusted in such a way that during operation 
the first and third masses oscillate in phase, and the second in antiphase, that is, when one cham- 
ber approaches in another, the jaw diverge. 

All masses of the crusher are connected in pairs by elastic systems 2 and 5 (Fig. 2.7) with 
stiffnesses and. The vibration machine is isolated from the foundation by vibration dampers with 
stiffness on which the second mass rests. Resistance forces act in the system, but these forces 
can be neglected to simplify calculations. 

The system of equations of motion for the model using the Lagrange equation of the second 
kind will take the form [5]: 


M,X, + 2x,C, — XpC, = F, sino; 
(m, + km,,)Xp + XoC, — X,C, + 205X» — 2C,X, + X96, = 0 


m, +km_, )X, + 2CX, —2C,X, = 0; 
(s+ hit Jia + 2s — BOX (2.31) 
ee S CX, 

x, >0; x, >0; x, <0; 

X5 + X35 = ore 
where F, — centrifugal force (F, = mr); r, — eccentricity, distance from the center of gravity of 
the mass m, to the center of the axis of rotation, m; « — angular speed of unbalance rotation, s7': 
m,, M,, M, — the first, second and third mass of the crusher, kg; c,, c, — stiffness of the springs 
connecting the masses of the crusher; x,, x, and x, — corresponding displacement of the masses 
of the crusher, mm; c, — stiffness of the connections of the system connecting the machine to the 
foundation; F,, — crushing force formula, N; ¢ — relative compression of the material; D_.. — the 
maximum size of the crushed material, mm. 
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The solution to the system of equations (2.30) and (2.31) is obtained in the form: 


x, =X, sinot; 
x, = X, sinot; (2.32) 
X, = X, sino, 


where X,, X,, X, — amplitudes of mass displacement. 

Substituting (2.32) into (2.18) and carrying out the appropriate transformations, let's obtain 
the dependences for determining the values of the amplitudes of the displacement of the mas- 
ses X,, X5, X3. 

The energy consumed by the electric motor is spent on the work on crushing the material, on 
the potential energy of deformation of the elastic system, on overcoming friction forces, on heat 
losses, on the deformation of the crusher parts and other losses. The greatest energy consumption 
during the operation of the crusher is spent on crushing the material, the kinetic energy of the 
crushing plates and the potential energy of deformation of the springs. 


F (wo? (—c, + w*m,)m, + ¢, (2c, — o’m,) + 2c, (c, — a (m, + m;))] 


X,= 
sf eee 
wo (c, — om, ) mM, + " +o (m,+2m,) = 
c; (2c, wm, )+ @*m, ~¢, + +C, 2c, + 
+26, 
+0° (M+ M;) ~2¢, ( + | 
+0 
+2(m, +m, ) 
— c,F, (2c, — w°m, ) 
a alee 
oO (c,-@ m, )M, + +o° (m,+2m, ) 3 
c; (2c, wm, )+ om, —¢, + +C, ~2¢, + 
+20, 
+0° (M+ M,) —2c,| (im, + 
+2(m, +m, ) 
X,= ot | (2.33) 


; a 
=O m, + 
© (—c,,+ @'m, )M, + +o°(m,+2m,)} ° 
c; (—2c, + @°m, )+ w*m, *. +0, 
Fees) 
—«* (mM, +M,) +C, 
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Thus, the total energy consumed by the crusher for crushing the material will look like this: 
Fy =O+P, (2.34) 


where C — kinetic energy of the system; P — potential energy of the system. 

Part of the potential energy is spent on the destruction of the material. Knowing the general 
displacements arising in the material during its destruction and applying the condition of the mini- 
mum potential energy of deformation, it is possible to write the following: 





Vv 1 
U, = ats Der eaae ee e+ Ae + Yoo + th) sinodpdedg. (2.35) 


Knowing the energy consumption for the destruction of the material, one can proceed to 
determining the energy consumption for the process as a whole. 

Let's consider two vibration-driven crushing systems shown in Fig. 2.8. 

The vibration jaw crusher (Fig. 2.8, a) is made in the form of a three-mass resonance system, 
in which rectilinear vibrations are realized and the dynamics of which occurs according to a three- 
mass scheme. The vibration exciter 1 with the base plate is the first or active mass, the movable 
body 5 is the second mass, the middle plate with the inner armor 4 is the third mass. All masses 
of the crusher are connected in pairs by elastic systems 2 and 3. 




















O Fig. 2.8 Design diagram of vibration crushers: a — vibration jaw crusher with 
three vibration masses; b — cone crusher with vibrators on the body 


Vibration cone crusher (Fig. 2.8, 6) consists of vibration exciters 1 installed on a stationary 
cone / of the crusher. The fixed cone 7 is connected to the crusher body 4 through the adjusting 
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ring 6 and the elastic system 2. The crusher body 4 is installed on the foundation through the 
elastic system 3. The movable crusher cone is driven through the drive 5. 

The energy consumed by the electric motors of the crushers (Fig. 2.8, a, 6) is spent on the 
work on crushing the material. Additionally, energy costs are spent on deformation of the elastic 
system, overcoming friction forces, heat losses, deformation of crusher parts and other losses. 
However, the greatest energy consumption during the operation of crushers is spent on the kinetic 
energy of crushing plates and the potential energy of deformation of the springs [2]. 

Based on the general equation of energy consumption for a vibration jaw crusher with three 
vibration masses (2.32), the following equation was derived: 


eee 
a p me + (cx +0,(3x> —Ax,x, + 2X5) + 4c, (Xp x; *)| (2.36) 


n=1 


where m, — respectively n-mass of the crusher, kg; c,, c, — stiffness of the springs connecting the 
masses of the crusher; x,, X, and x, — corresponding displacements of the masses of the crusher, mm; 
c, — the stiffness of the connections of the system connecting the machine to the foundation. 

When determining the total energy of a vibration cone crusher (Fig. 2.8, 6), let's assume that 
the movable cone is an absolutely rigid cylinder rolling on a solid undeformed body in the shape of 
an ellipse. In turn, let’s assume that the crusher body with vibrators is a solid body with mass m.,, 
and moment of inertia /, which is fixed on a fixed base by springs with stiffness c., [1]: 


.9 e 2 
le XO Vin ye || Sol) OO | 2.37) 
2 2 2 eC 


where 7, =1/2m,, xvi +v2)+U,, x O°)/2 — kinetic energy of rolling a movable cone [4]; 
I, = /2m,,v<, — kinetic energy of vibrations of a stationary cone; m,, — the mass of the mov- 
able cone; /_.— central moment of inertia of the cone; x,, y, — Cartesian coordinates of the center 
of the cone; @ — angle of rotation of the movable cone around its own axis; m,, — total mass 
of the stationary cone; v,. — speed of vertical displacements of a fixed cone; c,.— material stiff- 
ness coefficient; c, — the stiffness coefficient of the elastic system connecting the stationary cone 
and the crusher body. 

sandstone with a strength limit of up to 20 MPa was taken as the model of the material of 
destruction. Based on the equations of motion of a vibration jaw crusher (2.31), the speeds of 
movement of the vibration masses of the crusher shown in Fig. 2.9. 

Regarding Fig. 2.9, it should be noted that these speeds of mass fluctuations are valid at the 
frequency of the disturbing force f=19 Hz and the total mass of the laboratory unit m,,,=334 kg. 

Using the first part of equation (4) and the equations of motion of the masses of the vibration 
jaw crusher (2.31), graphs of the consumption of kinetic energy were built, which is spent on 
moving the vibration masses of the vibration crusher Fig. 2.10. 
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O Fig. 2.9 Speeds of movement of masses of vibration jaw crusher: 
1 — first mass; 2 — second mass; 3 — third mass 





O Fig. 2.10 Graph of changes in the kinetic energy of vibration masses of a vibration 
jaw crusher: a — first mass; b — second mass; c — third mass 


From the graphs in Fig. 2.10, it is clear that the total kinetic energy of the motion of the 
masses is within 120 J. However, equations (2.36) used to plot the graphs in Fig. 2.10 do not 
take into account the loss of dry and viscous friction force. 

To construct graphs of potential energy consumption for the crushing process, let’s use equa- 


tion (2.31). 
The movement of vibration masses of the crusher, determined numerically based on the equa- 


tion of motion, is shown in Fig. 2.11. 
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O Fig. 2.11 Change in the amplitude of vibration masses of a vibration 
jaw crusher: 1 — first mass; 2 — second mass; 3 — third mass 


From the graph in Fig. 2.11 it can be seen that the system setting is correct, that is, the sec- 
ond and third masses oscillate in antiphase. The total amplitude of displacement of the second and 
third vibration masses is in the range of 8-10 mm. The minimum required value of the amplitude of 
vibration of the crushing plates of a vibration crusher, at which destruction occurs, is within 3-9 mm. 

The graph of the total costs of potential energy for the destruction process in the crushing 
chamber of a vibration jaw crusher is shown in Fig. 2.11. This graph was constructed by numerous 
methods using equations (2.31) and (2.36). 

The graphs in Fig. 2.10 and Fig. 2.12 indicate an insignificant expenditure of energy for the de- 
struction process. It should be noted here that the system is tuned for resonant operation. That is, in 
the graphs in Fig. 2.10 and Fig. 2.12 does not reflect the energy consumption for the system to en- 
ter the resonance mode and maintain its operation in this mode. In addition, it should be borne in mind 
that, in order to simplify, the equation of motion (2.31) does not take into account energy dissipation. 

To assess the energy of destruction of a body containing a crack, graphs of the instantaneous 
force arising when the crushing plates of the second and third masses of the vibration crusher col- 
lide with the material in Fig. 2.13. To construct this graph, equations (2.31) and (2.36) were used. 

The given graph (Fig. 2.13) testifies to the adequacy of the calculation scheme to the accepted 
premises. So, the first mass moves in antiphase with the second mass. The third mass, accumulat- 
ing energy, is an auxiliary source of energy transfer to the material crushing process. 








O Fig. 2.12 Vibration jaw crusher potential energy graph 
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O Fig. 2.13 The value of the instantaneous force when the vibration 
masses of the vibration jaw crusher move: 1 — first mass; 
2 — second mass; 3 — third mass 


In equation (2.36), to determine the energy consumption of a vibration jaw crusher, the frac- 
ture resistance of the material is considered as the added mass. That is, the indicator of the stiff- 
ness of the material is the parameter that includes energy consumption at the micro level. In equa- 
tion (2.37) for a vibration cone crusher, this parameter is directly the stiffness coefficient of the 
material. Thus, the obtained analytical dependences (2.36), (2.37) reveal the physical essence of 
the grinding process, take into account the conditions of interaction between the crusher and the 
material to be crushed. The practical application of the obtained dependences requires the nume- 
rical values of the characteristics of a particular material, which is the subject of further research. 


2.4 DETERMINING THE OPTIMAL CHARACTERISTICS OF MATERIAL GRINDING PROCESSES 


otudy of the influence of the stiffness index on the amplitude. When designing vibration ma- 
chines, two modes are recognized as expedient: resonant (near resonance), with a frequency 
ratio w/w, = 0.85...0.95; overresonant (only for centrifugal vibration exciters) with a significant 
distance from resonance, when w/o, >7 [12, 13]. 

The advantage of the resonant mode is to achieve a given amplitude with less forced force. 
The resonant mode is characterized by the stability of the amplitude in case of random changes in 
resistance and frequencies of natural and forced oscillations. 

The operation of a vibration crusher in one of the listed modes is achieved by choosing elas- 
tic links that provide a vibration frequency, which may differ from the frequency of the forcing 
vibrations. The stiffness of the system depends on two main parameters of mass and vibration 
frequency. By manipulating these parameters, it is possible to determine the optimal stiffness value 
for the corresponding mode. 

To begin with, for given values of the masses of the system, let's change the value of the 
frequency. It should be noted that the stiffness of the system should help to ensure an effective 
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grinding regime for the material. To determine the stiffness of the system, let’s use the following 
relationship [5]: 


¢, = («,n)) (2.38) 





where w — natural angular frequency of oscillations (resonant frequency); n — coefficient of stiff- 
ness regulation in frequency. Substituting the corresponding values of the stiffness calculated by 
the formula (2.38) into the equation of the amplitudes of the displacement of masses (2.37), 
let’s obtain the graphs of the dependence of the amplitude on the coefficient n (Fig. 2.14) for the 
corresponding frequency modes [6, 8]. 





x,mm 














O Fig. 2.14 Graph of the dependence of the amplitude on the coefficient n 
at the following frequency modes: a — @,=104.67 rad/s; b — m-@,=314 rad/s 


The following conclusion can be drawn from the graphs in Fig. 2.14. The rational value of the 
coefficient nis within 0.6 < n <1. With an increase in the value of the coefficient (1< n< 2.5), the 
effective range of the crusher narrows and shifts towards higher frequencies. Graphs Fig. 2.14 
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show that the stiffness of the elastic bonds of the crusher ensures its efficient operation in the 
resonance and near-resonance modes, namely, between the second and third resonances. With 
the value of the coefficient 0.6 < n <1, the amplitudes of the masses have the greatest values and 
decrease as n approaches 1.4. With an increase in n>1, the change in the amplitude of the second 
mass is insignificant, and for n>1.4, the change in the amplitudes of the first and third masses 
is insignificant. The amplitude of the first mass reaches zero in the vicinity of the point n=2.92. 

Study of the effect of mass on the amplitude. Rationally selected ratios of vibration masses 
provide the required operating mode of the crusher (in-phase, anti-phase) and provide the required 
vibration amplitudes of these masses [13]. 

When choosing the ratio of mass m, and m,, the value of the mass m, is taken constructively, 
the value of the stiffness of the systems c,, c. is taken taking into account formulas (2.31), (2.38). 
To introduce the mass ratio coefficient depending on the equations of motion (2.37), let’s use the 
following logical transformations [7]: 


Nog = My +My; Koy = Mp / My. (2.39) 
Based on this: 
my = (Nogkog)/(1 gq); Mg = Npg/(1-+Kog). (2.40) 


On the basis of the above material, the graphs of the dependence of the amplitude on the 
mass ratio for different frequency spectra were plotted in Fig. 2.15. Analyzing the graphs in 
Fig. 2.15 the following can be noted. The effective value of the coefficient k,, for the resonant 
angular frequency w, = 104.667 rad/s is within the limits 0<k,, < 2.5, for w, =314 rad/s — 
1.5<k,, $2.45 [14-18]. 

Analyzing each graph separately, Fig. 2.15, the following conclusions were made. At the re- 
sonant frequency @, = 104.667 rad/s, the maximum value of the amplitude of the third mass is 
reached in the vicinity of the point k,, =1.3. And the maximum values of the amplitudes of the 
first and second masses are in the vicinity of the point k,, =0.79. For the resonant frequency 
@, = 9314 rad/s, the maximum values of the amplitudes of all three masses are achieved under the 
condition of the point k,, lying in the vicinity k,, = 1.9. 

Let’s use the same approach to find the optimal ratio of masses m, and m, — kz,. 

Construction of the amplitude-frequency characteristics of the vibration jaw crusher. Knowing 
all the unknown parameters of the equations of motion, graphs of the dependence of the amplitude 
of displacement of the crushing plates on the frequency of change in the disturbing force were 
built (Fig. 2.16, 2.17) [3, 5]. 

Based on the graphs (Fig. 2.16, 2.17), the following conclusions were made: 

— the effective mode of operation of the crusher is realized in the zone of the second resonance 
and in the range between the second and third resonances; 
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— with an increase in the calculated value of the resonant frequency of wheel vibrations ow, 
(2.38), the range of effective operation of a three-mass vibration jaw crusher shifts towards 
an increase in frequency; 

— the effective ranges of the crusher operation at the system stiffness values calculated at the 
corresponding resonant frequencies of natural vibrations @, (2.38) will be as follows: 


100 rad/s<w@ < 180 rad/s; w, = 104.667 rad/s: 
97 rad/s<@ <121rad/s; 186 rad/s<q@ < 234 rad/s; 


, = 314 rad/s; 291 rad/s< w< 322 rad/s; 365 rad/s<w < 400 rad/s. 














O Fig. 2.15 Graphs of the dependences of the amplitude change on the mass 
ratio K,, at different frequency ranges: a — @,=104.67 rad/s; b — @,=314 rad/s 


To ensure the operation of the crusher at the values of the stiffness of the system, calculated 
from the resonant frequency of natural vibrations of the w, = 104.67 rad/s (2.38), in the lower 
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frequency ranges, it is necessary to increase the masses of the vibration parts of the machine. 
The increased masses should not exceed the permissible limits, exceeding which will lead to the 
economic inexpediency of this machine design [3, 5]. 
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O Fig. 2.16 Amplitude-frequency characteristic of a vibration jaw crusher 
at @,=104.67 rad/s: a — crusher operation without material: 
6 — crusher operation with material 


The study of the effect of frequency on the operation of the crusher with changing parameters 
of the disturbing force and stiffness is presented in the form of vibrograms in Fig. 2.18, 2.19. 

The range at which the experimental studies of the operation of the three-mass vibration jaw 
crusher were carried out covers three modes: pre-resonant, near-resonant and over-resonant. 
As the work of the crusher approaches resonance, the difference between the values of the theo- 
retically calculated and experimentally determined amplitudes slightly increases, which is explained 
by not taking into account the theoretical calculations of dissipative resistances. When the crusher 
is operating at low frequencies up to f=/ Hz, that is, before the first resonance, there is a fluctu- 
ation of the masses in one phase. It follows from this that this mode is not effective. 
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O Fig. 2.17 Amplitude-frequency characteristic of a vibration jaw crusher 
at @,=314 rad/s: a — crusher operation without material; 
b — crusher operation with material 


When the crusher is operating in the vicinity of a frequency of 12 Hz, there is almost no crush- 
ing of the material. This can be explained by the fact that the vibrator does not create forces at 
a given frequency that could effectively destroy the material. 

In the frequency range 14.5 Hz<f<26 Hz, the material is efficiently comminuted. The efficien- 
cy of crushing is somewhat reduced when the crusher is operating in the vicinity of a frequency 
of 22 Hz. This can be explained by a decrease in the amplitude of oscillation of the second and third 
masses and a decrease in vertical oscillations of the crusher body, which increases the residence 
time of the material in the crushed chamber. 

When studying the influence of the stiffness of the elastic systems of the crusher on the 
crushing process, the following features were established. As the stiffness of the system increa- 
ses, the effective range of the crusher shifts towards an increase in the angular frequency. This is 
due to the shift of the resonance point towards an increase in the angular frequency of the crusher 
vibrations. This property can be used to create a crusher control system. 
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O Fig. 2.18 Vibrograms of vibrations of the masses of the crusher 
when working with material (F=4 111 N, c,=219 052 N/m): 
a — displacement of the second mass; b — displacement of the third mass 





























O Fig. 2.19 Vibrograms of vibrations of the crusher masses when 
working with material (F=4 111 N, c,=234 438 N/m): a — displacement 
of the second mass; b — displacement of the third mass 


However, too much increase or decrease in stiffness adversely affects the performance of 
the machine. So, with a stiffness N/m, in the frequency range, the second and third masses move 
in phase, and the first in antiphase. Thus, the operation of the crusher at low frequencies is not 
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efficient at a given stiffness. At very low stiffness, the efficiency of energy transfer from the vibrator 
to the third mass deteriorates, as evidenced by experimental data. The vibration amplitudes of the third 
mass, provided that N/m are insignificant. The features of the impact of stiffness are as follows. As 
stiffness increases, the crusher’s effective operating range shifts towards zero. If the stiffness is in- 
creased too much, the three-mass system turns into a two-mass one. With a decrease in stiffness, the 
range of effective operation of the crusher narrows, and in conditions of too large a decrease in stiff- 
ness, the system leaves the state of equilibrium, which has a negative effect on the crushing process. 


2.9 DISCUSSION OF RESEARCH RESULTS 


An efficient crusher operation can be achieved by operating it in the vicinity of resonance, 
i.e. at a high value of the coefficient of resonant amplification of oscillations. However, this mode 
is sensitive to changes in external influences on the workflow. It is possible to reduce external 
influences by automatically adjusting the operating mode of the machine. 

The control of the vibration machine is to ensure the optimal operating mode, in which the 
specified displacements, speeds and accelerations are achieved [6-8]. The set values of displace- 
ments, speeds and accelerations of a vibration jaw crusher can be achieved both in pre-resonant 
and over-resonant operating modes. 

The creation of a control system for the amplitude-frequency characteristic and taking into ac- 
count the dependences of the moment of resistance on the frequency of the driving force underlies 
the method of regulating the working process of the crusher. This task can be achieved by correctly 
calculating the parameters of the three-mass vibration jaw crusher. 

The second way to regulate the operating mode of a three-mass vibration crusher is to adap- 
tively change the stiffness of the system. The third way to regulate the operating modes of the 
crusher is to change the vibration frequency. 

On the basis of the above studies, algorithms were formulated for calculating the main para- 
meters and controlling the working process of a vibration three-mass jaw crusher (Fig. 2.20, 2.21), 
respectively. 

A physical model and a mathematical model of a vibration jaw crusher have been developed, 
taking into account the elastic — inertial properties, power and geometric characteristics of the 
«jaw — fractional material» system. 

One of the main research tools were the amplitude — frequency characteristics of the system 
movement. The studies have established the regularities of the change in the elasticity coefficients de- 
pending on the frequency, the coefficients that determine the rational ratios of the vibration masses. 

The optimal values of the amplitude of the oscillations of the masses have the greatest values, 
provided that the coefficient of regulation of elasticity in frequency n is in the range of 0.6<n<1. 
Within the frequency range of 14.59 Hz<f<26 Hz, material is effectively crushed. The numerical 
values of the coefficient k,, are established, which determine the rational ratios of vibration masses. 
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The studies carried out and their results indicate that the main stresses causing the destruction 
of the material are shear stress. With an increase in the size of the major semiaxis of the crack 
ellipse, the total stresses decrease, which indicates the destruction of the brittle body. However, 
the stresses op acting along the radius slightly increase. The graphs (Fig. 2.4) reflect the picture 
of the stress distribution in a spherical body relative to the point of application of the force along 
the angle ~. The graph (Fig. 2.5) shows the change in normal stresses depending on the distance r. 
Analysis of the graph (Fig. 2.6) makes it possible to conclude that at r—> O, the stresses increase 
indefinitely. As a result of this, the stresses and strains in the cracked body must be determined in 
the vicinity of the center of the body. In addition, it was found that the complete elliptical integral of 
the second kind has an insignificant effect on the change in the stress pattern in the body. It should be 
noted that when solving the problem using the Legendre functions, 15 constants appear in the stress 
equations, which together have a significant effect on the stress pattern in the body. These con- 
stants were a numerous method based on the equations of the potential energy of a deformed body. 

Equations (2.21), (2.23) and (2.30) fully reflect the picture of energy consumption for material 
destruction. Such studies are planned as a continuation of the topic under consideration through 
experimental studies and consideration of several members of the series and optimization of the 
parameters of crushing machines. The proposed approach to studying the energy characteristics 
of material destruction in the crusher chamber with a guaranteed zone of stability of parameters 
in the resonance zone can be used for other processes. Such processes include the destruction of 
materials during dynamic cutting of soils, processing of materials in mills and sorting of materials, 
which are gaining widespread use in various European countries. 


CONCLUSIONS TO SECTION 2 


1. A physical model and a mathematical model of a vibration jaw crusher have been developed. 

2. The stress-strain state of the material in the process of its destruction for a continuous 
medium and surrounded by an existing crack has been investigated. 

3. The efficiency of the crusher’s working process is proved by the practical achievement of its work 
in the vicinity of resonance. That is, at a high value of the resonant amplification factor of the oscillations. 

4. lt has been determined that the optimal values of the amplitude of the oscillations of the mas- 
ses have the greatest values, provided that the coefficient of regulation of elasticity in frequency n is 
within 0.6<n<1. Within the frequency range of 14.5 H7<f<26 Hz, material is effectively crushed. 

9. The proposed approach to study the energy characteristics of material destruction in the cru- 
sher chamber with the provision of a guaranteed zone of stability of parameters in the resonance zone. 

6. The studies carried out and their results indicate that the main stresses that lead to the 
destruction of the material are shear stress. 

7. Algorithms for calculating the parameters of a vibration three-mass crusher have been 
formulated. 
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ABSTRACT 


The studies carried out revealed the influence of the vibration frequency and the angle of in- 
clination on the efficiency of sorting. It was found that with an increase in the vibration frequency 
and vibration amplitude of the box, the separation of the material and the sorting speed increase. 
The proposed model made it possible to determine the transport speed in the case when the sieve 
screen carries out a circulating movement, which is the application of two independent vibrations 
with different amplitudes and frequencies. The number of throws of the material to be sorted is 
influenced by the transport time of this material over the surface of the sieve and the amount of 
vibration excitation to this material. The dependence of the used power on the design and techno- 
logical parameters of the vibration screen was revealed. An algorithm and method for calculating 
the main parameters of a vibration screen were developed. 


KEYWORDS 


Vibration screen, material, sorting, grain composition, physical, mathematical model, equation, re- 
sonance, parameters, experimental setup, sieve, particles, amplitude, vibration frequency, productivity. 


3.1 SELECTION AND JUSTIFICATION OF PHYSICAL AND MATHEMATICAL MODELS OF TECHNICAL 
SYSTEMS FOR MATERIALS SORTING PROCESSES 


As a technical system for the processes of sorting materials, let's take a vibration screen, 
which is widely used in industry [1—6]. Let’s imagine the sorting process as some ordered process 
of movement of a large number of different parts in a layer on a sieve (Fig. 3.1). Each share with 
size d characterizes its size, which does not change during movement, and is a value that can have 
a discrete set of values: d,, d,,..., d,. Since the number of particles in the layer is large, then d can 
be assigned arbitrary values in the interval of a given layer of fractions, and thus its value can be 
considered continuous. 

H.. ,,, — the height of the superlattice and sublattice product; S — sieve area; A/ — sieve 
section; D, d— average diameter of superlattice and sublattice products; A — thickness of the sieve 
wires; V,,. V,,, — volume of the superlattice and sublattice product. 

Based on this assumption, the particle system can be replaced with a continuous model, and 
the productivity P will be determined by the formula: 
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P=V/t = bhv = bhl/t = sh/t, (3.1) 
where V — volume of the material layer on the sieve, m®; t — time of transporting the material over 


the sieve, s; b — sieve width, m; v — sieve length, m; S — sieve area, m*; h — average layer height 
on a sieve, m; v— transport speed, m/s. 
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O Fig. 3.1 Layout of material on the screen: /, b — length and width of the sieve 


The volume can be represented as the sum of the volumes of the superlattice and sublat- 
tice materials: 


V = Veep 7 Viub =; Sr ee 7 Sap! lub: (3.2) 
If to divide the sieve along its length at the level of the section A/ (Fig. 3.1), then: 
Al = |/n = const, (3.3) 


where n — an arbitrary number. Then the volume V, will be distributed in equal portions V, by Al: 


Vo = AlN Psu) = Const, (3.4) 


c “sup Sup 


where c — the number of the segment Al. 
Using the main provisions of the theory of probability [3], the process of separation into frac- 
tions can be represented as the probability p of the passage of grains d through cells D: 
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p =(D-d¥/(0+8F, (3.5) 


where /=1,2,3...m — the serial number of the fraction, counted according to the characteristic 
of the size of the undersize product in the direction of growth d; d, — the average grain size of the 
i-th fraction, corresponding to the characteristics of the size of half of its output; A — sieve wire 
thickness. The volume V, distributed in this way along the length / forms a layer with the dimen- 
sions of a parallelepiped, averaging the inhomogeneity of its constituent grains. For a product with 
a volume of V,, this means that all the sieve cells through which it is obtained, for the case of close 
alignment, a rectangle is formed, the area of which is S,. 
According to (3.4): 


V.=A¥ hb, (3.6) 
k=1 


where h, and 6, is the average value of the height and width of the layers V,, depending on the 
probability p. 

Consequently, the grain composition of each portion volume differs from the grain composition 
of similar portions of the undersize product and is a value: 


v= SV, (3.7) 


where V,,— part of the volume of the /-th fraction V, sorted at the section Al. 
Under condition (3.4), the height of the /-th layer: 


hig = Ope | Dig (3.8) 


where P,,,,— probability of sieving the largest grains belonging to the composition V.,. 

Considering that there is a probability that parts of the material do not pass the parts of the 
material |...) U—/=h4.)) and condition that 1/p,,,,,=/,/A/, it is possible to obtain an expression for 
the volume of material passage, which will be proportional to: 


V,,=V, (y : PMs / Su -v,.)} 3.9) 


Expression (3.9) is dependence for determining the rational load of the sieve during sorting 
based on a continuous model obtained on the probability of sorting grains. 

For practical calculations, the considered model can be replaced by an adequate discrete model. 
For this purpose, let’s express 5,, taking into account S.,,, through the area occupied by the frac- 


tion of grains, which are located near D in size. Let's call this fraction «heavy», that is, boundary 
in the general fractional composition. Then there are the following relations: 
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=e ~~ 5,Ah,,,AC,,, 
5, =5,Ah,AC, (3.10) 
S55 a 5,Ah sub AC up 


where Ah — ratio of heights h,/h 
the yield of «heavy» grains). 
From the relation (3.10) it is possible to determine the area of «heavy» grains: 


; AC — ratio of C,,,/C, in percentage of the initial product (C, is 


sub} sup 


sup 


5 = (Say + S5y,)/(C,Ah +6,,Ah, + Cs). (3.11) 

Using the probability of passing (3.5) and assuming that with each throw V, of the layer V, 
onto the sieve, «heavy» grains are sown along its width n = b, /d,, and along the length of the sieve 
n=‘1/d,,, let's obtain the volume of limiting grains sorted in one throw of the layer on the sieve: 

V, =(nd,,b, -!)/6. (3.12) 

Then the capacity for sorting heavy grains will be: 

P =V,/T =0.08d, Clo. (3.13) 

Accordingly, the overall optimal power performance: 


P=100P,/C,. (3.14) 


Taking into account (3.11), let’s finally get the formula for determining the performance: 
p= (5,do[1 -d/ of) (ConA Mi +C,,Ah,, +6,)-(1-2)], (3.15) 


where 5, — area of the open section of the sieve; A — voidness of the material; Ah, = Area / Amin 
Ah, = uae 


3.2 INVESTIGATION OF THE DYNAMICS OF THE RESONANT SCREEN, ANALYSIS 
AND ASSESSMENT OF ITS PARAMETERS 


The proposed model [6] allows to indicate the speed of transport in the case when the sieve 
shaker carries out a circulating movement, which is the application of two independent vibra- 
tions with different amplitudes and frequencies. The design scheme for determining the rational 
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parameters of screens (Fig. 3.2) displays a single particle on the surface of the sieve, which 
vibrates according to the most general laws [7, 8]. 





O Fig. 3.2 Calculation diagram of the process of particle motion over the sieve 


Differential equations for the relative motion of a particle over the sieve surface in a movable 
coordinate system associated with it: 


V,, = —geosa + X,w sin(w,t + @,), (3.16) 
iy’ =V,, (3.17) 
V,. =gsina. + X,@° cos(w,t +@, ), (3.18) 
V=V,, (3.19) 


where V and V — projections of the visualization speed of the relative motion of the particle, and 
the last terms in (3.16) and (3.17) correspond to the inertial force of the translational motion: 


F, =mo,X, sin(w,t +, ), (3.20) 
F,=maX, sin(w,t +@, ). (3.21) 
Equations (3.16) and (3.18) do not take into account the forces of air resistance to the motion 
of the particle. A particle lying on the sieve surface passes into the station field predetermined by 


the field (3.16)—(3.19), provided that: 


X,a; sin(@,t + @,) > gcosa. (3.22) 
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Upon reaching the surface of the sieve in free flight (y’ = 0), the particle strikes the surface 
of the sieve. The change in speed from impact can be described by the ratios of inelastic impact. 


Vi" =-RV,, (3.23) 


y 


Ve =V TRS). (3.24) 
where FR — coefficient of speed recovery from the impact, f — coefficient of friction of the particle 
with the sieve surface, the subscripts «—» and «+» correspond to the states immediately before 
and after the impact. On the one hand, there is a significant uncertainty in the experimental de- 
termination of these coefficients, and on the other hand, in real conditions the particle strikes not 
against the sieve surface, but against the layer of particles on it, after which its relative speed is 
practically zero. Therefore, with an accuracy acceptable for practical calculations, it can be con- 
sidered V,, = V,’ = 0 with each impact. The coefficient of friction f will be defined as the coefficient 
of internal friction or the angle of natural slope of the bulk material. After a particle hits the sur- 
face of a sieve with layers of particles on it, several options for its further behavior are possible: 

— if X a sin(a,t + °,) > gcosa, then the particle is detached from the surface and conti- 
nues to move above the surface in accordance with equations (3.16)—-(3.19): 

— if at the moment of attachment of the particle to the surface and for some next period 
of time Xa sin(w,t +@,) < gcosa, then the particle remains on the surface until the sign of 
the roughness changes. Its movement along the surface of the sieve during this period of time is 
determined by the following conditions: 
if 

(gcose. —X,@, sin(o,t + 9, )) 





f >|gsino. + X,a, sin(w,t + @,)), (3.25) 


then 


V. =O, (3.26) 


xX 


and the particle is stationary on the surface of the sieve; 


f 
(gcoso - X,«; sin(w,t + @,))|F <|gsino. + X,«; sin(w,t + @, ) (3.27) 
then 
V,, = gsino. + X,o,° sin(w,t +9, ) f(gcosar — Xo," sin(,t + ,)), (3.29) 
ce (3,29) 
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Consequently, the particle moves along the surface under the action of the projection of the 
gravity force, the variable friction force due to the variable pressing of the particle to the surface, 
and the longitudinal exchange force of inertia of the translational motion. The system of equations 
(3.16)—(3.19) with nonlinear conditions (3.25)—(3.29) can be solved only by a numerous method. 

With an independent disturbance of the vibrations of the sieve in the longitudinal and trans- 
verse directions, the movement of the particle above and along the surface is rather complicated. 
In some practically important cases, which are necessary for evaluating certain parameters of 
vibration sorting, the system of equations of motion can be significantly simplified. Experiments with 
the model described above allow one to determine the characteristics of the motion of a particle 
over a horizontal surface, performing vertical oscillations, necessary for considering periodic sorting. 
Of these characteristics, the vibration amplitude and frequency of contacts are important, the 
vibration amplitude determines the increase in the layer height due to the effect of its «swelling». 
The frequency of contact of particles with the surface determines the conditions for their passage 
through the openings of the sieve. Calculations show that, regardless of the initial phase, a steady 
cycle of a share is formed rather quickly. From this, the amplitude of the particle tosses over the 
surface and the period of contacts with it are determined. 

Fig. 3.3 shows the effect of the angular frequency and the amplitude of surface vibrations on 
the amplitude of particle motion. 





x,mm 





O Fig. 3.3 Changing the vibration amplitude of x material particles 
from frequency w: 1 —x=4 mm; 2 —-x=6 mm; 3 -—x=8 mm; 4-x=10 mm 


With an amplitude of 5 mm, the separation of a particle from the surface of the sieve begins 
at a circular frequency of 45 rad/s, and for an amplitude of 2.5 mm — generally at a frequency 
of 63 rad/s. The effect of circular frequency and surface amplitude on the frequency of particle- 
surface contacts is shown in Fig. 3.4. 

Each amplitude has its own rotation frequency corresponding to the maximum frequency of 
contacts, that is, the fastest passage of particles through the holes. At a frequency of 50 rad/s, 
the effect of the surface amplitude on the contact frequency is not significant. 
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O Fig. 3.4 Influence of the vibration frequency of the screen on the 
frequency of contacts of particles with its surface at different amplitudes: 
1-10 mm; 2 — 2-8 mm; 3 - 38-6 mm 


The sorting process is influenced by the likelihood of the passage of the grains through the 
openings of the sieve. This probability depends on the size of the particles, the dynamic parameters 
of the vibrations of the screen, the design of the screen, and the shape of the holes in the screen. 
The average speed of passage of particles through the sieve is determined by the frequency of 
contacts of the particles with the surface and the ratio of the sizes of the particle and the hole. 
The frequency of particle-sieve contact is influenced by the amplitude and frequency of vibrations of 
the screen surface. An estimate of the probability of a particle passing through the sieve opening in 
one collision can be carried out on the basis of the design scheme shown in (Fig. 3.9). 








O Fig. 3.5 Scheme of passage of particles through the sieve openings 
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Taking into account the direction of impact of particles on the surface perpendicular to it, it can 
be assumed that only particles which centers are inside the square (/, —d) x (/, — d) will pass into 
the hole. In general, particles can reach the surface anywhere in the /x/ square. The probability of 
particles passing through the holes is determined by the relationship: 


p=nAtK,(1-d/I,) , (3.30) 


where n=n(X,,@) — the number of contacts of the particle with the sieve surface per time unit; 
© — coefficient of the open section of the sieve; d — particle diameter; /, — hole size; At — time 
of one transition in the matrix of transition probabilities g of the cellular model of the periodic 
sieving process. 

The angle of inclination of the sieving surface « determines the probability of particles pen- 
etrating through its holes p,, which is determined by the angle at which the particles attack the 
sieve surface. This angle depends on the parameters of the vibrations of the sieve, the trajectory 
of the screening surface and the angle of its inclination to the horizon; it can vary within wide limits. 

For experimental research, a model of a vibration screen has been developed and manufac- 
tured (Fig. 3.6) [9-12]. 
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O Fig. 3.6 Experimental device: a — diagram; b — general view; 1 — engine; 2 — feeder; 
3 — mechanism for regulating the supply of material; 4-6, 8 — sensors; 7 — measuring complex; 
9 — oscillator; 10, 13 — sieves; 12 — box; 14 — reactive mass; 11, 15 — elastic elements 


The model, in terms of its design and technological parameters, is analogous to a serial vi- 
bration screen, and the ratio of the length to width of the screen was taken according to the 
standard ratios of 2:1. 

In the course of the experiments, the search for the most efficient sorting modes was carried out. 

With a certain number of throws n,_.. and sizes of boundary grains, experimental 1 and cal- 
culated 2 graphs were built (Fig. 3.7). 
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O Fig. 3.7 Dependence of the probability of screening 
grains (the number of throws) on their size: 1 — calculated; 2 — experimental 


The position of curve 2 corresponds to the position of curve ‘1 in terms of the proportionality 
coefficient K,=0.99. The probability of sorting the boundary grains of the crushed mixture accord- 
ing to the formula (3.5) is determined by the ratio: 


I-Kd,\K 
_ (Kite) Kp ie) 9 (3.18) 
(/+d) 


sor 


In accordance with this correction, let’s perform calculations of the optimal productivity and 
the corresponding screening efficiency for experimental screening units. At the same time, let's 
neglect a similar correction for light grains, taking into account the insignificant size dl. 

Let's find the average value of the grain size composition of the material in each experiment and 
the empirical dispersion S° for each fraction. As an estimate of variances, let’s take the weighted 
average number of individual empirical variances. 

Let’s compare the obtained permissible measurement accuracy, having previously determined 
the nature of the variance discrepancies by the ratio. 

The critical value of these ratios for the degrees of freedom K,=49 and K,=13, F=2.3. This 
confirms the random nature of the variance discrepancies and the absence of a significant differ- 
ence in the measurement accuracy in individual experiments. Hence, with a reliability P=0.95 and 
an accuracy S’ for the actual grain size composition of the material, one can take the average, 
compiled according to the data of all experiments. As the result shows, the initial characteristic of 
the size of the material practically did not change during the experiments. Using the experimental 
and calculated points, let’s build the graphs of the dependence (Fig. 3.8). 

Analysis of experimental data on the number of throws n,_.. revealed the following. The va- 
lue 1, ,9 iS influenced by two main parameters: the time of transporting the material over the 
sieve t and the frequency of forced vibrations of the vibrator. The time of transporting the 
material ¢ increases the screening efficiency E even in the case when the limiting grain size d,,, is 
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large, i.e. d,,/D>0.9. On the other hand, an increase in the parameter t reduces its vibration 
transport speed v, which reduces the total amount of material passing through the screen per 
unit time. For example, at -=5 s, P=26...43 kg/min; for £=6 s, P=21...38 kg/min; for t=7 s, 
P=16...24 kg/min; for t=12 s, P=Y...21 kg/min; for -=16 s, P=3...14 kg/min; for t=26 s, 
P=4...9 kg/min. The graphical meaning is shown in Fig. 3.9. 
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O Fig. 3.8 Dependence of rational sorting productivity on the time 
of material transportation over the sieve: 1 — calculated; 2 — experimental 
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O Fig. 3.9 Dependence of the maximum sorting efficiency 
on the number of throws of the material layer on the sieve: 
1 — calculated curve; 2 — experimental curve 


As can be seen from the above examples, with an increase in the residence time of the material 
on the sieve ¢,, an inversely proportional decrease in the productivity of the vibration screen is 
observed. The size of the sorted material has a significant effect on the value of productivity. In this 
case, the productivity P can vary up to 42 %. For example, in the case of transportation during 
the time of material t=12 s at d,,=8.4 mm, at d,, /D =0.907 and S,=9.64 %, the productivity 
of the screen P=10.2 kg/min, and at d,=8.32 at d,, /D=0.924 and 5,=10.8 % screening 
capacity is 12 kg/min. 


o3 


DYNAMIC PROCESSES IN TECHNOLOGICAL TECHNICAL SYSTEMS 


Analysis of the experimental data (Fig. 3.9) on the number of material throws n,_.. over the 
sorted material and their influence on the main technological parameters showed that the dis- 
crepancy is 9.3-8.1 %. 


3.3 DISCUSSION OF RESEARCH RESULTS 


The studies carried out revealed the influence of the vibration frequency and the angle of incli- 
nation on the efficiency of sorting. So, with an increase in the frequency of vibration and the ampli- 
tude of vibration of the box, the separation of the material and the speed of sorting increase. The 
particle speed is of greater importance at the frequency of the box vibrations, w = 45...85 rad/s. 
When changing the angle of inclination from o = 2°... 204° it was obtained that with an increase 
in the angle of inclination of the sieve from 14° to 20°, the concentration, speed and power are 
more stable. The dependence of the used power on the design and technological parameters of 
the vibration screen is revealed. The number of throws to be sorted is influenced by the transport 
time t of this material over the surface of the sieve and the amount of vibration excitation to this 
material. These parameters are due to the magnitude of the forced vibrations of the vibration 
system, with an increase in both of these parameters [13-16]. 

An increase in the number of throws of material n,_... increases the residence time of the 
material on the surface of the sieve, reduces the speed of passage of the material and, in general, 
reduces the productivity of the screen. For example, for n,_,,=49, P=26...43 kg/min with an 
increase n,__.,. t0 6/7, the productivity of the screen P decreases to 16...24 kg/min, for n,_, =211 
it decreases to 3...14 kg/min, and with a value of n,_.. =36/, the productivity of the screen 
is only about 8.1 kg/min. 

The efficiency of the sorting process E is also affected by n,_.. as follows: with the number 
of throws n,__.. indicated in the previous paragraph, the efficiency of screening E is respectively: 
87.9 %, 89 %, 91.5 %, 94.7 % (Fig. 3.9, curve 2). The increase in the efficiency of screening F 
in accordance with the increase n, _.., is explained by the increase in the time of vibration exposure to 
the medium being separated. This is confirmed by the close placement of curves 1 and 2 (Fig. 3.9). 
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CONCLUSIONS TO SECTION 3 


1. The physical and mathematical model of the screen has been substantiated, on the basis 
of which the equations of material movement over the sieve are compiled. 

2. Distribution of the vibration amplitude of the material particles from the vibration frequency 
of the sieve has been determined. 

3. An experimental model of a vibration screen has been developed and experimental studies 
of its operation have been carried out. 
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4. Dependence of the used power on the design and technological parameters of the vibration 


screen has been revealed. 


9. The influence of the vibration frequency and the angle of inclination on the sorting efficien- 


cy has been determined. So, with an increase in the frequency of vibrations and the amplitude 
of vibrations of the box, the separation of the material increases and the particle speed is of great- 
er importance at a frequency of vibrations of the box of 45...85 rad/s. 
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RESEARCH OF TECHNICAL SYSTEMS OF PROCESSES 
OF MIXING MATERIALS 


|. Nazarenko, |. Bernyk, 0. Dedov, 
|. Rogovskii, M. Ruchynskyi, |. Pereginets, L. Titova 
ABSTRACT 


The process of a vibration concrete mixer has been investigated. Well-coordinated equations of 
motion of the system «mixer — concrete mixture» and their solutions have established the regulari- 
ties of changing the parameters of the mixer at different angles of inclination. An experimental design 
of vibration-gravitational mixing has been developed and manufactured, and studies of the working 
process of mixing efficiency have been carried out on it. Research carried out at three warehouses 
of concrete mixtures confirmed the assumptions accepted in the work, the assumption and the 
selected physical and mathematical model of the «Mixer — concrete mixture» system. The use of 
a vibration concrete mixer made it possible to reduce the duration of mixing by 2.5—3 times in com- 
parison with conventional free mixing concrete mixers. Provides the preparation of hard concrete 
mixtures and improving the quality of concrete mixtures due to the destruction of defective aggre- 
gates and a more uniform distribution of the binder throughout the volume of the mixed mixture. 

An algorithm and method for calculating the main parameters of a vibration concrete mixer 
have been developed. 


KEYWORDS 


Vibration mixer, material, physical, mathematical model, equation, resonance, parameters, 
experimental setup, vibration frequency, productivity. 


41 JUSTIFICATION OF THE CONSTRUCTIVE DIAGRAM OF THE TECHNICAL SYSTEM 
FOR THE PROCESSES OF MIXING MATERIALS 


The technical systems for the processes of mixing materials include mixers, which are classified 
according to the characteristics: the principle of mixing, the nature of work, the method of instal- 
lation [1, 2]. According to the principle of mixing the components, concrete mixers are divided into 
two classes: with free mixing of materials (gravitational) and with forced [3]. It is most expedient to 
use gravity mixers to prepare more plastic mixtures with lower resistivity when mixing them [4]. For 
the preparation of hard concrete mixes, mixers with forced mixing of materials are most effective. 
At the same time, modern requirements for the quality of concrete mixtures require the use of 
more efficient mixing equipment. Such equipment includes mixers using vibration effect [5]. In such 
a mixer, up to five component movements take place, and mixing is carried out simultaneously 
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by translational and rotational movements of the components of the mixture. A free-mixing vibra- 
tion concrete mixer (Fig. 4.1) includes a cone-cylindrical mixing drum 1, which is mounted on a base 
plate 3 with the help of elastic shock absorbers 2. 





O Fig. 4.1 Vibration concrete mixer 


This base plate 3 is rigidly connected in its central part to a hollow drive shaft 4, mounted 
in a bearing support 0. At the end of the drive shaft 4, a sprocket 6 is mounted, connected to 
an electromechanical drive (not shown). On the bottom of the mixing drum 1, a circular vibration 
exciter 7 is mounted, which is half mounted in the inner cavity of the mixing drum. In this case, the 
vibration exciter of circular vibrations 7 is connected to the drive motor (not shown in the Fig. 4.1) 
using a flexible shaft 8 passing inside the hollow drive shaft 4. Inside the mixing drum, the vibration 
exciter of circular vibrations 7 is closed with a cone-cylinder protective cap. drum. Blades 10 
are located inside the mixing drum 1. The bearing support 5 is mounted on the traverse of the fixing 
and tipping mechanism. 

Free mixing vibration concrete mixer works as follows. 

Mixing drum 1 is installed at an angle of 40—45 degrees to the horizon and is loaded with 
a portioned amount of sand, crushed stone and cement. Pour in the required amount of water. The 
mixer drum drive is included. In this case, the rotating mixing drum using the blades 10 provides 
mechanical mixing of the mixture. This preliminary mixing of the components of the mixture con- 
tinues for 15-20 s. 

Then, with the help of the fixing and tipping mechanism, the mixing drum is set in a verti- 
cal position and the vibration exciter of circular vibrations 7. Under the influence of vibration, 
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the mixing drum performs complex vibration movements that cause variable amplitude-frequency 
deformations in the concrete mixture. The vibration action of the applied simultaneous torsional and 
circular vibrations ensures efficient mixing of the concrete mixture. The duration of vibromechanical 
mixing, vibration treatment and vibration activation is 15-25 s, depending on the concrete mix. 
After the end of vibromechanical mixing, the vibration exciter of vibrations is turned off and, using 
the fixing and overturning mechanism, the mixing drum is transferred to its original position. The 
drum is installed at an angle of 40—45 degrees to the horizon and the process of mixing the mixture 
continues for 10—12 seconds until a homogeneous state is obtained. Then the concrete mixture 
is discharged from the mixing drum. 


4.2 DETERMINATION OF THE PARAMETERS OF THE MATHEMATICAL MODEL 
OF THE VIBRATION MIXER 


The main parameters influencing the choice of the mathematical model of the vibration mixer 
are the nature of the mixing drum vibrations, the coordinates of the center of gravity of the mixing 
drum, the physical and mechanical characteristics of the mixture, the geometric dimensions of the 
mixing drum and the volume of filling it with the concrete mixture. 

When vibration in the vertical direction on the bottom of the mixing drum from the side of 
the processed concrete mixture, support forces will act, which can be represented in the form of 
inertial forces F, and dissipative forces F,.,, i.e. forces of inelastic resistance F,.,. These forces 


can be described by the following dependencies: 


2 
=a ua (4.1) 
Fic = Dor ve (4.2) 
where m,,, — reduced mass of the concrete mixture interacting with the bottom of the mixing 


drum during vertically directed vibrations; 6,,, — reduced coefficient of inelastic resistance of the 
concrete mixture. 

In this case, the value of the reduced mass of the concrete mixture can be determined from 
the dependence, modified as a result of the refinement, given in [6]: 


Mon = tok, (4.3) 
1 


where p — dynamic density of the concrete mixture; F, — the area of the bottom of the mixing drum 
in contact with the concrete mixture; k, — the wavenumber; 
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k=— (4.4) 


a, =,{— (4.9) 
p 


E — dynamic modulus of elastic deformation of the processed concrete mixture; H, — estimated 
thickness of the processed layer, selected from the condition: 


a (4.6) 
4m \ p 


tise eee (4.7) 
4m \ p 4m \ p 


where fh — thickness of the processed layer. 
The coefficient of inelastic resistance of the concrete mixture is determined from the following 
relationship: 


Boy = HAF LOK, A, (4.8) 
where 1, — coefficient of dynamic viscosity of the concrete mixture; 


n, = £,H,JEp. (4.9) 


where f., — coefficient of internal friction of the concrete mixture. 
From the given dependences (4.2) and (4.8), let’s find the specific values of the reduced mass 
and the coefficient of inelastic resistance: 


Myo = [tok mH, (4.10) 
1 


b 


yor 1 


= 1)k,tgk,H.. (4.11) 
The specific values of the reduced mass and the coefficient of inelastic resistance of the con- 


crete mixture interacting with the cylindrical walls of the mixing drum during horizontal vibrations 
are determined from the following dependencies: 








m= — = Pt (4.12) 
me ew) Oke | 
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b 1),k 
bp = et Gtk (4.13) 
wre Oe) ew 








where wt — Poisson's ratio; L. — length of the processed layer in the horizontal direction. 
To determine the nature of the motion of the mixing drum in the vibration operating mode, let’s 
consider the design diagram of the dynamic system «mixing drum — processed medium» (Fig. 4.2). 


















































O Fig. 4.2 Design diagram of the dynamic system «Mixing drum — processed medium» 


To determine the value of the reduced mass and the coefficient of inelastic resistance of the 
concrete mixture interacting with the cylindrical walls of the mixing drum during horizontal vibra- 
tions in the direction of the coordinate, let’s conventionally divide the entire mixed concrete mix- 
ture located in the cylindrical part of the mixing drum into a number of elementary vertical volumes. 
the same width (Fig. 4.2). From Fig. 4.2: 


2h 


A (4.14) 


| 


The length of each i-th volume with the width A, located at a distance /A, from the drum shell 
is determined from the following relationship: 


L = /R? +A —iA,P. (4.15) 


Based on the obtained expressions (4.12), (4.13) and (4.15), let's determine the values of the 
reduced mass and the coefficient of inelastic resistance of the concrete mixture for each elementary 
volume interacting with the cylindrical walls of the mixing drum during horizontally directed vibrations: 
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hA/p 
m, tgk,L 
a 2k,{1+ uw) a 
= tok L. 
pre 21+ w) a 


4.3 EQUATION OF MOTION AND ANALYTICAL DESCRIPTION OF VIBRATION MIXER MOTION 


(4.16) 


(4.17) 


Based on the obtained expressions (4.10), (4.11), (4.18) and (4.19) of the calculation 
scheme shown in Fig. 4.2 and the method described above, vibrations of a mixing drum loaded with 


concrete can be represented in the form of the following system of differential equations: 


— vibrations in the horizontal direction along the coordinate axis: 


(m-+m oY 4 tb, ab se 


pro) pro) “+ C,y,, = Asinot; 


— angular (torsional) vibrations about the X coordinate axis: 


ey 


(m+M,,0) ea —+ +(b; + byl +C,X,, = Qcosat: 


— angular (torsional) vibrations about the X coordinate axis: 








(U, + Jy) Me +(n, + nl e +k w,. = Ql, sinot: 


— angular (torsional) vibrations about the Y coordinate axis: 


2 


dy. 
Fd) ie 








d 
+(N, + Ayo) = 3V yp = Al, coset; 


— vertical vibrations in the direction of the Z coordinate axis along the X axis: 


L yp = XW yp 


— vertical vibrations in the direction of the Z coordinate axis along the Y axis: 


Lup = YW es 


(4.13) 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


Where Xi. Vier Vans Z,p — linear displacements of the mixing drum in the direction of the coordi- 
nate axes X, Y and Z under the action of harmonic exciting forces Qcoswt and Qsinwt in the 
operating mode; y,,, y,, — angular displacements of the mixing drum relative to the coordinate 
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axes X, Y, respectively, under the action of the moments of exciting forces Qsin@t and Qcosae in 
the operating mode; /, — distance from the center of application of the exciting force of the vibra- 
tion exciter of vibrations to the center of gravity of the vibration dynamic system in the operating 
mode in the direction of the coordinate Z,,, 


| mi +m,.,(0.9h + a) (4.29) 
a M+ My | 


J,, — moment of inertia of the concrete mixture relative to the horizontal axis of gravity passing 
through the center of weight: 


1 
i= 5 Mora " qo Morel + Myo 


(0.5h+a-1); (4.30) 


— coefficient of inelastic resistance of the concrete mixture at angular displacements of the 
oscillating system relative to the horizontal axis: 


Ngo = O.9h0,,9. (4.31) 
The solution to the system of differential equations (4.18)—(4.21) for stationary oscillations 


describing the steady motion of a vertically installed mixing drum in operating mode is presented 
in the following form: 


y, (0) =A, sinlat — @,,J; (4.32) 
x,,(t) = A, coslwt + @,,J; (4.33) 
y,,(t) =, sinlot —§,_J; (4.34) 
y(t) = V5, coslat + €5,1, (4.35) 


where A,, — amplitude of harmonic oscillations of the oscillating center of gravity of the system in 
the direction of the coordinate Y in the operating mode; A,. — amplitude of harmonic oscillations 
of the center of gravity of the oscillating system in the direction of the coordinate X in the operat- 
ing mode; ‘¥,. — amplitude of angular (torsional) harmonic vibrations of the mixing drum relative to 
the axis X passing through the center of gravity of the oscillating system in the operating mode; 
‘Y,, — amplitude of angular (torsional) harmonic vibrations of the mixing drum relative to the axis Y 
passing through the center of gravity of the oscillating system in the operating mode; @,., @,,. — phase 
shift angle between the amplitudes of the exciting forces and the amplitudes of the forced vi- 
brations; €,., €,.— phase shift angle between the amplitudes of the moments of the exciting forces 
and the amplitudes of the angular forced vibrations: 
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A, = ——— (4.36) 
[c, —(m+ Mo) F + (b, + Byol @ 
A, = ——— (4.371 
Ic, —(m+ mo) + (b, + By) @ 
a (4.38) 
[k, —(J, +, J +(n, +b.) o® 
a (4,39) 
[k, (J, +d, J +(n, + B,.) 0° 
(05 Da) 
©, = arctg—2 (4.40) 
Cy —(m+ mM, 9) 
D5 D5) 
Oo, = arctg—2 ee (4.41) 
c, —(m+m,,9) 
(n, +B, Jo 
= arctg —*——*= __: (4.42) 
oP ep eet 
(n, + D,.J@ 
= arctg —> =>. (4.43) 
ee ary 


Using the obtained solutions (4.32)-(4.43) of the system of equations (4.18)—(4.28) of the 
considered dynamic system, let's successively determine the laws of motion of the bottom and 
walls of the mixing drum, as well as the protective cap, affecting the processed concrete mixture 
both in vertical and horizontal directions, and causing normal and tangential stresses in the pro- 
cessed medium, significantly affecting the destruction of internal bonds in the concrete mixture 
and, accordingly, the efficiency of vibration processing and mixing of the prepared concrete mixture. 

The law of motion of the bottom of the mixing drum interacting with the processed concrete 
mixture horizontally in the direction of the coordinate can be described by the following equation: 


X ,{t) = x,,()+0—ahy,,(t)=A,, coslwt + J+ ¥,,0/ — alcoslot + §,.J, (4.44) 


where X,.{£) — horizontal displacement of the mixing drum bottom in the direction of the coordi- 
nate X, which causes contacting stresses in the concrete medium. 

After transformations of expression (4.44), let's obtain a dependence convenient for analysis 
and computer simulation: 


X(t) = A,, coslat + 0,), (4.45) 
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where A,, — amplitude of harmonic vibrations of the bottom of the mixing drum in the direction of 
the coordinate axis, 


A, = VA, +[¥,,./-al +24, ¥,, / — allcoslo,. —&,.): (4.46) 


8, — angle of phase displacement between the amplitude of the exciting force and the amplitude 
of the forced vibrations of the bottom of the mixing drum in the direction of the coordinate axis X, 


A,.sing,, + Y,,.(/ — a)sin§g.,. 


(4.47) 
A,.cos@,, + ‘¥,.(/ — a)cosé,. 


0, =arctg 


The law of motion of the bottom of the mixing drum interacting with the processed concrete 
mixture horizontally in the direction of the coordinate can be described by the following equation: 


Y(t) = y,()+0-aly (()=A, sinlot —o,J+'¥,0/ —alsinfot —¢,), (4.48) 


where Y,.(€) — horizontal displacement of the mixing drum bottom in the direction of the coordi- 
nate Y , which also causes contacting stresses in the concrete medium. 

Transforming expression (4.48), let’s obtain a dependence convenient for analysis and com- 
puter simulation: 


Y, (t) =A, sinlot — 6,), (4.49) 


where A,, — amplitude of harmonic vibrations of the bottom of the mixing drum in the direction 
of the coordinate axis Y, 


= JA. +[Y, -al +24, (/—allcoslo,. —E,,); (4.50) 
d ‘Ir ir Ir ~ 1r Ar Ar 


8, — angle of phase displacement between the amplitude of the exciting force and the amplitude 
of the forced vibrations of the bottom of the mixing drum in the direction of the coordinate axis X, 


A, SIN@,, F Wop (¢ a a)sin€,. 


(4.51) 
A,.cOs@,, + W,,(¢- a)cosé,, 


0, =arctg 


The law of motion of the bottom of the mixing drum interacting with the processed concrete 
mixture in the vertical direction along the coordinate Y can be described by the following equation: 


Z,ly,t)=yw,.(t)= PY, ysintot —€,.) at -Asy<-rgersy<h. (4.92) 
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The law of motion of the bottom of the mixing drum interacting with the processed concrete 
mixture in the vertical direction along the coordinate X can be described by the following equation: 


Z,. (x,t) = xy, (E) = xw,, cos(@b+&,,) at -RSx<-r,r<x<FP. (4.53) 


Here Z,.(y,t) and Z,,(x,t) — displacement of the bottom of the mixing drum in the vertical 
direction along the coordinates X and Y, causing normal stresses in the concrete medium. 

The law of motion of the mixing drum shell interacting with the processed concrete mixture 
horizontally in the coordinate X direction can be described by the following equation: 


X (2,0) = x,,(0)+ zy,,(t) = 
=A,,coslot+@,,)+ ¥,.zcoslwt+é,.) at —(/—a)<z<l(H-/+ a), (4.94) 


where X,,.(z,¢) — movement of the shell of the mixing drum in the horizontal direction along the 
coordinate X , which causes normal stresses in the medium being processed. 

After transformations of expression (4.94), let's obtain a dependence convenient for analysis 
and computer simulation: 


X,.(Z,t) = A,,,(z)coslwt + 6,(2)], (4.59) 


where A,,, — variable amplitude of harmonic vibrations of the mixing drum shell in the direction of 
the coordinate axis, 


Asap a A, oT (wz) ue 2A; Wo,ZC08(@a, _ ex) (4.96) 


Q,(z) — variable phase angle between the amplitude of the exciting force and the amplitude of 
forced vibrations of the shell of the mixing drum in the direction of the coordinate axis: 


A, sing, + P,,ZSiNg,. 


(4.57) 
A,,COSQ,, + ¥, ZCOSE,, 


8,(z) = arctg 


The law of motion of the mixing drum shell interacting with the processed concrete mixture 
horizontally in the coordinate direction can be described by the following equation: 


rolZ ty, te 2w ltl 
=A,,sinlot —9,,J+¥,,zsinlot —€,,) at --alsz<(H-/+a), (4.58) 
where Y,,.(z,t) — movement of the shell of the mixing drum in the horizontal direction along the 


coordinate, causing normal stresses in the medium being processed. 
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Transforming expression (4.58), let's obtain a dependence convenient for analysis and com- 
puter simulation: 


Y 


obr 


(z,t)=A,,,(zsinlot — 0,(2)], (4.59) 


where A,,,, — variable amplitude of harmonic vibrations of the mixing drum shell in the direction of 
the coordinate axis Y, 


Ayn =A, +0Y,,2" +24, P,,2c08(@,, - &,,); (4.60) 


Q,(z) — variable phase angle between the amplitude of the exciting force and the amplitude of 
forced vibrations of the shell of the mixing drum in the direction of the coordinate axis Y, 


A, SiN, ES Y,,28iNG,, (4.61 } 


6, (z) = arct 
? d A,.cos@,. +P, Zcosé,, 


The law of motion of the shell of the mixing drum interacting with the processed concrete mixture 
in the vertical direction along the coordinates Y and X can be described by the following equations: 
— for the cylindrical part: 


Z,,,,(t) = Ry, (t) =, Asinlot - E,,): (4.62) 


obyr 


Z 


obxr 


(t) = Ry, (t) = ,,Rcosloot + &,,): (4.63) 


— for the tapered part: 


Z pot yt) = Ry, (0) = 
y-(h-£-a)] | 
=, |R a a sinlot +€,,) at (h-¢+a)<y<(H-£+a); (4.64) 
Lay (Xt) = Ay, (t) = 
=, c 2 101 east +€,) at (h-l+alsy <(H-I+al, (4.65) 
g 
where Z,,,,{¢) and Z,,,,{t) — displacement of the cylindrical part of the mixing drum in the vertical 


direction along the coordinates Y and X, causing tangential stresses in the concrete medium: 
Z4y(y,t) and Z,,, (x,t) — displacement of the tapered part of the mixing drum in the vertical direction 
along the coordinates Y and X. 
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4.4 DETERMINATION OF THE MAIN PARAMETERS OF THE WORKING PROCESS OF AN OPERATION 
OF CONCRETE MIXER 


Based on the obtained theoretical dependencies (4.1)—(4.65), describing the vibrations of the 
vibration drum in the operating mode, a program was created in the TurboPascal system designed 
to simulate the laws of motion of the mixing drum under the action of disturbing forces generated 
by the displaced vibration exciter of circular vibrations. 

Fig. 4.3—4.6 show the theoretical values of the vibration amplitudes of the surfaces of a vi- 
bration drum in contact with a concrete mixture of different consistency in the horizontal plane and 
in the vertical direction. Used cement-concrete mixtures with a cone subsidence OK=3.5—4 cm 
and a stiffness of 30 s, 60 s and 90 s. 

The data of the vibration amplitudes of the bottom of the mixing drum (Fig. 4.3, a) and the 
shell of the mixer drum (Fig. 4.3, 6) with horizontal (1) and vertical (2) movements in the operating 
mode are given. The calculations were performed for the preparation of a concrete mixture with 
a cone subsidence OK=3.5—4 cm. 



























































02 104 156 208 R, mm 0 30 100 150200 250 300 350H, mm 
b 






O Fig. 4.3 Amplitudes of vibrations of the surfaces of the mixing drum (mixture 
with cone subsidence OK=3.0—4 cm) 
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O Fig. 4.4 Amplitudes of vibrations of the surfaces of the mixing drum (mixture 
with a hardness of 30 s) 
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O Fig. 4.5 Amplitudes of vibrations of the surfaces of the mixing drum (mixture 
with a hardness of 60 s) 
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O Fig. 4.6 Amplitudes of vibrations of the surfaces of the mixing drum (mixture 
with a hardness of 9O s) 


When preparing hard concrete mixes, the distribution of amplitudes is somewhat different 
from mobile mixes. Yes, in Fig. 4.5, 4.6 show the results of calculations for mixtures with stiffness 
of 60 s and 90 s for the bottom of the mixer drum (a) with horizontal (1) and vertical (2) and 
the shell of the mixer drum (5) with horizontal (1) and vertical (2) movements. 

The data presented show that during the preparation process concrete mixtures are subjected 
to vibration action with sufficiently high amplitudes in the horizontal plane and in the vertical direc- 
tion at a frequency of angular vibrations of 292 rad/s, sufficient for effective mixing and vibration 
activation of the concrete mixture. 

Thus, the obtained expressions (4.1)-(4.65) make it possible to establish the law of mo- 
tion of all surfaces of the inner surface of the mixer interacting with the processed cement- 
concrete mixture during its preparation, to determine the main parameters of the proposed 
vibration device depending on the physical and mechanical characteristics of the processed 
mixture, to justify rational modes. actions on the concrete mix during its vibration processing 
and preparation. 

For experimental research, a vibration mixer was developed and manufactured (Fig. 4.7). 
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O Fig. 4.7 Design of the experimental device 


It is a cone-cylindrical mixing drum 1, which is mounted on a base plate 3 with the help of 
elastic shock absorbers 2. This base plate 3 is rigidly connected in its central part to a drive shaft 4 
installed in a bearing support 5. A worm is built into the cross-beam 6 a gearbox 7 connected to 
an electric motor 8 through a clutch 9. On the bottom of the mixing drum 1, a circular vibration 
exciter 10 is mounted, which is half mounted in the inner cavity of the mixing drum. Inside the 
mixing drum, the circular vibration exciter 10 is closed with a conical protective cap 11, which is 
hermetically attached to the bottom of the drum. Inside the mixing drum 1 there are blades 12. 
The traverse 6 and the locking and tipping mechanism 13 are mounted on the frame 14. 

The process of testing the operation of the mixer with the necessary instruments for mea- 
Suring parameters was carried out in accordance with the requirements of the interstate stan- 
dard (GOST 16349-85). 

— error of power measuring instruments — no more than 2.9 %: 

— rotational speed of the mixing bodies was checked under load: 

— weight of the mixer was determined by weighing on a balance; 

— total time of one mixing cycle was determined by a stopwatch and was counted from the 
beginning of loading the mixture to the end of unloading the mixed mixture: 

— composition, grade of the concrete mixture was determined at the test site: 

— during the tests, three batches of mixtures of the same composition were prepared. 
The amount of the mixture was sufficient for testing. A total of three mixtures were used. 

All indicators provided by the interstate standard, including vibration parameters (amplitude 
and frequency of oscillations), were recorded through the appropriate equipment on the computer 
display (Fig. 4.8). 

On the foundation on the side of the vibration mixer (Fig. 4.8) there is a fixed stand of the 
tripod, in which the sensor is fixed on a special clamp. The other end of the sensor is attached 
to the bottom of the mixing drum through a special bracket. The length of the sensor base and 
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its thickness are selected from the condition that the natural oscillation frequency of the sen- 
sor @,, exceeds the value of the oscillation frequency by 8...10 times. 

As a result of the experiments, the value of the vibration amplitude is recorded on the com- 
puter screen in the form of vibrograms. Based on the results of the analysis and processing of 
vibration records, both the mutual influence of the working body and the medium and the determi- 
nation of optimal modes are established. The analysis of the records was carried out in two stages: 
preliminary analysis of the oscillogram; complete control processing. 

By processing a series of typical vibrograms according to the above method, the numerical 
values of the vibration amplitudes of the bottom of the mixing drum and shell were obtained for 
mixtures with different cone draft and different stiffness. Based on the results of data processing, 
graphs were constructed with the repetition of each experiment in three dimensions. Fig. 4.9 
shows graphs comparing the theoretical and experimental values of the amplitude of oscillations of 
horizontal and vertical oscillations of the bottom and shell, depending on the radius, measured from 
the axis of rotation of the drum to the shell of the drum. 


Computer 














WZ 
O Fig. 4.8 Measurement scheme: X,,— horizontal movements 
of the mixing drum; X,— vertical movements of the mixing drum 




















O Fig. 4.9 Comparison of theoretical and experimental values of the amplitude: 
a — horizontal vibrations of the bottom and shell, depending on the radius: 
6 — vertical vibrations of the bottom and shell for the mixture G=30 s 
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An assessment of the performed studies has established that the discrepancy between the 
theoretical and experimental numerical values in the vibration amplitudes is 11-15 %. 


4.5 METHODOLOGY AND ALGORITHM FOR CALCULATING THE MAIN PARAMETERS 
OF A VIBRATION CONCRETE MIXER 


The initial parameters are the amplitude, the vibration frequency, set depending on the com- 
position of the mixture. An important indicator of the vibration process is the selected type and 
direction of vibrations and the nature of vibrations. If, according to the condition of the problem, 
the values of the amplitudes and frequencies of oscillations are absent, then their values can be 
determined based on the fact that the dominant parameter is acceleration. Then such properties 
as amplitude, speed, frequency, acceleration must be considered together. 

For harmonic oscillations with frequency, the oscillation amplitude is optimal [7]: 
(4-6)g 


OPT — D 
@M 


X (4.66) 
where g — acceleration due to gravity (g =9.8 m/s). For the frequency w = 3.14", the optimal 
vibration amplitude is determined by the lower and upper acceleration, i.e. X,,, =0.4—0.6 mm. 

The second parameter is the drive power of the concrete mixer, which is spent on vibrations 
of the mixture P, and friction in the bearings of the vibration exciter P.,. 

The vibration power is calculated according to the maximum value that the existing driving 
force can develop with an increase in the dissipative resistance of the system. 

Power per vibration, W, is determined by the dependence [7]: 


maxP, = na 0, (4.67) 


where F, — driving force necessary to maintain a given amplitude is calculated taking into account 
the effect of the concrete mixture on the working body of the machine; x’ — amplitude of displace- 
ments of the working body, taking into account only the reactive forces of the system: 


x= ig = Otte (4.68) 
ma+m k 


where w — frequency of forced vibrations. 
Friction force power in vibration exciter bearings, W, 


P,, = Fu— o, (4.69) 


where 1=0.005...0.008 — coefficient of friction in the bearings; d,— diameter of the bearing journal. 


72 


4 RESEARCH OF TECHNICAL SYSTEMS OF PROCESSES OF MIXING MATERIALS 


Motor power for vibration exciter, kW: 


_ max Peg Fe 


i | (4.70) 
1000n 


where wu — transmission efficiency. 
The main parameters of the gravitational action, as components of the general process of mix- 
ing the concrete mixture, is the drum rotation frequency, which is determined by the dependence: 


w = 29/0. (4.71) 


The algorithm for calculating the main parameters and characteristics of a vibration mixer 
based on gravitational and vibration components is shown in Fig. 4.10. 


Initial data 






The composition 
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a vibration mixer 
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by criterion K, 


Assessment 
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The final value of the 
parameters and characteristics 
of the vibration mixer 





O Fig. 4.10 Algorithm for calculating the parameters and characteristics of a vibration mixer 
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4.6 DISCUSSION OF RESEARCH RESULTS 


The analysis of the obtained expressions (4.382)-(4.65) shows that the mixing drum in the 
process of operation makes complex spatial vibrations, designed to provide an effective vibration 
action on the mixed concrete mixture. Moreover, each point located on the inner surface of the 
mixing drum has its own law of motion. As a result, the mixed concrete mixture will experience an 
alternating amplitude-frequency action, which causes multidirectional normal and tangential stresses 
in the mixture, which ensures the ultimate destruction of internal bonds and greater mobility in it, 
which contributes to an increase in mass transfer and effective mixing [8-13]. 

The use of a vibration concrete mixer makes it possible to reduce the mixing time by 
2.0—d times in comparison with conventional free mixing concrete mixers. Provides the preparation of 
hard concrete mixes and improving the quality of concrete mixes due to the destruction of defective 
aggregates and a more uniform distribution of the binder throughout the volume of the mixed mix. 
An experimental design of vibration-gravitational mixing has been developed and manufactured, and 
studies of the working process of mixing efficiency have been carried out on it. Research carried out 
at three warehouses of concrete mixtures confirmed the assumptions accepted in the work, the 
assumption and the selected physical and mathematical model of the «Mixer — concrete mixture» 
system. It was found that the use of vibration and gravitational action makes it possible to reduce 
the mixing time from 90 to 60 s. The strength of control concrete cubes is 17—18 % higher than in 
a gravity mixer. The discrepancy between the theoretical and experimental numerical values for the vi- 
bration amplitudes is 11-15 %. The proposed design and the studies carried out on it have confirmed 
the effectiveness of the vibration mixing method. As a result of the studies performed, it was found 
that the nature of the vibration of the mixing drum is close to harmonious [14]. The amplitude of verti- 
cal oscillations of the bottom varies in the range of 0.8—0.4 mm at a distance R=208—220 mm from 
the drum bottom. At a distance of R=120 mm, the vibration amplitude slightly decreased and its av- 
erage value is 0.2 mm. Horizontal displacements of the drum bottom have huge vibration amplitudes — 
0.02 mm. When loading with a concrete mixture, the amplitudes of vibrations at the same radius have 
significantly decreased, and their average value fluctuates between 0.2—0.25 mm and 0.14—0.16 mm. 


CONCLUSIONS TO SECTION 4 


1. Well-coordinated equations of motion of the system «mixer — concrete mixture» and their 
solutions established the regularities of changing the parameters of the mixer at different angles 
of inclination. 

2. The distribution of the amplitudes of vibrations of the working body of the vibration solvent 
has been obtained, depending on the distance of application of the energy source. 

3. An experimental design of vibration-gravitational mixing has been developed and manufac- 
tured, and studies of the working process of mixing efficiency have been carried out on it. 
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4. It has been determined that the use of a vibration concrete mixer makes it possible 
to reduce the duration of mixing by 2.0—-3 times in comparison with conventional free mix- 
ing concrete mixers. Provides the preparation of hard concrete mixes and improving the quality 
of concrete mixes due to the destruction of defective aggregates and a more uniform distribution 
of the binder throughout the volume of the mixed mix. 

0. The results of experimental studies are presented and a comparative analysis based on 
theoretical data has been carried out. 

It has been found that the use of vibration and gravitational action can reduce the mixing time 
from 90 to 60 s. The strength of control concrete cubes is 17—18 % higher than in a gravity mixer. 

6. An algorithm and method for calculating the main parameters of a vibration concrete mixer 
have been developed. 
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ABSTRACT 


The working process of concrete mixture compaction is investigated. The model of the pro- 
cessing medium is determined, which is characterized by three stages. At the first stage of 
compaction it is a free-flowing medium, at the second and third stages of compaction it is an 
elastic-viscoplastic medium. This result lies in the definition of the physical and mathematical 
model, the equations of motion. The vibration system «machine-medium» at any moment of 
movement consists of two components — one is capable of accumulating energy that passes 
from one form to another (reactance), and the second is energy dissipation (active resistance). 
The influence of active and reactive forces on the general movement of the system under study 
has been investigated and determined. By solving the equations, analytical dependencies were 
obtained, which established the laws of motion of the «Vibration machine — sealing medium» 
system. Practical implementation for technical systems for the formation of concrete products 
and soil compaction processes under the main pipeline has been carried out. An algorithm and 
method for calculating the main parameters of vibration systems for compaction of working 
media have been developed. 


KEYWORDS 


Vibration system, medium, compaction process, concrete mix, soil, physical, mathematical 
model, reactance, active resistance, equation, analytical dependences, parameters, experimental 
units, amplitude, vibration frequency. 


0.1 DETERMINATION OF PHYSICAL MODELS OF THE PROCESS OF VIBRATION COMPACTION 
OF MATERIALS 


The vibration effect on the concrete mix is of great practical importance and forms the basis of 
all modern technology for compaction of consistencies [1-3]. The essence of the vibration action 
lies in the fact that during vibrations, the concrete mixture acquires the properties of fluidity due 
to the violation of bonds between the particles. The compaction process is characterized by the fol- 
lowing movement of concrete mixture particles: individual, general (general) vibration of concrete 
mixture particles as a continuous system. The uniform movement of particles from one equilibrium 
position to another occurs predominantly in the surface layer. 
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Based on a general view of the compaction process, it can be noted that the concrete mix at dif- 
ferent stages obeys different laws. At first, it can be represented as a free-flowing medium, where 
the vibration movement of its particles occurs, that is, the forces of dry friction mainly affect. With 
Subsequent compaction, the mixture becomes more mobile and behaves like an elastic viscoplastic 
body. Along with the properties of a viscous liquid, the mixture is clearly elastic. Thus, the physical 
model of the mixture should take into account three stages of the compaction process (Fig. 5.1). 


Stages ofthe mixture 
compaction and its 
properties 


| property stage: 11 property stage: 111 property 
elastic, plastic elastic, viscous and stage: elastic and 
plastic viscous 





O Fig. 5.1 Stages of the compaction process 


The elastic properties according to Hooke’s law (Table 5.1) are presented in the form of 
a linear (linear model) or nonlinear (nonlinear model) layer of a concrete mixture. 

Viscous properties according to Newton's law (Table 5.1) are modeled by a damper consist- 
ing of a piston moving in a cylinder with a liquid. Other connections of elements with each other 
are possible, which is due to the need to take into account certain properties and composition of 
the components of a real mixture. The equations describing the propagation of a wave process in 
a mixture can be represented as follows: 


O*U 


9 b) p) 2 3 2 
du .»dU du_ »dU nN AU Sa ols SS. 
Z 


du _ pou. du nou | 5.1) 
oe OF a ae pda at 





The first equation (5.1) describes the displacement of the layer of the mixture ulz, t) as an 
elastic body with density at a speed c. The second equation (5.1) describes the motion of the mix- 
ture as an elastic-viscous body with density with a frequency-dependent dissipation of energy taken 
into account in the equation by the coefficient 7. The third equation (5.1) describes the motion of 
the mixture as an elastic — viscous body with density p with a frequency-independent dissipation of 
energy taken into account in the equation by the coefficient y. At different stages of the concrete 
compaction process (Fig. 5.1), all equations (5.1) can be applied. The first equation (5.1) is not 
applicable in the case of a possible resonance, since the absence of a dissipative component in the 
resonance zone gives a solution for the amplitude of oscillations A— oo. Therefore, the second and 
third equations (5.1) are left to choose from. At the present stage of development of vibration 
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technology for compaction of building mixtures, the third equation (5.1) is more applied due to the 
fact that for a specific vibration frequency, the stress state is determined by the vibration amplitude. 
Therefore, in further analytical studies, the third equation (5.1) is used. 


@ Table 5.1 [he main types of physical models 


Model 
properties 


Elastic 


Plastic 


Elastic-plastic 


Viscous 


Elastoviscous with 
stress relaxation 


Elastoviscous with 
delayed deforma- 
tions 


Elastoviscous with 
creep phenomenon 


Elasto-visco-plastic 


Elasto-visco-plastic 
with relaxation 


Elasto-visco-plastic 
in the process of 
shear deformation 


Physical model 
scheme 


-W—— 
At 1s 
pw 


Mathematical notation 


o=cE, tT=XxXG 


dV 

tT=XG, T=XG,+T,; 7 Ty +N 
dV 

= VSG SE C= ttl 


Hooke 


Saint-Venant 


Prandtl 


Newton 


Maxwell 


Voigt 


Kelvin 


Bingham 


Shvedov 


Generalized 
Kelvin and 
Shvedov 


When vibrations are transmitted from the working body to the mixture (Fig. 5.2), elas- 
tic waves propagate in the latter, which carry the energy that goes to the compaction of the 


mixture layers. 
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Energy Working Concrete 
El Oscillator |—> body = te 


O Fig. 5.2 «Vibration unit — concrete mix» block diagram 





Based on the above reasoning, the «machine — medium» vibration system by different signs 
of their behavior will be represented from the point of view of modeling as subsystems with dis- 
crete (machine) and distributed (medium) parameters. 

The next premise is based on the fact that the «machine — medium» vibration system at any mo- 
ment of movement consists of two components — one is capable of accumulating energy that passes 
from one form to another (reactance), and the second is energy dissipation (active resistance) [3]. 

Thus, the task of accounting for the medium to be processed is reduced to determining the 
active and reactive forces in the contact zone and introducing them into the equation of motion of 
the vibration unit. 

For this purpose, the «vibration unit — concrete mix» design diagram according to the me- 
thod [2] is shown in Fig. 5.3. 





‘uth, t) 





‘ulz, t) 





O Fig. 5.3 Design diagram of the «vibration plant — concrete mix» system: 

x — displacement; F(t) — external forced force; u(z, t) — longitudinal displacement 
of the current section of the concrete mixture column; A,_, — reaction of the 
concrete mixture in the section, provided that z=0; A,_, — reaction of the 
concrete mixture in the cross section of the column, provided that z=h 


80 


0 STUDY OF TECHNICAL SYSTEMS OF MATERIALS COMPACTION PROCESS 


The wave equation of displacement of the current section of the column of concrete mixture 
is taken in the form (9.1): 


OU 5, . OU 

— =c1+iy—. 

ot® Laer 

Here 0°u/dt* — acceleration of the concrete mixture layer; d°u/dz* — the second derivative 
of the displacement of the concrete mixture layer u along the coordinate of the mixture column z: 
c — speed of the mixture layer; y — coefficient of resistance of the concrete mixture, which 
characterizes the quantitative value of energy dissipation; / — imaginary unit showing the rotation 


between the elastic and dissipative components by the angle w/2. 
The solutions of the wave equation (5.1) will be used in the form: 


u(z,t) =(ue* +u,e" Je. (5.2) 


The displacement u is determined by the product of two functions, one of which depends on 
the arguments z(z) and T(t): 


z(z)=u,e" —u,e-”, (5.3) 
T(t)=e™. 


In solution (5.2) u, and u, are constants determined by the boundary conditions, k — complex 
wavenumber: k =(o.+ iB), where o., 8 — coefficients [1-3]: 


o=n: aes (5.4) 
C» C» 
where u = oa) aed — vity +1 
eater)" Y atey) 


The coefficient o determines the attenuation of the wave propagated in the layer of the me- 
dium, and the coefficient 8 determines the attenuation of the same wavelength. Indeed, if there 
is no resistance (y=0), then n=0, B = /c,, that is, the wave under this condition propagates 
without extinction, which corresponds to the first equation (5.1). 

The main task is to provide an analytical description of the dynamics of the «vibration unit — 
concrete mix» hybrid system. 

The design diagram of the vibration system (Fig. 5.4) is described by the equation obtained 
on the basis of the use of the reaction method [4-6]: 


(m+m,a)X + (b+m,@d)x + ¢,x = Fo-sinot. (5.5) 
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The solution to this equation gives an expression x = x, sin( ot + y) for determining the vibra- 
tion amplitude of a vibration machine, taking into account the entire complex of resistance forces: 


= —————————— (9.6) 


| (c —mw*)+ mca} + (boo + mwd) 


where a — wave coefficient of the reactance of the concrete mixture; d — wave coefficient of the 
active resistance of the concrete mixture: 


ash2oh + Bsin2Bh . cosin2ebh — Bsh2ah 
sf (5.7) 
h (x? + B*)[ ch2och + cos2Bh | h (a? + B°)[ cheoh + cos2Bh | 

Because of the task at hand, it is important to determine not only the influence of the nominal 
resistance forces on the amplitude value of displacement in the contact zone, but also to find an 
expression for determining the nature and magnitude of the distribution of the amplitudes of dis- 
placement over the layers of the mixture z = z(h). 

Using the solution of the wave equation, let's find an expression for finding the displacement 
of any layer of the mixture: 


x(Zt)\= 


x ,ch(h — z)(o, +B) A ie) oa 
chI h(«. + iB,) | 

The modulus of expression (5.8) gives a formula for the amplitude of the displacement of the 
mixture layer: 


x(z,t) = 


x,ch(h—z)(o, +/B,)e"* 


(3,9) 
ch| h( (a , + iB,) )| 

When to enter an additional element, for example when loading, the limit conditions are 
changed. Now the stress on the upper contacting surface will not be equal to zero. It is obvious 
that the surcharge, like the working body, in its motion will also perceive the influence of the 
active and reactive resistances of the medium. By analogy, the force of resistance to surcharge 
fluctuations can be expressed as: 


R(ht) = E5(1+ in 2 ). : 


iwt 
= pis| & Jafue —u,e”° = eit = Tell (ue — ue"), (5.10) 





The notation is accepted here k =w/c; m,=pSh; 6=a+i8. 
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0.2 INVESTIGATION OF THE PARAMETERS OF VIBRATION MACHINES FOR THE COMPACTION 
OF CONCRETE MIX 


To study the parameters of vibration machines for the compaction of a concrete mixture according 
to the obtained formulas, calculations were made according to the following variable parameters: the 
height of the mixture column h = 0.1...1.5 m; the frequency of the oscillations w = 100..300 1/c; 
the speed of propagation of vibrations in the concrete mixture C, = 30...50 m/s; the loss factor 
y =0.05...1.0; the ratio of the masses of the concrete mixture m, to the mass of the vibration 
machine m: m,/m=0.4...1.0; static pressure. Fig. 5.4 shows the results calculations. 



























O Fig. 5.4 The results of calculations of changes in the parameters of vibration machines 
for the processes of compaction of a concrete mixture: a — the influence of the ratio of the 
masses of the medium and the machine on the amplitude of oscillations; b — the dependence 
of the amplitude of the oscillations on the height of the mixture column; c — the influence 

of the wave number on the amplitude of the oscillations; d — the influence of the loss factor 
on the vibration amplitude of the working body 


The effect of the mass ratio of the medium and the machine on the vibration amplitude depends 
on the specific height of the mixture column. So, in zones far from resonance, these are the heights 
of the mixture column h=0.2 g and h=0.9 m (Fig. 5.4, a, curve 1 and curve 2), the amplitude 
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of the oscillations decreases. In the resonance zone h=O0.3 m (Fig. 5.4, curve 3), the amplitude 
of oscillations increases, which is confirmed by the graph (Fig. 5.4, 6). A noticeable effect on the 
amplitude has the speed of propagation of oscillations (Fig. 9.4), the value of which is included in 
the determination of the numerical values of the wave coefficients a and d (5.7). So, an error in 
determining the vibration speed by 10 % can lead to a discrepancy between the calculated and 
actual vibration amplitudes by 150 % (Fig. 5.4, curve 2). 

The degree of influence of active resistance forces (Fig. 9.4, d) is determined by the height 
of the mixture column. In the zone of the greatest value of the vibration amplitude for a height of 
h=0.3 m (curve 1, Fig. 5.4, d), an increase y leads to a natural decrease x,.For the height (curve 1, 
Fig. 5.4, d), the picture is the opposite: here an increase y leads to an increase in the amplitude 
of oscillations x... 

The change y from 0.1 to 0.3 led to an increase in the x, amplitude from 0.156 to 0.173 mm. 

The considered examples clearly confirm the taken into account the reactive and active resis- 
tance of the concrete mixture. 


0.9 STUDY OF THE PARAMETERS OF SOIL COMPACTION PROCESSES UNDER THE MAIN PIPELINE 


Main pipelines are rightly considered the energy arteries of Ukraine, because they provide en- 
ergy resources not only for the needs of the state's own economy, but also for the transportation 
of oil and gas to more than 15 European countries. The current state of ensuring reliable and safe 
operation of trunk pipelines requires the urgent development and implementation of new equipment 
and technologies for capital repairs of pipelines, the use of which will significantly accelerate the 
pace of repair work with a decrease in material and labor costs for their implementation [7-9]. 

A feature of the new technology is that after a complex of works on replacing insulation under 
the repaired pipeline, soil is poured and compacted under the pipe to such an extent that it excludes 
subsidence of the latter after repair. 

A schematic diagram of soil compaction under the pipeline by two surfaces symmetrically 
located relative to the longitudinal axis of the pipeline, synchronously moving against each other 
along the required trajectory, is shown in Fig. 5.5. In the process of compaction, there is a mutual 
influence of stress diagrams in the soil from both soil-compaction surfaces, which leads to the for- 
mation of a generalized stress diagram in the zone of the vertical soil mass under the main pipeline. 

The stress state of the massif depends on many factors, primarily on the magnitude of the ap- 
plied external loads, load modes, physical and mechanical characteristics of the compacting soil, the 
size of the soil mass under the pipe, the shape and size of the trench, and the design parameters 
of the soil compacting working equipment. 

The values of technological factors in this study are taken to be fixed. The diameter of the 
pipeline is taken D,=1220 mm, as the most common in Ukraine. The weight of one running meter 
of the pipeline is assumed unchanged and maximum in value, for the case when the pipe is filled with 
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oil g,=14.3 kN/m. The height of the soil mass under the pipeline (the height of the placement of the 
lower point of the pipe relative to the bottom of the trench) is taken no more than H=820 mm. 
The geometrical dimensions of the trench for pipeline repair D,.=1220 mm, the following: depth 
h=3000 mm, width of the trench at the top b,=4400 mm and at the bottom b,=2620 mm (Fig. 9.5). 

The speed of application of loads to the soil-packing surfaces V, reflects the mode of operation 
of the soil-packing equipment. It must be consistent with the overlap factor of the compaction 
zone [10, 11]. Soil moisture W has a significant effect on the characteristics of its compaction. 
Let’s believe that the compacted soil has an optimal moisture content W, [12, 13]. 

The structural and force factors that significantly affect the soil compaction process under 
the pipelines include the geometric dimensions of the sealing surfaces, their relative position in 
the trench relative to the pipeline and the pressure they create on the compacted soil mass. Flat 
undeformed surfaces (blades) of rectangular shape are taken as working surfaces. 

Determining the approach to the study of the stress-strain state of the soil mass under the 
pipeline, let’s note that in the elastic region the stresses are directly proportional to the value of 
the applied loads. With a known value of stresses in the soil mass from the action of each fixed 
pressure of the sealing blades on the soil, it is possible to determine the value of stresses for 
another value of the applied pressure to the soil mass [14, 15]. Thus, with a known stress in the 
soil mass under the action of a fixed pressure P created by the working surfaces, it is possible to 
find new stresses in the soil mass for other values of the applied pressure of the blades on the soil: 

os 5.11) 

Pp 
where P — fixed pressure created by the blades on the soil mass; o — known value of stresses in 
the soil mass from the action of a fixed pressure. 





b, 
“D,=355...1420 mm Oil, gas, ammonia 





O Fig. 5.5 Scheme of compaction of the soil mass under the main pipeline 
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It follows from this that the pressure P created by the blades on the soil can be excluded from 
the number of variable factors. 

since the pressure on the soil from the blade is distributed evenly both along the length and 
along its width, the problem of determining the stress-strain state of the soil under the pipe in 
the direction along the pipeline can be considered as flat. Consequently, the stresses in the soil, 
when they are determined within the framework of flat surrender, will depend on the value of the 
pressure under the blades, regardless of the width of the blades B. On the other hand, under all 
equal loading conditions (fixed pressures), when the width of the blade changes, the stresses in the 
soil mass will change proportional to the width of the blade B: 

ae: (5.12) 

B 
where o’ — stress at an arbitrary blade width; o — stress determined for a fixed blade width. 

Therefore, such a geometric parameter as the width of the blade B is also referred to as 
constant factors. 

Combining formulas (5.11) and (5.12), let’s obtain the dependence of the stress o’ for ar- 
bitrary P’ and B’ at a known stress o arising in a flat body with thickness B and a given load P: 
P’B 
PR’ 


, 


0’ =0: (5.13) 





Thus, on the basis of consideration and analysis of the factors affecting the stress-strain 
State of the soil under the pipeline, it was established that the unchanged (constant) ones should 
include: h= 3000 mm, 6,=4400 mm, 6,=2620 mm, H=820 mm, D,=1220 mm, g,=14.3 kN/m, 
hy=H+0.5D,, P=10, w=0.37. 

Consequently, the constant factors, the value of which must be determined experimentally, in- 
clude: the compression rate V,, the compaction delay time C,, the compaction overlap coefficient K,, 
the elastic modulus F, the optimal moisture content Wo. 

In the process of studying the stress-strain state of the soil under the pipeline, it is necessary 
to determine the effect on the characteristics of the stress state of the following parameters: 
blade angle o, blade length L and width B, blade pressure P on the ground, distance / from blade 
to pipe, as well as distance S between blades (Fig. 5.6). Determination of the rational values of 
these factors and their interaction in the process of soil compaction is the main task of the study. 

The essence of the study is to determine the effect of changes in the linear parameters L and 
B of the sealing blades, the pressure P they create on the ground, on the provision of the neces- 
sary stresses in the soil mass under the pipeline, taking into account the parameters a, / and S, 
characterizing the position of the blades in the trench relative to the pipeline of a given diameter. 

To construct a mathematical model of the stress state of the soil mass under the pipeline, 
let’s select the values of the variable factors characterizing the compaction process for the zero 
level and the limits of change of these values are the intervals of variation of the factors. The value 
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of the variation intervals is chosen based on the following considerations and requirements. The 
compression angle was chosen w=90°...70°, based on the values of the angle of the prism of the 
soil massif under the pipeline, at which a sufficient stability coefficient of the massif will be ensured. 
When choosing the intervals for changing the linear parameters of the blades / and B, we leave, 
first of all, from the position of ensuring guaranteed non-damage of the pipeline insulation coating 
during soil compaction L=450...650 mm. 

Based on the results of studies of soil compaction with stamps, the recommended wicth of the 
Stamp should be equal to half of its length, because it is with this ratio that the maximum depth 
of the active compaction zone, B=220...325 mm, is achieved during the compaction process. 
In further studies, let’s take the value of the blade width to be constant and equal to B=300 mm. 

Since it is known that the depth of the core during soil compaction [13] is equal to two 
dimensions of the smaller side of the stamp [11], and taking into account the predicted increase 
in the active zones as a result of their overlap as a result of the oncoming motion of the sealing 
blades, let’s take S=1200...1600 mm. The levels of the selected factors and the intervals of their 
variation are presented in Table 5.2. 


@ Table 5.2 Factor levels and intervals of their variation study of the stress state of the MP soil bed 


Independent variables 6 Variation . 
Upper +1 Basic 0 Low -1 intervals A/ 

X1 a, degree 90 80 70 10 

Xe S, mm 1600 1400 1200 200 

X3 L, mm 650 000 450 100 

X4 /;mm 200 200 150 00 


The design diagram of soil compaction under the pipeline is shown in Fig. 5.6. The geometrical 
dimensions of the investigated soil massif were determined on the basis of the dimensions b,, 0,, 
h of the trench profile, built by the pipeline opening machine D,=1220 mm. The height of the soil 
massif under the pipe H is taken in accordance with the size of the developed soil massif under the 
pipe H=820 mm. 

The soil mass, compacted under the pipeline, is bounded by two inclined surfaces symmetrical 
about the vertical axis passing through the center of the pipeline, the positions of which are regu- 
lated by the values of the compression angle X1 and the distance between the blades X2 (Fig. 5.7). 
It is accepted that the value of the parameter X2 is calculated horizontally, located relative to the 
lower point of the pipeline at a distance of 1/3H. 

In this case, the mass of compacted soil is limited from the top of the generatrix of the pipe of 
a given diameter D,, and from the bottom by the bottom of the excavated recess, which is located 
at a distance H from the lower point of the generatrix of the pipeline. Since the pipeline wall and 
trench walls have significantly higher elastic moduli than the compacted soil mass, they can be 
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considered undeformed surfaces. Therefore, in the zones where the compacted soil borders on the 
generatrix of the pipeline and the walls of the trench, it is legitimate to impose boundary conditions 


for the movement of the soil perpendicular to the undeformed surfaces. 





O Fig. 5.6 Design scheme for determining the stress state 
of the soil mass under the pipeline 
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O Fig. 5.7 Results of calculating the stress state of the soil 
mass under the pipeline for typical experimental plans 


Taking into account the above, for further research let's take soil compacted isotropic in volume. 
The sizes, which are determined (reviews) in the study, are equivalent stresses at the characte- 
ristic points of the soil mass along the axis under the pipeline. Typical points are considered (Fig. 5.7): 
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Y1 — the magnitude of stresses in the soil mass in the area under the pipeline; Y2 — the value of 
stresses in the middle part of the soil massif at a distance of 1/2N from the pipeline; Y3 — the 
magnitude of stresses in the lower part of the massif at the level with the bottom of the trench: 
Y4 — the value of stress in the soil at a distance of 1/4H from the pipeline. 

Determination of the rational values of the parameters X1...X4, B and P, their functional 
interdependence and influence on the stress state of the soil compacted under the pipeline made 
it possible to provide the necessary values of the soil compaction coefficients and the uniformity 
of its compaction. 

The visual presentation of the results of calculating the stress state of the soil [8] under the 
pipeline using the finite element method made it possible to evaluate the influence of the interaction 
of compaction nuclei in the soil mass created by the blades of the working body, under the action 
of various combinations of values of the initial parameters that determine the compaction process. 
Typical results of calculating the stress state of the soil indicate that there may be combinations 
of initial parameters at which the expected effect of superposition of compaction nuclei is not 
observed, the stresses in the soil mass at different depths under the pipeline are insignificant, 
therefore, the coefficients of soil compaction remain minimal. The combination of the initial parame- 
ters (Fig. 5.8) provides a more uniform soil compaction under the pipeline, which indicates the 
presence of the effect of the interaction of compaction nuclei in the soil. 

Thus, it can be argued that by changing the combinations of the initial parameters that determine 
the process of soil compaction under the pipeline, it is possible to control the interaction of soil com- 
paction nuclei, as well as to influence the value of the final indicators of this process [11, 12, 16-18]. 

As a result of processing the data of numerical experiments [10], in accordance with the 
folded matrix, mathematical models were obtained to determine the stress state of the soil at the 
characteristic points of the massif under the pipelines, which are predictive models. Their adequacy 
was confirmed by the Fisher test with a confidence level of 0.95 [11]. 








pg (0.08637 —0.01696- x, +0.02101- x, - 
ate, 7 (5.14) 
' PB’ | -0.02524- x, —0.00306- x? + 0.00367: x? 
5, = ~~ -(0.09741— 0.0027. x, -0.00798-x, +0.01312- x, +0.0024-x5); (5.15) 
PB’ 
S, = Spr (0.0725 -0.00508-x, — 0.01162: x, -0.00276- x?); (5.16) 
5, = —”..(0.09889 — 0.00642. x, +0.00605- x, -0.01498- x, +0.00306-x5). (5.17) 
PB’ 


Thanks to the predictive models (5.14)—(5.17), it became possible to control the process of 
interaction of the compaction nuclei by changing the linear dimensions of the blades, the parameters 
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of their installation, namely: the compression angle «=X1, the distance between the blades S=X2, 
the blade length L=X3 mm, the distance from the blade to the pipeline /=X4, blade width B. 

Based on the solution of the system of equations (5.14)-(5.17), a rational combination of 
the output parameters of the soil compaction process under the pipelines was found: w=/0°, 
S=1400 mm, L=600 mm, /=180 mm, B=300 mm, as well as the dependence for the charac- 
teristic points 1...4 of the change in the values of the compaction coefficients on the pressure of 
the blades on the ground is built (Fig. 5.8). 
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O Fig 5.8 The value of the pressure of the blades on the soil 
from the coefficient of its compaction 


This made it possible to obtain the value of the pressure of the blades on the soil P=0.2...0.35 MPa, 
at which the necessary values of the soil compaction coefficients under the pipeline will be provided 
within the range of K,=0.9...0.99. 

Rational values of the parameters of the working body and operating modes under which the 
interaction of soil compaction nuclei under the pipeline is guaranteed, adopted for subsequent 
research and for use in the design of working equipment for soil compaction under the pipeline. 


0.4 DISCUSSION OF RESEARCH RESULTS 


A certain model of the processing medium, based on a general view of the compaction pro- 
cess, is characterized by three stages. This is at the first stage of compaction of a granular me- 
dium (Fig. 5.1), at the second and third stages of compaction, an elastic-viscoplastic medium. 
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This important result requires a clear definition of the physical and mathematical model, which 
is suggested by equations (5.1). Further studies are required to determine the temporal as- 
sessment of the course of these stages. Application of the idea of considering the «technical 
system — medium» system as a system that at any moment of movement consists of two com- 
ponents. One is capable of storing energy that passes from one form to another (reactance), 
and the other is capable of dissipating energy (resistance). A certain significant influence of 
the components of resistance on the movement of a technical system (Fig. 5.4) testifies 
to the hypothesis put forward about the need to take into account the forces of resistance 
of the medium. A technical system is a system with discrete parameters, and a medium — 
with distributed parameters. In this case, this general system is reduced to a computational 
system in the form of a system with discrete parameters, in which the wave phenomena of 
the medium are preserved and are represented by a contact force. The obtained wave coef- 
ficients (5.7) have general properties. They can be used in any equations of motion for such 
technical systems. This is the novelty of the method of analytical study of vibration systems for 
compaction of mixtures. 

The presentation of the results of calculating the stress state of the soil under the pipeline 
using the finite element method made it possible to assess the influence of the interaction of 
compaction of nuclei in the soil mass. By changing the combinations of the initial parameters that 
determine the process of soil compaction under the pipeline, it allows to control the interaction 
of soil compaction nuclei, as well as to influence the value of the final indicators of this process. 
Based on the results of studies of soil compaction with stamps, the width of the stamp was de- 
termined (B=225...325 mm). Its value is equal to half of its length, since according to this ratio, 
the maximum depth of the active zone of compaction is achieved during the compaction process. 
A rational combination of the initial parameters of the soil compaction process under the pipe- 
lines (a=75°, S=1400 mm, £=600 mm, /=180 mm, B=300 mm) has been established under 
which the interaction of soil compaction nuclei under the pipeline is guaranteed. 


CONCLUSIONS TO SECTION 5 


1. The model of the processing medium, characterized by three stages of compaction, has 
been determined. 

2. The review and analysis of physical and mathematical models of the process of compaction 
of technological media has been carried out. 

3. A model of the «vibration plant — medium» hybrid system is proposed, in which active and 
reactive resistances to wave coefficients are taken into account. 

4. The results of calculations of changes in the parameters of vibration machines for the 
processes of compaction of concrete mixture have been obtained, confirming the implementation 
of accounting for the reactive and active resistance of the concrete mixture. 
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0. Predictive mathematical models have been created to determine the stress state of the 
soil during compaction at the characteristic points of the massif under the pipelines, which make 
it possible to obtain rational values of the parameters of the working bodies for soil compaction. 

6. Regularities of movement of technical systems for compacting concrete mixtures and soils 
and their rational parameters have been established. 

7. An algorithm and methodology for calculating the main parameters of vibration systems for 
compaction of working media have been developed. 
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RESEARCH OF THE PROCESSES OF ACOUSTIC CAVITATION 
TECHNOLOGY FOR PROCESSING DISPERSED MEDIA 


|. Nazarenko, |. Bernyk 
ABSTRACT 


The operation process of acoustic treatment of technological media is investigated. The 
assessment and substantiation of methods for studying the parameters of acoustic treatment of 
technological media is done. The functional relationships between the acoustic parameters of the 
cavitation apparatus and the rheological properties of the processing technological media have 
been revealed. The process of staged acoustic treatment of technological media is described. 
A number of criteria and key parameters were determined, the use of which was carried out 
as an assessment in the calculation algorithm, depending on certain known initial data of the 
cavitator and the medium. The values of the input resistance of the compensator are determined 
and the condition for the maximum transfer of energy from the acoustic apparatus to the tech- 
nological medium is obtained. To calculate the parameters of the «cavitator — technological me- 
dium» ultrasonic cavitation system, an algorithm and method for creating a synergistic system 
were developed. 


KEYWORDS 


Acoustic system, cavitator, medium, bubble, stage, cavitation process, physical, mathematical 
model, rheological properties, equations, analytical dependences, parameters, amplitude, vibration 
frequency, intensity, energy, wave resistance. 


6.1 ASSESSMENT AND JUSTIFICATION OF THE CHOICE OF METHODS FOR STUDYING 
THE PARAMETERS OF ACOUSTIC TREATMENT OF TECHNOLOGICAL MEDIA 


Cavitation technology occupies a leading place in a number of advanced and effective me- 
thods for processing liquid media and creating new materials. This technology is being imple- 
mented for a number of different technological processes, including dispersion, emulsification, 
homogenization, degassing, and others [1-6]. Cavitation refers to vibration processes widely 
used in food [1-3], chemical [4, 5], mechanical engineering [6], construction [7], and other 
industries. The generally accepted definition of cavitation as a physical phenomenon is the pro- 
cess of the formation of cavities (of a certain size of bubbles) in a fluid flow under the influence 
of external pressure. Depending on the method of pressure formation in the medium, cavitation 
can be hydrodynamic or acoustic. 
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In hydrodynamic cavitation, a decrease in pressure occurs due to an increase in the local 
speeds of the liquid flow [5]. If a decrease in pressure is due to the passage of acoustic waves 
of a sound or supersonic spectrum of vibration frequencies, cavitation is considered acoustic [6]. 

During cavitation, the energy density of the sound field is transformed into a high energy 
density inside and around the bubble, which collapses over time. Consequently, energy is spent 
on the formation of shock waves, thermal energy, local electrification of bubbles, excitation of 
sonoluminescence, and the formation of free radicals [8]. In this case, during cavitation, there is 
a change in the rheological properties of the technological medium, such as viscosity, plasticity and 
elasticity [9]. Therefore, the determination of the parameters of acoustic treatment of technolo- 
gical media and energy requires the use of a computational model [10], which takes into account 
these changes. The key parameter of the evolution of gas and air bubbles in the acoustic field is 
the pressure that is formed in the contact zone of the «cavitator — medium» system [11-14]. 
The specified parameters and characteristics (Fig. 6.1) [15, 16]. 





Parameters of the 
cavitation process 


Energy 
parameters 





Power Medium 
parameters parameters 


p Energy of the Density, viscosity, 
ae cavitation elastic modulus, 
frequency apparatus, energy wave propagation 
viratian: of cavitation, speed, cavitation 
amplitude energy of shack volume 
waves 


Associated process parameters 


Action time, Bubble Cavitation 
natural frequency radii index / 


O Fig. 6.1 Parameters of the development of the cavitation process 
when processing technological media 
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The dominant influence on the sequence of the cavitation process is exerted by the following 
acoustic parameters and properties of the medium [17, 18]: 


F=f{AA.fFvW,P.L¢1,p.0,c,£, pu}, (6.1) 


where F — functional (integral criterion for evaluating a process); A — amplitude of oscillations of 
the contact zone «acoustic apparatus — medium» zone; A, — flow amplitude of medium oscillations 
at a distance x, from the boundary of the contact zone «acoustic apparatus — medium»; f — vibra- 
tion frequency of the acoustic apparatus; v — speed of vibrations of the contact zone «acoustic 
apparatus — medium»; W — energy, P — power; L — intensity; t — time; /— characteristic size of 
the medium in the direction of propagation of the acoustic wave in it; p — the pressure on the 
medium; o — stress in the medium; c is the speed of propagation of an acoustic wave in a medium: 
E — modulus of elasticity of the medium; p — density of the medium; wu — coefficient of viscosity 
of the medium. 

Integral parameters of function (6.1) are energy, intensity and power. Ultrasonic energy 
from the emitting surface of the ultrasonic apparatus to the medium undergoes several stages 
of transformation: 

— accumulation of energy in cavitation bubbles during their formation and expansion: 

— transformation of energy into energy of a shock wave during a burst of a cavitation bubble. 


6.2 DETERMINATION OF FUNCTIONAL DEPENDENCIES BETWEEN THE ACOUSTIC PARAMETERS 
OF THE CAVITATION APPARATUS AND THE RHEOLOGICAL PROPERTIES OF THE PROCESSING 
TECHNOLOGICAL MEDIA 


Let’s consider the influence of the parameters of the cavitation process on the formation of the 
volume of the cavitating technological medium. By analogy with the general functional (6.1), the rate 
of change in the volume of bubbles V can be represented by the following functional dependence: 


V(t) = ft, y,/,0,P),by,O,PysbtysPyoP,: (6.2) 


Taking time t, bubble radius A,, density p, as the main independent parameters in terms 
of dimension and using the basic provisions of the theory of dimensions [19], let’s reduce all pa- 
rameters (6.2) to a dimensionless form with respect to the formed bubble volume: t — time, IT]; 
| — characteristic length in the direction of the wave motion of the bubbles, [L]; A, — initial bubble 
radius, [L]; n — the number of bubbles per unit volume, [1/L°]; p, — density of the medium, [M/L']: 
1, — viscosity, [M/LT]; o — surface tension, [M/T*]; p, — gas density, [M/L*]; 11, — gas visco- 
sity, [M/LT]; p, — external pressure, [M/LT*]; p, — pressure inside the bubbles, [M/LT*]. Then the 
functional dependence (6.2) can be represented in dimensionless form: 
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J = finn He a eee er ee } (6.3) 
4 eR? /t an p, Ry /t pase P, p, Ry /t p, Re /t* p, Ry /t* 
3 


The next analysis procedure (6.3) is to assess the importance of taking into account the com- 
ponents of this dependence on the rate of formation of the volume of the cavitating technological 
medium. If to accept the condition of bubble bursting time t, in the form [20]: 


t, =0.915R, |e, (6.4) 





then dependence (6.3) can be re-formed into the following form: 


2 
="! ny; — ay at nee =| (6.5) 
3A | Te 
For media with low viscosity at the stage of bubble burst, according to the data of [20], the 
viscosity and surface tension can be neglected. Then dependence (6.5) for such conditions can be 
simplified to the form: 


2 
——-((! N: | (6.6) 
4 oR? |e, P.-D; T 
3 


The obtained dependence testifies to the importance of the influence of pressure, bubble diame- 
ter and burst time on the rate of formation of a volume of a cavitating medium with low viscosity. 
For other technological media, dependence (6.3) is used. The calculations of the numerical values 
of the volume of the cavitation zone depending on the intensity are given in Table 6.1. 


@ Table 6.1 Numerical value of the formation of the volume of the cavitation zone, depending on 
the intensity 


Intensity, W/cm? 5.0 10.0 20.0 30.0 40.0 50.0 
The volume of the cavitation zone, m?-10~ 0.10 0.30 0.40 0.42 0.44 0.46 


The importance of the data in this table lies in the fact that the ratio of the volume of 
non-bursting AV, (long-lived bubbles) to the volume of bubbles that burst AV, is an estimate 
of the efficiency of the cavitation process: 


K,=AV/AV,. (6.7) 
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In general, a number of criteria and key parameters have been determined, the use of which is 
proposed as an estimate in the calculation algorithms, depending on certain known initial data of the 
cavitator and the medium (Table 6.2). of the cavitation process of processing technological media The 
main set of criteria involves the assessment of the energetics of the process as a key characteristic, 
and taking into account possible changes in key parameters from the stage of nucleation of bubbles 
to the final stage — a burst of bubbles. Criteria that directly determine the reliability of the choice of 
the mode (harmonic or pulsed mode) are extremely important. These include criteria 8, 9 and 12. The 
importance of criterion 8 lies in the reliability of the selected speed of the cavitation apparatus, and 
criterion 9 is a direct indicator of the ratio of the burst time and the period of the technological pro- 
cess. The most numerical value of this criterion determines the provision of the minimum value of the 
ratio (6.7) and the possible expediency of using a pulsed mode of loading on the technological medium. 


@ Table 6.2 Criteria and key parameters for evaluating the effectiveness of the cavitation process of 
processing technological media 








No. Criterion, parameter Analytical dependence 
1 2 3 
1 The intensity of the cavitation process by the type of vibration 
1.1 Sinusoidal vibrations, /, W/cm? : 
inusoidal vibrations, /, W/cm pe ova 
Cpc 
1.2 Two-frequency oscillations, /, W/crn* oxA (1+kn)’ 
7 an 
1.3 Multifrequency oscillations, /; W/cm ox A : (1 ‘i ¥ kn) 
~ An’T 2 
1.4 Non-linear (non-sinusoidal) fluctuations, axA 
1, W/em? 7 ArT 
2 oystem synergy coefficient (efficiency), k;  k,=E,/E,.. 
3 The rate of change in the stress state of do 
the medium in time v,,, N/m’s ss Ee 


4 Coefficient of the ratio of the volumes of K,=AV,/AV, 
long-lived and cavitating bubbles, K, 


0 Coefficient of energy change over time for u=1/TUnE, /E,) 
the entire period of bubble processing, pu, u/s 





6 Energy absorption coefficient, ap, m~' Nostell 
7 Coefficient of the ratio of wave resistances 1/2 
eae ‘ 
of the cavitation medium, k 0,C 
; k =the | KB 
P pc, 1+ : 
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@ Continuation of Table 6.2 


17 


18 


19 


20 


21 


2 3 
speed ratio, k, K=V,/0) 


Coefficient of the ratio of the burst time to k,=1/T 
the time of the oscillation period, k, 


Contact zone work, A, J A,, = TM, X,O°U 
Work on the formation of cavitation A. 
bubbles, A,,, J Ay =a TAF 
Specific work of the cavitator, A, J A=P.-t 
Specific work of volume, A, cavitating A, =E,AV 
medium, J 
Coefficient of conversion of acoustic energy Ea 
into energy of shock waves, n, nes 

k 


Productivity of the cavitation process, AP, m*/s AP=kP/E, 








Bubble displacement amplitude, A,, ium 4 [9] 
(a oO pc 
Bubble speed amplitude, v,, m/s By 
v= pe 
Bubble acceleration amplitude, a,, um/s* By 
a, =@,/— 
pc 
Natural bubble oscillation frequency: 5 2 
discrete model, w,, rad/s oO, = 22 (R + “4 man ss |-(2) 
pF, Ry py pA, OY) 
Natural bubble oscillation frequency: contin- n \(El 4 iE" )(t ) 
uous model, w,, rad/s Oo, =—,/———— 
al) p(t.) 


Relative bubble size, kr K,=R.,/Ry 


The application of this or that criterion is due to the initial information and the set goal formu- 


lated before the researcher to select the model and parameters, or by the engineer to create the 
necessary acoustic wave system for a specific technological processing process or create a new 
technological medium. 


A similar example can be given for creating the required acoustic wave system. To implement 


the criteria-based assessment of the justification of the rational choice of the structural-parametric 
system «acoustic apparatus — technological medium», an algorithm has been developed (Fig. 6.2). 


The essence of the algorithm lies in the ability to vary not only the initial parameters and the 
layout of the cavitator relative to the processing medium, but also to determine the influence of 
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variable parameters on the maximum value of a particular criterion (Table 6.2). Depending on the 
formulation of the problem, a criterion is selected, which is fixed in block 8. Further, in block 2 
«Formation of the technological scheme», preliminary calculations are made for the mode of ener- 
gy transfer to the medium in accordance with the selected criterion. An important stage of the 
algorithm is blocks 3 and 5, which determine the physical and mathematical models based on the 
research results [7—9, 12-18]. It is in blocks 4, 6, 7 that the initial data are formed to determine 
the numerical values of the impact parameters and the limits of their rational use. The completion 
of the calculation is the parameters that serve as the initial information for making decisions to 
improve the rational design and technological parameters of the acoustic cavitator. 


1.Output design diagram, 
parameters of the 
cavitator and medium 


2.Formation of 
a technological scheme 
3. Physical model of 
the «K-S» system 













4. Formation of the 
parameters of 
the physical model 





6. Setting the boundary 
values of parameters 


7. Determination of parameters 
for process evaluation 





5. Mathematical model 
of the «K-S» system 













8. Criterion type and 


definition of its meaning: K,, 


10. Formation of an 
algorithm for calculating 
the acoustic system 


O Fig. 6.2 Scheme of criteria-based assessment of the justification of the rational 
choice of the «acoustic apparatus — technological medium» structural-parametric system 
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This technology of the proposed algorithm is a fundamentally new system of automatic enume- 
ration of parameters to fulfill the condition for ensuring the criterion of block 9. The second impor- 
tant aspect of this algorithm is the creation of a control system for the processing of objects by 
an ultrasonic cavitation unit in the optimal operating mode. The energy absorption coefficient (cri- 
terion 6), which characterizes the change in the specific energy P,, and the ultrasound intensity / 
in the irradiated medium, serves as a criterion for evaluating the efficiency of the acoustic treat- 
ment of this medium. 


6.3 DETERMINATION AND ASSESSMENT OF THE INFLUENCE OF THE WAVE IMPEDANCE 
OF THE MEDIUM ON THE ACTION OF THE ACOUSTIC APPARATUS 


An important parameter of the cavitation process is the wave resistance of the medium of the 
acoustic apparatus. The characteristic impedance is the ratio of the acoustic pressure P, to the 
vibration speed of the particles in the medium v: 


Z,=+P /v. (6.8) 


The importance of knowing the coefficient Z, is that it is an important characteristic of the 
interaction of the «acoustic apparatus — medium» system, not only as a resistance, but also as 
a characteristic of wave motion. Therefore, the obvious task of an effective acoustic device is its 
design, which ensures the transfer of maximum energy for the implementation of the cavitation 
process. That is why the signs should be taken into account in formula (6.9): the upper (plus) iden- 
tifying the wave transmitted from the radiation surface into the medium (incident wave), and the 
lower (minus) for the wave traveling in the opposite direction, i.e. from the border of the system 
«apparatus — technological medium» (wave). Obviously, for efficient transmission of the energy of an 
acoustic wave, the properties of the apparatus and the medium must be consistent with each other. 

This consistency consists in the fact that their wave supports must be the same, which is 
practically impossible in real conditions of the cavitation process. Indeed, the acoustic resistances 
for solid materials from which the emitting surface of the acoustic device is made, and the resis- 
tance of the technological medium, differ significantly from each other. 

For example, such a hard material as sapphire has an acoustic resistance for longitudinal waves: 
Z,=44.3-10° kg/(m*:s), and for water Z, =1.5-10° kg/(m*-s), for air Z, =4.27-10° kg/(m*. s) [21]. 
Comparing these digital values, it becomes obvious that in order to solve the problems of rational 
design of an acoustic apparatus and determine its parameters, it is necessary to analytically inves- 
tigate the sequence of acoustic matching of the apparatus and the technological medium in order to 
obtain the conditions for maximum energy transfer to the processing medium. For this, the process 
of propagation of an acoustic wave from the radiation surface of the apparatus to the processing 
medium is considered (Fig. 6.3). 
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reflected wave 


boundary 





O Fig. 6.3 Passage and reflection of the incident wave at the interface of the 
«emitting surface of the cavitation apparatus — process medium» system 


Let’s assume that a plane acoustic wave along the X-axis from the apparatus to the boundary 
with the medium propagates with acoustic resistance Z,, and in the medium at the border with the 
radiation surface of the apparatus along the same X-axis, due to the action of this resistance, wave 
resistance Z. arises. It is quite obvious that, in addition to the wave transmitted into the medium, 
a wave moving in the opposite direction arises at the interface between the radiation surface and 
the medium. Thus, a complex wave field is formed, which can be represented by the transmission 
and reflection coefficients of the wave in the form: 


K et = P,,.-/P. ainc } oe / ae : (6.9) 


By analogy with solving the equation of motion, it is possible to write expressions for the inci- 
dent and reflected waves of pressure and speed: 


P(x) = Pe + Pre) D(X) = Dying + Vpges (6.10) 
where k, = @, /c, — the wavenumber; / — imaginary unit. 

Applying the condition of continuity of motion parameters in the contact zone of the «appara- 
tus — medium» system (no violation of continuity): 


P,=P,; 0, = V,, (6.11) 
let’s obtain: 

Pie Paat Fos, (6.12) 

Z, Z, Z 
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Dividing the left and right sides by the amplitude of the acoustic pressure P._ and using (6.9) 


ainc 


let’s obtain expressions for determining the reflection and transmission coefficients of the wave: 
Ko= (Z, -Z,)/(Z, +Z,); Kas = 2Z,/(Z, +Z,). (6.13) 


The condition for determining the wave impedance for its maximum transmission can be rea- 
lized by introducing a compensator between the boundary of the radiation surface of the apparatus 
and the medium, which plays the role of resistance to the load and the reflection of acoustic waves 
exclusively in the zone of the apparatus (Fig. 6.4). 














A 1 
Z 
a 
| Zz . 1 — cavitator; 
V4 : 7 2 — compensator; 
3 3 — medium 
¥ Z. ; / 


O Fig. 6.4 Passage and reflection of the incident wave in the system «radiation 
surface of the cavitation apparatus — compensator — technological medium» 


In this case, the load resistance must be equal to the wave resistance of the medium: Z, = Z,. 
Using solution (6.10) with the replacement x = /, where / — the length of the radiation surface 
of the apparatus, where the wave propagates, let’s obtain the resistance in the section (x = —/): 





ik,| ay 
Ae ee ee fede (6.14) 
v{—/) eM _k 5 
Taking into account expression (6.14): 
Fee ey ac (6.15) 
Z,cosk,! + iZ, sink,| 


If to assume that a compensator with a length is installed between the radiation surface of 
the apparatus and the medium / = A/4, where A — the wavelength in the compensator, then 
let’s obtain the following dependence: 


kl = (2n/A)/(A/4) = 1/2. (6.16) 
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substituting (6.16) into (6.15), let’s obtain an expression for determining the input wave impedance: 
Z ANA) SZ, 12. (6.17) 


So, choosing the value of the input resistance of the compensator with a length of A/4, let’s ob- 
tain the condition for maximum transmission, according to which the characteristic impedance of the 
apparatus and the compensator is in agreement. Indeed, by placing an auxiliary layer of material with 
such an acoustic impedance between the boundary of the apparatus and the medium, it is possible to 
ensure the equality of the acoustic impedance of the apparatus and the equivalent transmission line. 
Then the waves of reflection from both boundaries of the additionally installed layer will be equal to the 
amplitude and will move in antiphase direction, which will lead to their mutual compensation. And the 
shift between the phases of movement for 180° is provided by the difference in wave path between the 
boundaries at half the wavelength. Equality of amplitudes is ensured by the optimal choice of the wave 
impedance of the compensator. If the compensator has its own parameters: resistance Z,, wavenum- 
ber k, and thickness /k, then provided that /, = 2, /4 is the input resistance at the boundary x = —I, : 


Z,=2°, /Z,. (6.18) 


To fulfill the condition under which there were no reflected waves from the section x = —I., 
the input impedance (6.17) must be equal to the acoustic impedance of the apparatus Z, : 


L272. (6.19) 
From (6.19) let’s obtain an expression for calculating the acoustic resistance of the compensator: 
ee aay ae (6.20) 


To apply dependence (3.72), which allows matching the acoustic impedance of the apparatus 
and the medium, it is necessary to determine the real numerical values of the wave impedance of 
technological media under conditions of their cavitation treatment. 

One of the most important parameters, on which the wave impedance of the technological 
medium depends, is the wave propagation speed c. In terms of physical content, the dependence 
of the wave propagation speed determines the ratio of elastic (E) and inertial (mass p) characte- 
ristics of the technological medium: 


c= JE/p. (6.21) 


An experimental device was created to measure the speed of wave propagation in a techno- 
logical medium (Fig. 6.5). 
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O Fig. 6.5 Experimental device for measuring the speed 
of propagation of waves in a technological medium 


The essence of the device is as follows. An electrical signal of a fixed ultrasonic frequen- 
cy is referenced from the pulse generator to the transmitter. The sensor, which in this case 
is a transmitter and receiver, converts the incoming signal into elastic vibrations of the same 
frequency (inverse piezoelectric effect) passing through the layer of the technological medium. The 
receiver converts elastic vibrations of ultrasonic frequency into an electrical signal of the same 
frequency (direct piezoelectric effect) and sends the signal to the computer screen. The computer 
screen displays the generator impulse and the signal transmitted through the tube with the pro- 
cess medium (Fig. 6.6). 





O Fig. 6.6 Vibrograms for measuring acoustic parameters: a — water: f=2.98 MHz; b — f=52 kHz 


The distance between them h characterizes the time t of passage of elastic vibrations from the 
source to the reflector and back. Consequently, the speed of propagation of waves in the techno- 
logical medium is determined by the formula: 


c=h/t. (6.22) 
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6.4 METHODOLOGY AND ALGORITHM FOR DETERMINING THE RATIONAL PARAMETERS 
OF THE WORKING PROCESS OF ACOUSTIC TREATMENT OF TECHNOLOGICAL MEDIA 


An algorithm was developed to calculate the parameters of the «cavitator — technological 
medium» ultrasonic cavitation system (Fig. 6.7). 


Algorithm for 
calculating 


|. Output data 


ll. Calculation of acoustic} |IIl. Calculation of rheological 
parameters characteristics 












IV. Intensity 





V. Process parameters 
LiacVetab 


O Fig. 6.7 Algorithm for calculating the process parameters 
of the «cavitator — technological medium» ultrasonic cavitation system 


The initial parameters used for the calculation are formed depending on the specific processing 
process or the creation of a new technological medium. The first step in the process of selecting 
the initial parameters is to substantiate the general scheme of the structural-parametric system 
«acoustic apparatus — technological medium» (Fig. 6.1) and the mode of action of the transducer 
on the medium. The structural-parametric scheme selected according to the algorithm (blocks 2—4, 
Fig. 6.1) requires specification of the limiting values of the parameters (block 6, Fig. 6.1) and the 
choice of a mathematical model of the «K-S» system (block 5, Fig. 6.1). In the future, the proce- 
dure for formulating the initial data, the calculation of acoustic parameters and the calculation of 
the rheological characteristics of the medium are carried out [4, 7-9, 11-15, 18-20, 22-26]. 
The criterion for the reliability of the calculation is the intensity (block IV, Fig. 6.7). When the 
numerical value of the intensity is provided, the procedure for calculating the process parameters 
is completed (block V, Fig. 6.7). 
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The value of the input resistance of the compensator of length A/4 is determined and the condi- 
tion of maximum transmission is obtained, according to which the wave impedance of the apparatus 
and the compensator is in agreement. By placing an auxiliary layer of material with such an acoustic 
impedance between the boundary of the apparatus and the medium, equality of the acoustic impe- 
dance of the apparatus and the equivalent transmission line is ensured. Then the waves of reflection 
from both boundaries of the additionally installed layer will be equal to the amplitude and will move in 
antiphase direction, which will lead to their mutual compensation. And the shift between the phases 
of movement for 180° is provided by the difference in wave path between the boundaries at half the 
wavelength. As a result, the maximum transfer of energy from the acoustic device to the techno- 
logical medium is created. For the implementation of the criterion assessment of the justification 
of the rational choice of the «acoustic apparatus — technological medium» structural-parametric 
system, an algorithm is proposed (Fig. 6.1). The essence of the algorithm lies in the ability to vary 
not only the initial parameters and the layout of the cavitator relative to the processing medium, but 
also to determine the influence of variable parameters on the maximum value of a particular crite- 
rion (Table 6.2). The necessary degree of influence and consideration of stages, rational modes and 
parameters for the development of techniques that implement an increase in the efficiency of the 
cavitation process in conditions of minimizing energy costs have been determined. 


CONCLUSIONS TO SECTION 6 


1. The estimation and methods of research of parameters of acoustic treatment of techno- 
logical media have been carried out. The functional relationships between the acoustic parameters 
of the cavitation apparatus and the rheological properties of the processing technological media 
have been revealed. The process of staged acoustic treatment of technological media is described. 

2. The analytical dependences of the determination of the main parameters of the cavitation 
apparatus have been obtained and the regularities of the processes of staged acoustic treatment 
of technological media in conditions of energy minimization have been established. 

3. A number of criteria and key parameters have been determined, the use of which is carried 
out as an assessment in the calculation algorithm, depending on certain known initial data of the 
cavitator and the medium. The values of the input resistance of the compensator have been deter- 
mined and the condition for the maximum transfer of energy from the acoustic apparatus to the 
technological medium has been obtained. 

4. An algorithm is proposed to substantiate the rational choice of the «acoustic apparatus — 
technological medium» structural-parametric system, which allows to vary not only the initial pa- 
rameters and the layout of the cavitator relative to the processing medium, but also to determine 
the effect of variable parameters on the maximum value of a particular criterion. 
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7 STUDY OF RELIABILITY OF TECHNICAL SYSTEMS RELIABILITY 


|. Nazarenko, M. Delembovskyi, 0. Dedov, 
A. Onyshchenko, |. Rogovskii, M. Nazarenko, |. Zalisko 
ABSTRACT 


The main states of reliability of a technical system and its elements are investigated by the 
example of vibration platforms for compacting concrete mixtures. Experimental studies on the de- 
velopment of vibration platforms for failure have been carried out. Complex indicators were used, 
which were the coefficient of technical utilization, the coefficient of availability and durability. The 
received data of malfunctions was made by fixing by groups of prefabricated units, parts and ele- 
ments to determine the data on their operating time. Based on these data, the analysis of the ope- 
rating time of the main elements to failure and the most frequently out of order was carried out. 
The conducted research identified the main prefabricated units and failing parts: engine, gearbox, 
synchronizer, vibration exciter, propeller shafts, couplings. At the same time, propeller shafts and 
couplings most often failed. In some cases, the destruction of bearings in vibration exciters has been 
evidenced. The parameters of the Weibull distribution law have been determined and graphs have 
been constructed for the model of reliability and failures of propeller shafts of vibration platforms. 
A graph of the dependence of reliability on the operating time and graphs of the regularities of the 
distribution of resource indicators and the distribution function, which served as information for the 
development of recommendations, were built. 


KEYWORDS 


Technical system, reliability, vibration platform, prefabricated units, complex indicators, distri- 
bution laws, technical utilization factors, availability and durability. 


7.1 STATEMENT OF THE PROBLEM AND ANALYSIS OF METHODS FOR DETERMINING 
THE RELIABILITY INDICATORS OF VIBRATION PLATFORMS 


Ensuring the reliability of vibration platforms in the construction industry [1-3] is an im- 
portant problem of increasing their efficiency. The process of operating technical systems is 
impossible without performing technical services, the most important task is to restore the 
working condition of machines or units by eliminating failures during operation, as well as 
overhaul or refurbishment at manufacturing enterprises. Repair planning assumes the avail- 
ability of information on reliability indicators in order to perform high-quality repairs with the 
least amount of time. 
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In vibration platforms, when performing work, there is a probability of failure (random and 
fluctuations) in the process of deviation of specified fluctuations, considered as a result of changes 
in time of the parameters of the system (amplitude of fluctuations, accelerations, permissible 
stresses in the structural elements of vibration platforms). The analysis of random events of oscil- 
latory processes [1, 2] makes it possible to assess such deviations associated with the processes 
of accumulation of failures, which are the subject of the theory of reliability. 

Assessment of the reliability of technical systems generally implies an integrated approach 
to solving the procedure for assessing reliability indicators at all stages of the existence of such 
systems. Methods for assessing such indicators are based on the use of mathematical models of 
the existence of a system, depending on the stage of the life cycle of a technical object. The physical 
essence of the existence of a system is due to the presence of a number of factors and factors 
that must be taken into account when creating a model. 

So at the design stage, methods are used to predict reliability characteristics based on general 
methods of probability theory, analytical mechanics, theory of machines and devices, resistance of 
materials, vibration theory, etc. 

structural methods of calculation are characteristic at the stages of manufacturing and opera- 
tion of technical structures, since they involve the consideration of the system by means of a struc- 
tural-functional diagram. In such a logical scheme, the transition of a system from one state to 
another, the relationship between such states, the interaction of system elements, taking into ac- 
count their functional purpose, manufacturing technologies, and conditions of use are investigated. 

Physical methods are used to determine reliability, durability and safety of technical objects for 
which degradation mechanisms are known under the influence of external and internal factors, which 
leads to failures — the limiting state during operation. Such methods are based on the representation 
of processes using mathematical models capable of determining reliability indicators taking into ac- 
count the design features, manufacturing technology of its individual elements, operating modes, etc. 

Each of the methods has its own advantages and disadvantages, therefore, in specific cases, it 
is more efficient to use a combination of such methods to solve the problem. This approach allows 
one to obtain a result that has not only an imaginary understanding of the existence of a technical 
system, but also has a physical essence, confirmed by experimental research. 

One of the methods that satisfy an integrated approach to solving the problem is the probabilis- 
tic-physical method of reliability analysis. This method of reliability analysis involves the use of probabilis- 
tic-physical reliability models, the parameters of which have a specific physical interpretation and relate 
the value of the probability of failure and the physical (determining) parameter, which leads to failure. 
Among the advantages of this method, there is a quantitative assessment of the reliability indicators of 
systems (elements) both on the basis of statistical information about failures and statistical information 
about physical processes leading to failures. In this case, all data is used in aggregate and separately. 

Thus, the application of the above method will make it possible to assess the necessary reliabi- 
lity indicators, including the complete characteristic — the distribution law based on the processing 
of statistical data on the failure of technical structures. 
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Prediction of the reliability of vibration platforms occurs at all stages of design, operation 
and repair [8-5]. At the same time, it is important to obtain the design parameters of parts and 
assemblies depending on predicting reliability indicators: mean time between failures, reliability, 
durability and maintainability. 

Today in this area the works of such scientists as: |. Kravchenko, V. Bolotin, A. Barash, A. Pro- 
nikov, K. Frolov and others are known. As a basis for research for assessing the reliability of 
vibration platforms, a model was adopted in which the object under study is in one of the following 
states: serviceable, faulty, efficient, disabled and extreme. The structural and logical diagram of 
the technical system and the path of transition from one state to another is shown in Fig. 7.1. 


Working condition 


{F] 
Operable condition 


Inoperable condition 
Boundary condition 





O Fig. 7.1 Scheme of the main states of object reliability: 1 — damage; 
2 — failure; 3 — transition to the limiting state; 4 — recovery; 5 — repair 


For the implementation of research, the following states of the system and their evaluation 
criteria were adopted: 

— working condition — a property of the vibration platform, in which it meets all the require- 
ments of the normative, technical and design documentation: 

— faulty condition — a property of a vibration platform, in which it does not meet at least one 
requirement of regulatory, technical and design documentation; 

— operable condition — a property of a vibration platform, at which the value of all parameters, 
characterizing the properties of performing specified functions, meet the requirements of norma- 
tive, technical and design documentation; 
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— condition — a property of a vibration platform, in which the value of at least one parameter 
characterizing the property of performing a given function does not meet the requirements of 
regulatory, technical and design documentation; 

— boundary condition — a property of a vibration platform [5], in which its further use for its 
intended purpose is unacceptable and impractical, or the restoration of its serviceable or service- 
able state is impossible or impractical (Table 7.1). 


@ Table 7.1 Criteria for the limiting state of the components of the vibration platform 
Unit or unit name Limit state criteria 


Frame structures 1. Tiring transverse cracks along the perimeter of more than 40 %. 
2. Concavity, twisting and other damage exceeding permissible limits, for elimination of 
which it is necessary to completely dismantle the frame structure 


Gear transmission 1. Limiting wear of teeth. 
2. Crumbling of the working surface of the tooth with a total area of more than 25 % 


Clutch coupling Worn splines or bearing seat on the shaft or breakage of the shaft, in which it is 
necessary to replace it 


Propeller shaft 1. Crack in pipes, forks along welded seams. 
2. Limiting wear of spline joints. 
3. Limit wear of the bores for the bearing shells 


Failures are divided by reason of occurrence, nature of manifestation, interrelation, groups of 
complexity and method of detection. In addition, there are resource and degradation failures. The 
classification of failures depending on the nature of their occurrence and the course of development 
processes leading to the onset of failure are shown in Fig. 7.2. 


Due to the By the nature es By difficulty By detection 
of the detection | [] By relationship method 


independent 
Dependent Hide 


Operational Displacin Third 








O Fig. 7.2 Classification of failures 


Thus, improving the quality and reliability of a vibration platform [5] depends on many reasons, 
the main of which is the quality of vibration equipment. An integral part of quality is the accuracy 
class of the vibration platform or the error with which the vibrations are simulated. Therefore, 
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the task of the study is not only to assess the reliability of existing structures, but also in new 
technical systems of complex structure. To solve such a problem, a set of criteria is proposed, the 
assessment of which will improve the quality of such technical systems. 


7.2 ENSURING THE RELIABILITY OF VIBRATION PLATFORMS IN THE CONSTRUCTION INDUSTRY, 
TAKING INTO ACCOUNT THE METHODS OF ANALYSIS 


Vibration platforms are widely used in the construction industry in the manufacture of concrete 
and reinforced concrete products. The effectiveness of their work largely depends on a fairly speci- 
fic accounting of the operating forces of the system and the reliability of the parts of the vibration 
platforms. Currently, there are practically no works on the development of methods for assessing 
the reliability of vibration platforms in various modes of their operation. However, the practice of 
operating vibration platforms in the construction industry testifies to the frequent failure of parts 
and assemblies of these machines, which significantly reduces their performance in the operating 
modes provided for by the technology. As a result of this situation, formed goods may be defective. 
Therefore, the problem of reliability methods, determination of its criteria and improvement of the 
design of elements of vibration areas is urgent. 

In accordance with DSTU 2861-94, reliability analysis for vibration platforms can be carried 
out in order to: 

— checking the feasibility of the established requirements and assessing the likelihood of achiev- 
ing the required level of reliability of the components and the object as a whole: 

— verification of the effectiveness of the proposed (implemented) measures for finalizing the 
design, manufacturing technology, maintenance and repair strategies to improve reliability; 

— predicting reliability and choosing the best ways to ensure or improve reliability. 

Considering the above, the methods of analyzing the reliability of vibration platforms [1—6] are 
used to predict the reliability, maintainability, readiness and measures to guarantee safety, as well 
as to fulfill the consequences of forecasting with the given requirements. 

Consequently, the tasks of reliability analysis and their scope depend on the stage of the life 
cycle of vibration platforms, the depth of development and reliability, the consequences of failures 
and their limiting states, and other factors. 

The DSTU 2861-94 standard establishes two main approaches to the analysis of the reliability 
of all objects, including vibration platforms: 

— analysis of the reliability of vibration platforms based on the results of measures and methods 
of ensuring reliability at the stages of design, production and operation in accordance with the 
reliability assurance program: 

— quantitative methods for analyzing the reliability of vibration platforms based on the analysis 
of operating conditions, causes and mechanisms of failures, reliability indicators of elements, main- 
tenance and repair strategies. 
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The analysis of the reliability of vibration platforms presupposes a reliability assurance 
program (Fig. 7.3), according to the results of the implementation of which the desired result 
can be achieved. 


Obtain a qualitative Predict and optimize reliability 
assessment of the achieved in accordance with 
level of the vibration ie established requirements, 
platform and establish the / RE Y \ maintenance and repair 
possibility of meeting | PROGRAM | strategies, and measures 
reliability requirements taken to irnprove reliability 


Prepare the necessary data Determine quantitative Carry out a comparative 
for the use of quantitative indicators of reliability and analysis of the reliability of 
methods for analyzing the establish the ability to meet the design options for 


reliability of a vibration the requirements of the vibration platforms 
platform 





O Fig 7.3 Block diagram of the reliability assurance program 


The analysis of the reliability of the vibration platform system is carried out using the fol- 
lowing methods: 

— qualitative analysis (determination of types of malfunctions, failure mechanisms of elements 
and their consequences for the system, analysis of the functional diagram of the system, analysis 
of the maintenance and repair system, construction of structural diagrams of system reliability): 

— quantitative analysis (construction of mathematical models of the reliability of elements and 
systems according to the indicators under consideration, obtaining quantitative indicators of reli- 
ability by calculation or modeling, analysis of the importance of failures and sensitivity, assessment 
of the possibility of improving the characteristics of the system based on backup subsystems and 
maintenance and repair strategies). 


O Fig. 7.4 Diagram of the stages of the methodology 
for analyzing the reliability of vibration machines 
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Accordingly, the assessment of the analysis results involves comparison with the mandatory 
reliability indicators and recommendations for measures to comply with the mandatory reliability 
indicators, which may include a design review, identification of weak points, imbalances, modes, 
replacement of parts with a high risk of malfunctions, development of alternative ways to improve 
reliability, implementation trade-off analysis and cost estimation of vibration platform design op- 
tions as a whole. 


7.3 EXPERIMENTAL STUDIES OF RELIABILITY ASSESSMENT BASED ON FUZZY LOGIC 


The solution to the problem of increasing the reliability and efficiency of vibration machines is 
achieved by studying the parameters and characteristics of failures of the elements of a vibration 
machine in the operation stage to develop recommendations for determining the reliability indica- 
tors of units and parts of vibration machines with an increase in their efficiency. 

To a large extent, the study of reliability indicators requires the involvement of a large number 
of observations, which in fact takes quite a significant amount of time. Taking this into account, it 
is necessary to carry out studies of reliability indicators using modern information technologies and 
alternative solutions to the problem. 

When performing work on modeling the assessment of the reliability of vibration machines 
based on a fuzzy model, it is necessary to represent them in the form of fuzzy networks, the ele- 
ments and sets of elements of which implement various components of fuzzy models and stages 
of fuzzy inference. 

In his work, Lotfi Asker Zadeh presented a fuzzy productive model as follows [4, 7]: 


(): Q; P; A= B: S: FN, (7.1) 


where Q — field of application of fuzzy products; P — condition for the activation of the fuzzy produc- 
tion core; A — condition of the kernel (antecedent); B — kernel output (consequent); S — method 
for determining the quantitative value of the degree of truth of the kernel output; F — coefficient of 
confidence of fuzzy products; N — post condition of the productive rules. 

A fuzzy imitation of the relationship between antecedent and consequent can be viewed as 
fuzzy production: 

if xe A, then ye B; 

X — domain of definition of the antecedent; 

A — fuzzy set visible on X: 

uu, (E) <[ 0,1] -membership function of an odd set A: 

Y — area of definition of the consequent: 

B — fuzzy set defined on Y; 

1, (E) [0,1] — membership function of an odd set B. 
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If the ratio function of odd sets A—j1,(E), is known, then for odd sets B, the ratio function is 
determined by the composition rule: 


U,(C) = SUp{T (Wa(X)sthe(XY)) (7.2) 


where sup — operation of determining the upper limit of the sets of parts; T — T-norm operation. 
When modeling the failures of vibration machines as a rule for calculating fuzzy implication, the 
classical fuzzy implication of L. Zadeh: 


I, (x,C) = max min[ i, (x), us (y)}.[4-bs (x) |}. (7.3) 


For the purpose of inference, the main methods in a fuzzy production model are forward and 
backward inference. Direct inference is based on such a concept as fuzzy modus ponens. 

In the process of developing a fuzzy production model, the assessment of failures of vibration 
machines requires the formulation of certain conditions X ={E,},i = An - factors that are the 
source of failures, and therefore the entire range of conclusions Y = {Vy = 1,m is an indicator 
of failures of various elements and components of vibration machines. 

Based on the previously established objects and concepts of the systems under consideration, 
namely vibration platforms, it is necessary to formulate the corresponding linguistic variables: 

M — electric motor; 

S — synchronizer; 

C — clutch: 

R — bearing (roller); 

P — propeller shaft; 

SP — support; 

F — frame (vibration plate). 

Let's take B as the initial linguistic variable — the level of failure. 

Let’s define terms for input and output linguistic variables. For input linguistic vari- 
ables P, R, M let's introduce term-sets {critically low, low, average, high, critically high}, 
for S, C, SP, F — {low, medium, high}. The term sets of the initial linguistic variable describe the 
failure detection rate — {absent, unlikely, low, high enough, critically high}. 

The level of failure of a vibration platform is estimated in the process of fuzzy inference, using 
a large number of fuzzy rules that make up the aggregate knowledge base of a given subject area, 
represented in the form: 


If (x, is A) and (x, is B), then Cy is C), 


where A, B, C — linguistic values, identified in a fuzzy way through the corresponding membership 
functions for the variables x,, x, and y. 
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As a result of the study, a fuzzy knowledge base was determined, presented in the form of 
certain rules [7]. 

For further modeling of this kind of system, the Fuzzy Logic Designer software toolkit is used, 
which is a MATLAB b2020 extension package containing tools for designing a fuzzy logic system [6]. 

When determining the shapes of the curves of membership functions, it is taken into account 
that they are constructed subjectively according to the results of a survey of experts, therefore, 
in a sense, they are «approximate». The most common membership function is of the Gaussian 
type [8, 9], but in practice, the shape of the membership function curves is chosen based on the 
complexity of the calculations. The most widespread are triangular and trapezoidal membership 
functions due to their versatility and lower requirements for computing resources in their hardware 
implementation. 

To display the selected fuzzy subsets of linguistic variables, let's them. The parameters of the 
input membership functions are shown in Fig. 7.9. 


03 04 05 05 O7 
input variable "Heuryn™ 


Cepequa 


on 


03 04 05 06 07 Os 
input variable *Pana(Bid poryi6a)" 





O Fig. 7.5 Parameters of input linguistic variables 


The Mamdani algorithm, which has received the greatest practical application in fuzzy modeling 
problems and consists in the application of a mini-maximum composition of fuzzy sets, is used as 
a mechanism for the fuzzy logical inference being developed. 
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The process of processing fuzzy inference rules in this case consists of four stages: 

1. Fuzzification, which consists in determining the degree of truth, that is, the value of the 
membership function for the premises (left-hand sides) of each rule. 

2. Fuzzy inference, which consists in applying to the conclusions (right side) of the rules of 
the calculated truth value for the premises of each rule. As the inference rules in the Mamdani 
method, the minimum (min) operation is used, which «cuts off» the membership function of the rule 
conclusion by height, which corresponds to the calculated degree of truth of the rule prerequisite. 

3. A composition that combines, using the maximum (max) operation, all fuzzy subsets defined 
for each variable output, and forms one fuzzy subset for each variable output. 

4. Defasification, which realizes the scalarization of the composition result, i.e. transition from 
fuzzy subset to scalar value. 

The scheme of the described fuzzy inference system implemented in the MATLAB medium 
is shown in Fig. 7.6. 


< Fuzzy Logic Designer: JocnigKenna Haaitnocti - oO x 







File Edit View 





Reliability study 


> (mamdani) 


Frame (vibration plate) 








Aggregation 


Defuzzification 


O Fig. 7.6 Scheme of a fuzzy inference system in MATLAB 
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This system makes it possible to assess the failures of vibration platforms based on the input 
linguistic variable, the parameters of which are shown in Fig. 7.7. 





O Fig. 7.7 Output linguistic variable parameters 


When performing the appropriate calculation, it is possible to preliminarily assume based on 
certain studies and observations that have input data, namely: 

1. Engine — probability of failure 0.2. 

2. Bearings — probability of failure 0.7. 

3. Propeller shafts — probability of failure 0.7. 

4. Frame (vibration plate) — probability of failure is 0.1. 

9. Clutch — probability of failure is 0.4. 

6. Synchronizer — probability of failure 0.25. 

Thus, taking into account the corresponding coefficients, let’s obtain the results of construct- 
ing fuzzy logic, which in this case will be based on certain 27 rules, according to which the overall 
failure rate of the vibration platform will be 0.624 (Fig. 7.8). 

Accordingly, based on this range of values, it can be argued that the level of failure of the 
vibration platform is at a sufficiently high level. 

Conducted similar studies of the presented fuzzy model for different sets of initial data also 
showed acceptable results (Fig. 7.9-7.12). 

Failure to operate the vibration platform leads to downtime of the entire technological line and 
damage to the concrete mixture that has already been supplied to the line. In turn, this approach 
as a whole gives a clearer idea of which constituent elements need to be paid more attention to 
during the repair and maintenance of vibration platforms. 

The use of this approach can significantly reduce the cost of maintaining vibration platforms 
and, accordingly, reduce the downtime rate of the production lines of the shop. 
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O Fig. 7.8 Results of calculating fuzzy logic for 27 rules 
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O Fig. 7.9 Indicators of the effect of motor and bearings on failure 
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O Fig. 7.11 Indicators of the effect of propeller shafts on failure 
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O Fig. 7.12 Impact of bearing failure 











7.4 CALCULATION OF VIBRATION RELIABILITY OF A TECHNICAL SYSTEM 


When calculating the vibration reliability of optimal vibration protection, a vibration platform [7] 
was considered as a system described by the equation of motion: 


i+ 2e-U+05-u = —a, lt). 
Accordingly, the quality conditions are reflected in the form: 
lutt)| < uz, jale]| <a., (7.4) 


where a= a, +U — absolute vibration acceleration. As a criterion for the optimality of the system 
in terms of reliability, the condition is selected so that the reliability functions (7.4) at the time 
instant t = 7. acquire a maximum value. 


sup 


P(t)=P, * 
Vm Vee: te[0t| 








ax, D|<a.,~ |alx,,2)|-<u. | (7.5) 
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Taking into account the high reliability of systems for stationary oscillations, this condition, 
according to: 


t 
Pi=1- } Adc, 
0 


is equivalent to the condition A(T’) = min, €, w,, where the average number of emissions per unit 
time is calculated for a rectangular region of available states (7.9). 

It is proposed that the process a,l(t) at the output of the system is narrow-linear with 
a carrier frequency @. In Fig. 7.5 the results of the corresponding calculation data A(T) are carried 
out at u. = 10m, a, = 20 m/s. 

In Fig. 7.13 a graph is built according to the calculation A(I) performed for the case of an 
exponentially correlated process at the input with a correlation Ke (t) = 0% -exp(—a|t|). Con- 
sidering the above, calculations have been made for the parameter values m e€/@, = 0.5, a, = 20,, 
U. =(0,/a°)-10~. Curves 1 in Fig. 7.5, 7.6 correspond to band emissions |a|<a,, curve 2 — 
band emissions U| < U,, curves 3 are plotted for the total number of emissions. 


-4 
100005. 0.10 0.15 0,0 






O Fig. 7.13 Optimization in terms of reliability with 
a narrow-linear initial impact 


A 
a 
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O Fig. 7.14 Optimization for reliability with an exponentially 
correlated initial action 
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Basically, vibration loading is accompanied by a process of accumulating fatigue damage. The cal- 
culation of damages characterized by a one-dimensional process v(t) was carried out accordingly. 
Accordingly, this process is rather slow, therefore, it requires a transition from v,,, =v, +flv,,s,] 
to a continuous analogue — a kinematic level: 

m = glv,s). (7.6) 

dt 

In such a case, the measure of damage v(t) is an unfortunate function, 0<v <1. In the above 
equation (7.5) s(t) — some characteristic value for the load process f(t), for example, amplitude. 
Taking this into account, the quality condition was chosen in the form v <v, =1, $< 8,(v), where 
s.{v) — the limiting value of the bearing capacity, taking into account the accumulation of damage. 
In this case, the reliability function: 


Plt) = Piv(t) < 1st) <s.),r€[0,¢ |}. 


As arule, as the damage caused by the fatigue process takes place: 


where N(S-) — the number of cycles leading to failure with a characteristic S,, then the condition 
of fatigued failure takes the form: 


1 


~ NS) 


To assess the durability, a value 7, equal to the duration of the load, which is necessary to 
achieve the condition v =v, = 1, is used. Therefore, with a low dispersion of durability, the value 7, 
is close to the mathematical expectation (T). The value 7. is called the expected durability. Accord- 
ingly, if do not take into account the probability of instantaneous destruction due to the process 
bursting beyond the level s.{v), then the expression for the expected durability takes the form: 


so —' 
=7,|(2=| 77 
* |« Ms) 


where 7, — some effective cycle period; pls) — distribution of the parameter s — the loading 
process f(t). 

Having considered the suitable processes, it can be argued that for stationary normal narrow- 
linear load processes, the maximum value m s(t) =f__ (t) is taken as the characteristic value m s(t). 


max 


In this case, the distribution of the maximums of the normal process was found by the formula: 
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: 2(B° — 
v.-4 (v, — al’ Vv. -a 
+ J2n -exp| — =— |: P|] ————_ ], (7.8) 
0, 20, Be > lo, 
where 
B- 0; 0, Vivi 
o = ALO 


Calculations made according to the formula (7.7) with a power-law approximation of the fa- 
tigue curve: 


N=N, (2 iss], 
S 
N >, (s<s,) 
give this expression: 
N,-X"-T, 


_j—____— = ____: 
ane r(t4 2) put m+2) 


where X = s,/o,; N,, 5, m — empirical constants; P(X*,7) m function X* — Pearson distribution. 
The same expression, expressed due to an incomplete gamma function Ta, x), is: 


N,-X"-T, 
xe): 
ome Tl 4 a 
Gea 


Accordingly, for wide-linear load processes, there are different ways to select the load pa- 
rameter, leading to different estimates of durability. To give an example of a cycle, the following 
maxima, maximum values, half-difference of subsequent maxima and minima, maximum process va- 
lues at a certain characteristic time interval, etc. are considered. To schematize a given process by 
a load of some equivalent narrow-linear process, methods have been proposed: outliers (the num- 
ber of cycles coincides with the average number of process zeros), maxima (the number of cycles 
is taken as the average number of process maxima), swing (a cycle is characterized by an amplitude 
equal to half the process increments between adjacent ones), full cycles (a method consisting in 
the sequential elimination of the process of intermediate cycles with increasing amplitudes), etc. 


= 
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Having considered the expression for the probability density p(s) of the load parameter s(t) in the 
case of a stationary normal process f(t), the distribution of the maxima p,s) is given by formu- 
la (7.8). By emission method: 


s a 
Pts)=—exp| -2- |: 
(s Sen a 


s 


by calculation method: 


2 2 ae 
pis) =P ‘-e(-P =} 
fe) 








2 
S 


by full cycles method: 


2g 2 92.62 
FF eof (s<s,J; 


0; 
© > exp _S (s<s,) 
B of 20° ~ 


where 8 — coefficient of wide linearity of the process; a, c, s, — some constants. 
The nature of the dependence of the expected durability (7.7) on the nature of the wide 
linearity B is shown in Fig. 7.7. 





10 15 20 25 30 3.4 


O Fig. 7.15 Comparison of estimates of expected durability when schematizing 
by different methods: O — bursting; 1 — maximums; 2 — swings; 3 — full cycles 


The expected durability corresponding to the probability density Ps) is indicated 7, accordingly. 
The calculations took into account that s, =0,m=9. 

The maximum method gives a lower estimate for the expected durability, the swing and full 
cycles method overestimates this estimate. For narrow-linear processes, all of the above methods 
lead to practically the same results. 
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7.9 JUSTIFICATION AND APPLICATION OF THE PROVISIONS OF THE MATHEMATICAL THEORY 
OF RELIABILITY 


Let's consider the exploitation of a number of N elements over time ¢. Let there be NV, able-bodied 
and n failed elements by the end of the service life. In this case, the relative number of failures Q(t): 


Q(t)=—. (7.9) 


If N is large enough, it can be viewed Q(t) as the probability of failure. 


Then the probability of failure-free operation is P(t): 
N n 
P(t)=—- =1-—— =1-Q(£). (7.10) 
(=e 1-21-00 


It follows from (7.10) that: 
P(t)+ Q(t) =1. (7.11) 


The distribution of conditions in time is characterized by the distribution density function f(t) of 
the operating time to failure: 


dQ(t) 1dn_ aP(t) 


(js oe (7.12) 
dt =Ndt dt 


In this case, the probability of failures and failure-free operation in the density function f(t) is 
expressed by the dependencies: 


Q(t) =|F(t)de. (7.13) 
P(t) =1—Q(t) =1- | f(t)de = [F(t)de. (7.14) 


Failure rate A(t): 





f(t aP(t aP(t 
Me) =F = ff) (7.15) 
Let's rewrite (7.15) as: 
-A(e)dr= Fe) (7.16) 
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Integrating (7.16), let's obtain: 


t a 
InP(t)=-[A(t)dt or P(t)=e” (7.17) 
0 


During the period of normal operation, failures have a constant intensity, since worn out (gra- 
dual) failures do not appear yet, but only sudden failures occur, i.e. 


A(t) =A =—=const, (7.18) 
M, 
where m, — mean time to failure, hours. 
Probability of failure-free operation: 


t 
— fade 


P(t)=e*® =e”, (7.19) 


It obeys the exponential law of distribution of failure-free operation t and is the same for any 
equal period of time (Fig. 7.16). 





P(t), A(t), fl) 





O Fig. 7.16 Function of probability P(t) of failure-free operation, 
probability density f(t) and failure rate A(t) according to exponential distribution 
during normal operation 


The significant advantage of the exponential distribution is its simplicity: it has only one parameter. 
Distribution density in general: 


f(t) =- at 0). ie (7.20) 


7.6 RESEARCH OF INDICATORS OF RELIABILITY OF TECHNICAL SYSTEMS 


Experimental studies were carried out on the basis of the collected information on the deve- 
lopment of vibration platforms to failure [1, 3]. The received data of malfunctions was carried out 
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by fixing by groups of prefabricated units, parts and elements, in order to separately determine the 
data on their operating time. Based on these data, an analysis of the operating time of the main 
elements between failures and the most frequent failures was carried out. 

The studies carried out [2, 7, 10-19] identified the main prefabricated units and parts that 
fail: engine, gearbox, synchronizer, vibration exciter, propeller shafts, couplings, photographs of 
such failures are shown in Fig. 7.17. 





O Fig. 7.17 Photo fixation of failures of mechanical assemblies and parts of vibration platforms: 
a — failure of the coupling; b — failure of the electric motor; c — failure of the bearing assembly; 
d — failure of the vibration block; e — failure of the propeller shaft; f — failure of the synchronizer 


According to the results of the study, it was found that most often propeller shafts and 
couplings failed. Destruction of bearings in vibration exciters was recorded only in a few cases. 

In the investigated technical objects, depending on the design, the working body can consist of 
both separate elements such as a vibration plate or a form, and of a welded frame. In the course of 
research, in this group of elements, frequent defects were found as a result of contact of vibration 
plates with a form (Fig. 7.18). 

Based on the results of the performed structural examinations, the main elements were 
established (Fig. 7.19), the failure of which has a direct impact on the reliability and safety of 
vibration machines. 

The research methodology for the operating time of parts to failure consisted in conducting an 
experiment in real operating conditions of VB-10V vibration platforms. 

Due to the frequent fixation of the identified failures of the propeller shafts of vibration 
platforms, the results of which are shown in Table 7.2. The construction of a statistical series 
for studying the distribution patterns of continuous random variables is possible with a sufficient 
number of observations. In this case (if N>30) the values obtained during observation in the 
variation series were grouped, and the entire series of values of the observed trait from ¢, to ¢, 
was divided into ordinary intervals. 
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O Fig. 7.18 Photo fixation of failures of the working bodies of vibration platforms: a — failure 
of vibration plates; b — deformation of the limiters; c — damage to the contact platform; d — damage 
to the protective casing; e — damage to vibration plates; f — blockage of the vibration platform 


The number of intervals was determined from the condition of detecting the regularity of the 
distribution of the indicator values depending on the sample size WN. 
The number of intervals was determined by the formula: 


K=N, (7 21) 


where N — the sample size. 

For the convenience of further calculations, all intervals are assumed to be equal, targeted and 
without gaps. However, in some cases, when processing statistical data that are rather unevenly dis- 
tributed, it is sometimes convenient in the area of the highest distribution density to choose narrower 
intervals than in the area of the smallest. The value of the interval was determined by the formula: 


p= fo Io, (7.22) 
K —1 
where 7, and 7_.. are the highest and lowest values of reliability indicators in the summary table 


of information. 
The left and right boundaries of the distribution area are shifted by 0.9h and taken accordingly: 


t, =T,,, -0.5h, (7.23) 


min 


t, =T,,. +0.5h, (7.24) 


where €, — the left border of the first interval; ¢, — the right border of the last interval. 
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O Fig. 7.19 Main structural elements of vibration 
platforms affecting reliability 
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@ Table 7.2 Operating time to failure of propeller shafts 


Operating time to failure, h 


125 148 126 131 127 
130 159 135 164 133 
150 15/7 145 166 141 
162 178 175 2 163 
167 184 195 200 182 
186 188 202 242 23g 
199 219 209 202 242 
pee 226 268 261 271 
225 236 272 2/9 288 
28/7 269 280 288 290 


In subsequent calculations, the average value of the interval was taken: 


eens (7.25) 


avi 2 


where t, — value of the right limit of the n-interval. 

After determining the boundaries of the intervals of the grouping, the frequencies /, the number 
of observations that fell into each of the intervals [h], i=1, K. 

The sum of the frequencies of all intervals is equal to the sample size: 


k 
Yn, =N =50. (7.26) 
i=1 


Another characteristic of the statistical distribution is the frequency P, (relative frequency), which 
is determined for each interval by the ratio of frequency n, to the total number of observations NV: 


pa, (7.27) 
N 


Empirical distribution density of resources F(t) in the i-th interval / = 1,2,3... 


F(t) = —_. (7.28) 
| N-h, 
Empirical function: 


k ~~ 
Alt=P +P; F=>P. (7.29) 
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The average value of the resource 7..., hours was found by the formula: 


r.av.3 


T= (7.30) 





(7.31) 





The results of processing observational data using the sum method are summarized 


in Tables 7.3, 7.4. 


@ Table 7.3 Numerical values of distribution parameters 


Interval 


Interval Frequen- ,» 


silaee Interval limits mid point f cy ny . Ty; Tn, 

1 111.25-1388.75 125 / 1(9201.6 971.25 134 761.2 
2 138.75-166.25 152.5 10 a7iea Je 1 662.5 D7)e ase} 

d 166.25-193.75 180 8 g7.939.1 1 590 300 312.8 
4 Wavio—ee lea peuy a 6 Mees | Beye 208 7096 
5 221.20-248.75 235 7 61876.6 1741.25 433 136.2 
6 248.75-276.25 525 7 eee oda 934 198.7 
7 276.25-303.75 580 is 92 264.1 1518.75 461 320.5 

ell ~=10705 Y=2 433 830 


@ Table 7.4 Calculation of characteristics by the method of sums 


Interval number = Interval mid pointt,,, Frequencyn, K,=110 K,=84 

1 125 ] 7 - 

2 152.5 10 = = 

3 180 8 33 — 

A 207.9 6 20 61 

Q 239 / 19 36 

6 920 / |2 ly 

] 080 a) 9 0 
le! ie 
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Auxiliary intervals are determined: 
u, =k, -A,. (7.32) 


The average value of the resource T 


r.av.3 


hours, is found by the formula: 


_ Mey 


ls av.i , (7.33) 
: “ON 


The mean square deviation of the resource o, taking into account the auxiliary intervals: 


= 2 
se pte (7.34) 


Variation coefficient: 


V=—. (7.35) 


Actual value of the criterion A for extreme points: 


d,=—lt—t,,), f= 4 (7.36) 


p 
t 


where t and i, are adjacent information points. In this case, the condition A, <A, is fulfilled. 
Thus, as a result of checking the information for outliers, it was found that the minimum and 
maximum values of the sample are ¢,,,=129 and ¢,,..=126 with probability y = 0.95. 


The histogram of relative frequencies (Fig. 7.20) indicates that in this case the process is 
subject to the Weibull distribution law: 


b-' b 
f(t) = (2) eol-(') (7.37) 
a\a a 
where a and b are distribution parameters a> 0; be O. 


It was found that Pearson’s criterion y° with a large number of observations reduces errors 
to a minimum, which compares favorably with other goodness-of-fit criteria: 


«| Fle.) - FC) | 
Te Heh 


2 
(7.38) 


where N — the total number of observations; fh, — width of the /th interval; k — the number of 
staking intervals; f(t,), f(t,) — respectively empirical and theoretical probability density of distribution; 
t.,,— midpoint of the /-th interval for a continuous random variable or possible values for a discrete one. 
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—~— f(t)t. —i fltle. 




















0 
138.75 166.25 193.75 221.25 248.75 276.25 303.75 Tr, year 





O Fig. 7.20 Concordance plot of empirical and theoretical distribution 


Therefore, it was found that P{r,x°) = 0.99 exceeding the accepted significance level of 0.95. This 
gives grounds to assert that the hypothesis that the empirical distribution belongs to cannot be rejected: 


P(t) = vol-(5)} (7.39) 


The Weibull distribution law is characterized by the dependences: 
t b 
F(t) =a(t)=1- eol-(*) 7.40) 


In relation to these calculations, the possibility of failure and trouble-free operation was ob- 
tained (Fig. 7.21). 


138.75 166.25 193.75 221.25 248.75 276.25 303.75 
O Fig. 7.21 Probability of failure and failure-free operation 
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7.7 DISCUSSION OF RESEARCH RESULTS 


The results of the studies carried out indicate a low level of reliability of the propeller shaft, 
which leads to long downtime during the unplanned repair time, which is 0.8—1.2 hours, if there 
are spare parts in stock. As a result of processing the statistical information of the propeller shaft 
reliability indicators, it was found that the average resource is 262 hours, the standard deviation 
is equal to Go = 900.2 hours, the coefficient of variation is V=0.27. 

A hypothesis is put forward about the error of the empirical distribution of resource indicators 
in comparison with the theoretical Weibull’s law. The parameters of the Weibull distribution law 
were determined and graphs were constructed for the model of reliability and failures of propeller 
shafts, the graphs of the regularities of the distribution of resource indicators P(t) and the distri- 
bution function F(t) used in the developed recommendations. 

Summarizing the presented calculation method, it is possible to build an algorithm for calculat- 
ing the determination of the reliability indicators of vibration platforms (7.22). 


CONCLUSIONS TO SECTION 7 


1. The main states of reliability of a technical system and its elements have been investigated 
by the example of vibration platforms for compacting concrete mixtures. 

2. The model for assessing the reliability has been developed based on fuzzy logic, implemented 
in the MATLAB b2020 medium. 

3. The obtained results of modeling, taking into account the normalized indicators of reliability, 
show a «sufficiently high» level of failure. 

4. When processing statistical data on the operating time to failure of the propeller shafts 
of vibration platforms, the Weibull distribution law has been established and the possibility of 
failure-free operation of such structural elements has been determined. 

9. The obtained results of experimental data processing can be used in the development of new 
designs of similar technical systems. 
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ABSTRACT 


The distribution of stresses and strains in the elements of vibration structures of machines 
has been investigated. Vibrograms of changes in the stresses of the forming surface were ob- 
tained for individual elements located near the application of a dynamic external force. Comparing 
the obtained results of the stress-strain state of the forming structure, a new effect of using 
high frequencies of the operating mode when compacting concrete mixtures was discovered. 
Due to the application of a spatial forced force to the forming structure, a complex stress-strain 
State of metal structures is created. And direct contact with the concrete mix helps to reduce 
energy consumption for the compaction process. The study of the motion of systems related 
to block structures, the adopted model, which is a two-mass vibration system. The identified 
transient process is intended to be taken into account when determining the parameters and 
locations of the vibrators. In such modes, the forms of natural vibrations of the system are 
realized with large vibration amplitudes and, accordingly, a lower frequency. And this opened up 
a real opportunity to reduce the energy consumption of vibration machine drives. The stress-strain 
state of the frame and forms of a vibration unit with spatial vibrations has been investigated. 
The distribution of stresses in the frame elements is uneven. The concentration of stresses in 
the welding places of elements having small values in comparison with the ultimate strength of 
steel has been determined. The static and dynamic loads of the slewing ring of a truck crane 
have been investigated. The positions of the system and its elements with the highest stresses 
are established. The results obtained were used in the design of the metal structures of the 
machines under study. 


KEYWORDS 


Metal structure, vibration system, stress-strain state, form-generating structure, model, nu- 
merical calculations, natural vibration modes, vibrator, spatial vibrations, truck crane, slewing ring. 


8.1 STATEMENT OF TASKS AND RESEARCH METHODS 


Studying the distribution of stresses and deformations in structural elements of vibration ma- 
chines will allow analyzing and identifying places of stress concentration, which allows optimization 
of design solutions. 
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Based on the existing methods and a certain mathematical description (see Section 2), the 
problem arises of specifying the constituent elements of the pre-static and dynamic action of the 
load determined in the study of the machine-medium systems. 

Thus, the formulation of the problem of computer modeling is reduced to the construction of 
a system that replaces and reflects the actual system, and in cognitive processes is in the corre- 
sponding similarity with it. 

In the general case, any existing or created object can be considered in three-dimensional 
space both under static and dynamic action on it [1, 2]. Since the entire system is three-dimen- 
sional [3, 4], then all elements belonging to it also have three coordinates of displacements and 
three angles of rotation relative to the axes, that is, six degrees of freedom (DOF). 

In accordance with the principles of modeling, the determination of the structural parameters 
of the elements is based on the ratio of the geometric dimensions and the type of arising (possibly 
arising) stresses in the structural elements. 

So, to simulate a structural element located only in the field of tensile — compression forces (pro- 
vided that it works in the zone of elastic deformations) and its length with respect to the cross section 
differs by an order of magnitude, it is enough to consider a one-dimensional problem. Since the pro- 
cesses occurring within the section can be neglected. The model that will reflect the behavior of such 
a system will look like a bar, the ends of which are clamped by a fixed and movable hinge, respectively. 

For structural elements that are in the field of complex types of forces and the ratio of overall 
dimensions to plane to height are within several orders of magnitude, a system of two-dimensional 
elements with a thickness — plates — should be considered. Accordingly, structural elements that 
have overall dimensions of the same order and (or) are in a complex stress-strain state are con- 
sidered volumetric elements and, accordingly, the system (model) is considered three-dimensional. 

Knowledge of the stresses in the sealing layers under the action of a force load allows one to 
determine the sealing zone and, as a result, to determine the parameters of the sealing equipment 
with sufficient reasoning. 

Improvement of the structures of vibration machines is possible by simulating the operation of 
this structure under loading, starting from the stage of an unloaded structure and ending with its 
full load [5, 6]. The study of displacements and deformations of structural elements of technolo- 
gical machines will allow identifying places of stress concentration and optimizing design solutions 
using finite element analysis [7]. 

At all nodes, generalized coordinates 2, are specified, called nodal offsets, the totality of which 
for a given element is written in the form of a matrix: 


T 


eae, (8.1) 


where N/ — total number of node displacements of the element; the sign 7 means the transposition of 
the matrix. Nodal displacements can represent the components of the vector of displacement of the 
nodes along the coordinate axes, as well as the angles of rotation of the element at the nodal points. 
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Within each element, for the components of the displacement vector of any point V/, an ap- 
proximation is set in terms of nodal displacements, which are unknown quantities: 


u,=®,(M)A,, i= 12,3, k= 12...) (8.2) 


| 


the same in matrix notation {u}={@${2 and vector form: 


i= 0,6, ={®} {A}. (8.3) 


where ®,,(M) values are called functions of the shape of the element and express the relationship 
between the nodal displacements and the displacement of the point M of the body; polynomials 
are usually used as a shape function; outside the element, these functions are set equal to zero; 
©, = 0,8, {6} ={,,6,,..,6,}, & — unit vectors. 

Relations (8.3) are substituted into the body equilibrium equation, from which the nodal dis- 
placements {a} for each element are determined. 

When using the finite element method, it is most convenient to obtain the body balance equa- 
tion based on the principle of possible displacements. Let wu is the field of displacements of points 
of a deformable body under the action of external loads applied to it. By setting each point a small 
displacement du allowed by the constraints imposed on the body (movement is possible). Accord- 
ing to this principle, the increase in the work of internal forces is equal to the work of external 
forces on possible displacements, i.e.: 


6U = dW. (8.4) 
Denoting as qg the external load distributed over the volume of the body V, and as p — the load 


distributed over its surface S. 
Let’s obtain: 


SW = [q-d0dV +) p- dudS. (8.5) 
V S 
The expression for the work of internal forces has the form: 


SU = J o- Sed V, (8.6) 
V 


where 6 = 0,6, — stress tensor, ¢ = €,6@, — strain tensor, é, — unit unit vectors, /, j=1, 2, 3. 
Then relation (8.4) takes the form: 


o-dedV = | q-dudV +] p- dudS. (8.7) 
fo sao 
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In the case of small deformations of the body: 


e= VU, (8.8) 
where 
 ¢O 
vi=—| 4S lee 
2 ax, ax, Ji! 


— tensor operator; i, /=1, 2, 3; X,,X.,X, — coordinate axes directed along the unit vectors 6,,é,,6,. 
Substituting (8.2) into (8.8), let's obtain an expression for the components of the strain 
tensor in terms of the nodal displacements: 


 a@, 
| Pee ep (8.9) 
OX, OX, 





or in matrix form: 


fel = {Bl al, (8.10) 


where 
1{ 0b, 0, 
B eee 
2} =(vs}=|3( 22a 
— matrix associates deformations with nodal displacements. The relationship between the compo- 
nents of stress and strain tensors for an elastic body is expressed by a law: 








Gy = DE y (8.11) 
where D,,, — the elastic constants of the body, i, j, k, /=1, 2, 3, or in matrix form {o} = {De}. 


Substituting expression (8.30), let’s find the dependence of the stress tensor on nodal dis- 
placements: 


fol={DMBMal. (8.12) 


After transformations, let’s obtain the equation of equilibrium of an elastic body containing the 
movement of its points: 


[DOVa-3(Va)dv = [q-dadv +) 5-805. (8.13) 
V V § 
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Let’s apply relation (8.13) to a finite element with some volume V, bounded by the surface S., 
and find: 


Ve 


Sh, [V®,-DV®,-2,dV - [GbdV - os} =0, (8.14) 
Ve 


where /,/ = 1,2,...,N. 

since 6A, — arbitrary nonzero values, the last equality requires all expressions in curly braces 
to rotate to zero. From these conditions, let's obtain a system of linear algebraic equations ex- 
pressing the equilibrium conditions for a finite element: 


Kabah, (8.15) 


where K, = = |V%,. DV® dV — stiffness matrix of the element, which, using relations (8.10) and 
(8.12), can also be written in the form: {K} = {By {D\{B}: 


=[9-OdV+] p-O4d5 
Ve 


— vector of nodal forces of the element, where /,/ = 1,2....,N. 

The set of equations (8.15) for all elements is supplemented by the equations of ligatures 
imposed on the body (limiting conditions), which are a system of equations for the equilibrium of 
the body under consideration, written in a form similar to (8.35): 


IK HAS = {Ff (8.16) 


where 1K} is called the global body stiffness matrix; 1a} and iF — vectors of nodal displace- 
ments and forces of the whole body. 

Equations of the type (8.16) are used for strength analysis of structures under static load. 
From their solution, the vector of nodal displacements is determined, then, according to rela- 
tion (8.2), one can find the displacement of the points of the body, and according to (8.8)-(8.11) 
or (8.13) — deformations and stresses. 

From (8.15) it is easy to obtain the equation of motion of the element. By introducing, ac- 
cording to the d'Alembert principle, the volumetric forces of inertia into the integral for the nodal 
forces (8.15): 


o> 
a a ee (8.17) 


let’s obtain a system of equations: 
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{MHA} +{K} {a} = {Fh, (8.18) 


where M, =| p@,- dV — mass matrix of the element; p — material density; A — the second 
time derivative of the vector of nodal displacements. In the presence of viscous resistance forces 
in the system, proportional to the speeds of the points, the matrix of damping coefficients {B} is 
introduced into (8.38), after which the equations of motion take the form: 


{MHA} +{B}Al+ {KT {A} = {F}. (8.19) 


When studying the problems of elastic resistance of structural elements, the equilibrium equa- 
tion is added taking into account the change in the geometry of the body in the deformed state: 


{K}-B{K,} {a}=o. (8.20) 


where by means of a matrix of geometric stiffness 4K, +, called differential in MSC.vNOW, the 
work of external forces due to a change in the geometry of the body is taken into account; 
8 — loading parameter. 

Equating to zero the determinant of the system (8.40) det| {kt — B{K,} | = Q, let’s find the 
value of the load parameter B,,f,,..., at which there are nontrivial displacements for nodal dis- 
placements {a}, i.e. new forms of body balance appear, different from the original. Such values B, 
called critical, show how many times the critical load F, at which the loss of resistance of the 
original form of body balance occurs, is greater than the current load F: 


F° =BF, i=12.... (8.21) 


The first lowest critical load F,” is usually of practical interest, since only it will be realized 
during the operation of the structure. To study the working process of energy-saving machines 
and their structural elements for the implementation of specific technological work processes, the 
following research sequence was established: 

— justification and development of calculation schemes for the machine: 

— analysis of the basic forms, amplitudes and frequencies of vibration of the structure; 

— determination of modes, rational parameters to create a new highly efficient design. 


8.2 INVESTIGATION AND DETERMINATION OF STRESSES AND STRAINS IN THE FORMING 
SURFACE OF THE STRUCTURE OF VIBRATION MACHINES 


The search for constructive solutions for the creation of machines with a variable amplitude-fre- 
quency oscillation mode is to search for an increase in the efficiency of the working process. 
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Revealing new phenomena in the operation of compacting machines and taking them into account 
when modeling work processes. Improvement of models that adequately meet the real conditions 
of motion of a vibration machine. An effective approach in this direction is the use of continuous 
models, which make it possible to take into account the propagation of waves, both in the design of 
a vibration machine and in the compacting mixture. This approach is the basis for determining the 
real distribution of amplitudes and frequencies of oscillations and the use of multi-mode effects. 

So, it is assumed that the design, subject to research and further creation, allows to trans- 
fer the maximum amount of energy from the working body to the medium due to the frequency 
spectrum with different energy components and the use of internal resonance phenomena of the 
system by choosing an operating mode that is coordinated in elastic-inertial and elastic-viscous 
properties of the machine and medium. 

Thus, on the basis of the created finite element model, studies of the structure were carried 
out in terms of natural frequencies and vibration modes. The calculation of frequencies and modes 
of vibration was carried out in two stages. At the first stage, according to the nonlinear theory, 
Static prestressing of the structure by all acting forces is performed, then at the second stage, 
the frequencies and, accordingly, the vibration modes are calculated. 

Studying the distribution of stresses and deformations in structural elements of vibration ma- 
chines will allow analyzing and identifying places of stress concentration, which allows optimization 
of design solutions. 

To study and determine the stresses and strains in accordance with the goal, the task of cal- 
culating the supporting structures of the system in its more critical states is posed. Therefore, for 
research, the design of a vibration unit was selected, which is a form-forming surface for compact- 
ing a concrete mixture. The structure consists of a welded box-section frame mounted on rubber 
resilient supports and a metal sheet as a shaping surface. Loads of auxiliary structural elements 
and concrete mix were included in the dynamic component of forces. 

The total number of finite elements turned out to be 09 856, the number of nodes — 100 998. 
It is assumed that the structure frame is rigidly fixed on rubber supports and the materials of all 
elements are deformed only in the elastic stage. 

The nature of the stress-strain state of the structure under the action of external forces and 
gravity is illustrated in color in Fig. 8.1. 

The scale of displacement values corresponding to this palette is represented on the left in the 
form of a colored column. As follows from Fig. 8.1 the maximum values of the equivalent stres- 
ses of the forming surface are 2.0 MPa. Higher values of stresses are concentrated in the area 
of application and distribution of the driving force. 

The frame design accepts a maximum stress of 11.7 MPa (Fig. 8.2). The area of action of such 
stress values is limited by the contact zones of the frame with the supports and is of a local nature. 
In general, the frame of the structure is in a state of uniform stress distribution. 

In order to study the behavior of the forming surface under dynamic action, the distribution of 
stresses in time is analyzed. In particular, Fig. 8.3 shows vibrograms of changes in the stresses 
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of the forming structure for individual elements located near the application of a dynamic external 
force. As seen from Fig. 8.3, at the beginning, a transient process was detected, caused by the 
exit of the system from a state of rest into oscillatory motion. An increase in the amplitude of 
stress changes at this stage is a consequence of the transition of the system through resonance 
modes at lower vibration frequencies. 


2027494, 
2213119, 
1898743, 
1584368, 


1269992, 
959617, 
641241, 
326866, 

12490, 





O Fig. 8.1 Equivalent stress distribution (SolidVonMisesStress) of the forming 
surface (time 0.0748 s) 


11673281, 
10216176, 
8759071, 
7301966, 
9844862, 
43877957, 
2930652, 
1473547, 
16442, 


O Fig. 8.2 Equivalent stress distribution (SolidVonMisesStress) of the structure 
frame (time 0.0748 s) 


In such modes, the forms of natural vibrations of the system are realized with large vibration 
amplitudes and, accordingly, a lower frequency. When the operating mode is set, negative stress 
arise in structural element 18 857. The nature of the change in these stresses is periodic and 
antiphase in comparison with the applied spatial force. The stress of element 48 345 are positive 
and vary with an offset of 2/2 relative to the applied force. This dependence is explained by the 
Spatial arrangement of the elements in relation to the point of application of the force. But the 
nature of the change in this force of the force load, since it changes along the two axes of space. 
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For structural elements 48 114 and 17 290, which are located symmetrically on the opposite 
side of the applied force, the nature of the change is shown in Fig. 8.4. In contrast to the previous 
vibrograms, one can see the absence of in-phase changes in stress in relation to the forced force. 
The significant difference can be explained by the dissipation phenomena occurring in the structure. 
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—— Solid XY Shear Stress (Element 48345) vs Set Value 


—— Solid XY Shear Stress (Element 18857) vs Set Value 
— sin(360*133.5*!x+sin(360*133.5*!x+90)(6); ScaleFactor: 100000) 


O Fig. 8.3 Dependence of shear stresses (Solid XY ShearStress) of the forming 
surface (elements 48 345 and 18 857) and the vertical component of the driving force 
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—— Solid XY Shear Stress (Element 48114) vs Set Value 

—— Solid XY Shear Stress (Element 17290) vs Set Value 

— sin(360*133.5*!x+ sin(860*133.5*!x+90)(6); ScaleFactor: 100000) 


O Fig. 8.4 Dependence of shear stresses (Solid XY ShearStress) of the forming 
surface (elements 48 114 and 17 290) and the vertical component of the driving force 


It is quite obvious that the change in the stress-strain state of the system under study is 
rather difficult. Such a change is associated with the implementation of complex spatial vibra- 
tions and the multi-mode of the sealing process. An important criterion for evaluating a design 
from the point of view of the efficiency of the compaction process is shear stress. Indeed, 
in the presence of such stresses in the medium, there is an intensive movement of particles 
and compaction. 
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To analyze the shear stresses of the forming surface, time intervals are selected that cor- 
respond to 1/4 of the period of oscillation of the driving force. Fig. 8.5 shows the distribution of 
stresses in the XY plane at a time of 0.0748 s, which is taken as a relative origin. 
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O Fig. 8.5 Distribution of shear stresses (Solid XY ShearStress) of the forming 
surface (time 0.0750 s) 


This is evidenced by the zones with the maximum positive stress values, in which the force is 
applied. The entire surface is in a complex sign-like stress-strain state. 

When the force is turned by an angle of 7/4, the stress distribution changes (Fig. 8.6). The 
appearance of additional zones of maximum stresses indicates the transfer of the load in the hori- 
zontal direction and the presence of bending and torsional vibrations. In this case, the middle of the 
structure is at relative rest. 
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O Fig. 8.6 Shear stress distribution (Solid XY ShearStress) of the forming 
surface (time 0.0768 s) 


When the force returns to an angle of 2/4, almost the entire surface of the forming surface 
is disturbed (Fig. 8.7). Deformation and stress have a sign-like character indicating the presence 
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of wave phenomena. The transfer of stresses occurs in the direction of the diagonal connecting the 
points of application of external dynamic forces. 
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O Fig. 8.7 Shear stress distribution (Solid XY ShearStress) of the forming 
surface (time 0.0786 s) 


In the period of time 0.0806 s (3/47), the propagation of oscillations continues, as evidenced 
by a change in the stress distribution (Fig. 8.8). Deformation and stress are alternating in nature, 
as in the previous period of time, but their stress distribution is more uniform and orderly. It should 
be noted that there are zones in the corners of the forming surface. In these zones, the change in 
stresses is less intense and most of the surface perceives positive stress values. At the same time, 
the number of surface areas with negative stress values is insignificant. 
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O Fig. 8.8 Shear stress distribution (Solid XY ShearStress) of the forming 
surface (time 0.0806 s) 


Comparing the obtained results of the stress-strain state of the forming structure, the follow- 
ing can be noted. The use of high frequencies of the operating mode will allow obtaining a new effect 
when compacting concrete mixtures. Due to the application of a spatial forced force to the forming 
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structure, a complex stress-strain state of metal structures is created. And direct contact with 
the concrete mix helps to reduce energy consumption for the compaction process. 

To study the motion of systems related to block structures, a model was adopted (Fig. 8.9), 
which is a two-mass vibration system. The lower mass of the system consists of four vibration 
blocks resting on elastic elements calculated according to the vibration isolation conditions. 

The upper mass simulates the shape with the medium to be processed. Between the upper and 
lower masses there are elastic elements, the stiffness of which is taken from the condition for the 
realization of the resonant mode of motion and providing the technological process. 

In order to carry out research, the following overall dimensions of the structure were adopted: 

— total length 3.30 m; 

— width 0.50 m: 

— height 0.34 m. 

The model was investigated by the finite element method. The finite element model of the 
structure was compiled by approximating all load-bearing elements, including the forming surface, 
with two-dimensional finite elements of the PLATE type, elastically deformed under the action of 
a longitudinal force, bending moments in two planes and a torque. 

Vibration isolation supports and elastic elements of the model are adopted as three-dimen- 
sional CE of the SOLID type, since the processes occurring in such structural elements are more 
complex in terms of energy dissipation. 

The total number of finite elements was 19 258, the number of nodes was 19 912, and the 
total number of searched variables was 20 928. 

It is assumed that at the extreme points, the supports rest on the foundation and are fixed, 
therefore, in the model, the extreme nodes of the supports are restrained in the X, Y, Z directions, 
and rotations in all three axes are also prohibited. 

It is believed that the materials of the entire structure are deformed only in the elastic stage. 

When checking the strength and determining the functioning of a system that can be subjected 
to loads that are different in their nature and applied in various combinations, it is advisable to first 
analyze the behavior of the structure under the simplest loads, and then complicate them with the 
addition of other disturbances. 





O Fig. 8.9 Model of block vibration unit 
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To check the model and assess its functioning under the action of static and dynamic loads, the 
following research stages were taken: 

— static calculation under the action of the weight of structural elements and the equivalent 
load of the medium being processed; 

— modal analysis: 

— dynamic analysis of the structure under the action of a forced force with a vibration frequen- 
cy close to the natural vibration frequency. 

The studies of the unit (Fig. 8.10) were carried out in several stages. at the first stage, 
a modal analysis was carried out in order to determine the main forms and the corresponding 
vibration frequencies. 





O Fig. 8.10 General view of a 3D model of a block vibration unit with directional vertical vibrations 





O Fig. 8.11 First vibration mode (Mode 1, f=1.14 Hz) 


50, at an oscillation frequency f=1.14 Hz (Fig. 4.12), f=1.32 Hz (Fig. 8.12), the vibration 
unit will vibrate in the first and second modes of vibration in the horizontal plane (XY). This mode 
of vibration is due to the value of the stiffness of elastic supports, which in the horizontal direction 
have an insignificant stiffness due to geometric features and physical and mechanical characteristics. 

Fig. 8.13 shows the operating mode of the unit at an oscillation frequency f=2.97 Hz. In this 
mode, vibrations are performed in the horizontal plane and have a twisting character. 

When the mode is realized at a frequency of f=4.10 Hz (Fig. 8.14), the oscillations are di- 
rectional in the vertical direction. This mode corresponds to the resonant mode of vibration, when 
the vibration frequency corresponds to its own frequency on vibration isolation elastic elements. 
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O Fig. 8.12 Second vibration mode (Mode 2, f=1.32 Hz) 


O Fig. 8.13 Third vibration mode (Mode 3, f=2.57 Hz) 


O Fig. 8.14 Vibration mode (Mode 8, f=4,10 Hz) 


The implementation of oscillations according to the mode at a frequency of f=16.61 Hz (Mode 23) 
is a consequence of the influence of the rigidity of the form (Fig. 8.15). Oscillations occur due to the 
bending of the structure in the YZ plane. Such modes are insufficient for the process of compaction 
of concrete mix. Purely unidirectional oscillations with low frequencies do not cause stresses in the 
medium that do not exceed the voltage loss of continuity. The operating mode, in which the phase 
displacement will be realized, is possible when increasing to a frequency of f=27.40 Hz (Fig. 8.16). 
This is due to the fact that with increasing frequency, the numerical values of the compressive 
stresses of the mixture layers and the rupture stress of the layers increase significantly. Also, this 
mode implements vibrations of the shaping surface, which are realized on elastic elements. 
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O Fig. 8.15 Vibration mode (Mode 23, f=16.61 Hz) 


O Fig. 8.16 Vibration mode (Mode 29, f=27.40 Hz) 


A further increase in vibration frequencies and the implementation of the corresponding vibration 
modes leads to vibrations of the vibration unit form (Fig. 8.17), which is explained by the sufficient 
stiffness of the form frame and the corresponding stiffness of elastic elements (rubber gaskets). 

Since the operating mode of such structures is considered as close to resonance, the choice 
of the operating mode must be carried out in such a way that the value of the operating frequency 
is 0.90—0.95 of the natural frequency. 

Thus, the operating mode is adopted with a vibration excitation frequency of 29 Hz. 

The implementation of the dynamic analysis was carried out in order to obtain the distribution 
of the vibration amplitudes of the shaping surface along the length of the structure with the appli- 
cation of a driving force on each vibration block. 


O Fig. 8.17 Vibration mode (Mode 31, f=35.72 Hz) 
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To analyze the dynamics with directional vibrations, the values of displacements in time 
were obtained. 

To study the behavior of the shaping surface, a list and place of points on the surface of the 
form were determined with the determination of the distribution of vibration amplitudes depending 
on the location of blocks with vibration exciters (Fig. 8.18). 

Thus, the implementation of oscillations in one period is shown in Fig. 8.19. 





O Fig. 8.18 Arrangement of control points on the surface 
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O Fig. 8.19 Distribution of vertical vibrations along the length of the forming surface 


As can be seen from the distribution, the oscillations are carried out in the vertical plane, which 
fully confirms the classical approach to modeling such systems with discrete parameters when con- 
sidering a one-dimensional problem. The maximum value of the vibration amplitude is 0.22—0.25 mm. 
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50, the implementation of oscillations for one period is shown in Fig. 8.20. As can be seen from 
the distribution, oscillations occur in the vertical plane, which fully satisfies the classical approach 
to modeling such systems with discrete parameters when considering a one-dimensional problem 
and thus confirms the adequacy of the created model within the framework of the research. The 
maximum value of the vibration amplitude is 0.22—0.25 mm. 

The analysis of vibrations in the horizontal direction was carried out relative to the 
Y-axis (Fig. 8.20), showed the practical absence of vibrations and displacements of the forming 
surface. The vibration amplitudes are 2-10°’ m, which is three orders of magnitude less than in 
the vertical direction and it is quite obvious that such values can be neglected. 
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O Fig. 8.20 Distribution of horizontal vibrations along the length of the forming surface 


Likewise, the absence of oscillations along the X axis, since the forcing force is directional and 
does not act in this direction. 

The considered mode of operation of the block structure of the vibration platform with direc- 
tional vertical vibrations made it possible to analyze the distribution of amplitudes along the length 
of the forming surface and confirm the adequacy of the adopted model. 

To implement the idea of a polyphase oscillation mode, studies of a finite element model have 
been carried out. Modeling of the vibration unit movement was carried out under two conditions 
of application of external forces: the first option provided for the application of a driving force with 
a shift angle of 90° on each vibration unit (Fig. 3.2, o,=0; o, =90°; @, =180°; @,=270°); the 
second variant of the load with angles ,=0; », =60°; @,=120°; @,=180°, respectively. The 
frequency and value of the driving force were similar to the case of directional vertical vibrations. 
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SO, as a result of dynamic analysis, distributions of vibration amplitudes along the Z (Fig. 8.21, a) 
and Y (Fig. 8.21, &) axes were obtained for the angle @ =90°. 
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O Fig. 8.21 Distribution of vibration amplitudes along the length of the structure 
for one period of vibration (angle=90°): a —X axis; b — Y axis 


The maximum value of the vibration amplitude within the framework of the studies performed was 
obtained at the boundary of the form and, respectively, are 0.6 and 1.1 mm in time intervals from 3/47 
to 7/4r for oscillations of the Z axis, and for oscillations in the Y direction — 22-10 and 5-10 mm. 
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It should be noted that there are points at which vertical vibrations occur to a much lesser extent, but at 
the same points there are vibrations with a horizontal component, therefore, at these points vibrations 
are present in both directions and there are all the conditions for compaction of the building mixture. 

Distribution of vibration amplitudes at angle @ =60° is shown in Fig. 8.22 in practice, it is 
possible to see a similar picture with the previous study, although there is a certain difference, 
which consists in the fact that there are no points with less vibration intensity. 


0.000800 
0.000600 
0.000400 


0.000200 
e 


N 


0.000000 
~0.000200-—Z 
~0.000400 
0.000600 








~0,000006 
~0.000009 
~0.000012 
& 9.000015 
~0.000018 
-0.000021 


—0.000024 














O Fig. 8.22 Distribution of vibration amplitudes along the length of the structure 
for one period of vibration (angle=60°): a —X axis; b — Y axis 
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So there are places where the displacements pass through the zero mark, but compared to the 
previous mode of operation, these points do not have a clear position along the length of the form, 
but shift with each 1/4 time interval. Thus, transferring vibrations to a larger volume of concrete 
mix near the contact zone. It should be noted a certain increase in the amplitude of oscillations in 
the horizontal direction in comparison with the shear angle @ =90°, as for vertical oscillations, that 
is, a decrease in the values to 0.6 and 0.8 mm. 

Comparing the three operating modes of the block platform (@ =0°, ~=60°, » =90°), it is 
possible to note a positive contribution to the qualitative characteristic of oscillations when changing 
phase angles other than 0. And when comparing the general distribution of amplitudes at an angle 
of ~=90° and m=60°, it is obvious that it is the shift angle @ =60° that has the advantage. 

For a more detailed study of the motion of the shaping surface with this implementation of the 
model (Fig. 8.10), let's consider modeling the motion in time. 

The movement of the surface at points 1, 2, 16, 17 (Fig. 8.18) in time indicates the presence 
of different values of the vibration amplitudes (Fig. 8.23). So, at point 16, the amplitude of oscilla- 
tions reaches 0.6 mm, and at point 2, at the same time, it has a value of 0.2 mm, which is evidence 
of the presence of a wave process. 

At other points (Fig. 8.24) vibrograms of movement are identical in character to those dis- 
cussed above. In terms of magnitude, the amplitudes of the oscillations are of lesser importance. 
At point 14, the magnitude of the amplitude reaches 0.4 mm, and at point 3, in the same period of 
time, the amplitude of the oscillations is 0.1 mm. At points 14 and 15, there is practically no phase 
displacement between the vibration amplitudes. 
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O Fig.-8.23 Vibrograms_of movement_of points 1,2, 16,17 in time 
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Vibrograms of movement at points 0, 6, 12, 13 in time (Fig. 4.25) are similar in nature, as for 
other points (Fig. 8.24, 8.25). 
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O Fig. 8.24 Vibrograms of movement of points 3, 4, 14, 15 in time 
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O Fig. 8.25 Vibrograms of movement of points 5, 6, 12, 13 in time 


However, there are also certain differences. So, it was recorded that at points 6, 12 and 13 
the vibration amplitudes have similar numerical values within 0.18-0.23 mm. The change in the 
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oscillation amplitudes at points 7, 8, 9, 10, 11 in time (Fig. 8.26) is characterized by a certain 


deviation in numerical values and in phase. 
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O Fig. 8.26 Vibrograms of movement of points 7, 8, 9, 10, 11 in time 


At these points, the amplitudes of oscillations reach values of 0.27—0.40 mm, that is, the limit 
of change is the smallest in comparison with the movement of points 1, 2 and 16, 17 (Fig. 8.26). 

The research results indicate the presence of vibration amplitudes different in shape and nu- 
merical values over the area of a vibration unit with a polyphase vibration spectrum. 


8.3 INVESTIGATION OF THE STRESS-STRAIN STATE OF THE FRAME AND FORMS OF A VIBRATION 
UNIT WITH SPATIAL VIBRATIONS 


8.3.1 MODELING THE STRESS-STRAIN STATE OF VIBRATION PLATFORMS 


Determination of the stress-strain state of the vibration platform frame is not an easy task. 
The complexity of the research is associated with various factors: an increase in the number of 
finite elements, the use of various types of finite elements in one numerical model, problems of 
mechanical properties of rubber platform supports, modeling of various types of loads (static and 
dynamic), modeling of welded joints of a space frame [8, 9]. 

The following types of finite elements can be used to calculate frames and shapes: 

a) the frame elements are modeled with SOLIDS5 elements. The element is a three-dimensio- 
nal hexagon with eight nodes located at the corners. This type of element allows the use of 
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isotropic, plastic and other nonlinear types of materials, with the exception of hyperelastic ones. 
Elements of this type can take degenerate shapes. Degenerate elements have at first a quadran- 
gular or hexagonal shape, but by combining two or more nodes that take the shape of a triangle for 
flat elements or a tetrahedron, prism or wedge for volumetric ones; 

b) welded seams are taken as idealized and modeled from volumetric finite elements in the form 
of a six-node triangular prism with the same bases equal to the leg; 

c) to model rubber supports, a finite element SOLID185 is used — an element similar to 
SOLID55, but allows the use of hyperelastic materials. 

The displacements of the lower planes of the supports are limited in the X, Y, Z directions, 
which simulates a rigid clamping of the structure in the foundation. The vibration exciters on the 
frame are shown conditionally, and their effect on the frame is determined by the driving force. 
The load on the frame from the form and the concrete mix is simulated by applying a uniformly 
distributed load. 


8.3.2 VIBRATION PLATFORM WITH VERTICAL VIBRATOR 


A vibration platform with a horizontal vibrator is divided into finite elements (Fig. 8.27). For it, 
a displacement distribution diagram was obtained (Fig. 8.28) in the X, Y, Z directions and a stress 
distribution diagram (Fig. 8.29). 

The largest total displacements of the frame elements (Fig. 8.28) are observed in the outer 
frame elements, their maximum value is 0.0615 mm. 
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O Fig. 8.27 Final element diagram of a frame with a vertical vibrator 
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O Fig. 8.28 Distribution of total displacements: 
a — frames with a vertical vibrator; b — frames with 
a vertical vibrator (the top sheet is conventionally not shown) 


Movement in the vertical direction — along the z-axis, has both positive and negative va- 
lues (Fig. 8.29), which indicates the presence of small deflections (0.0569 mm) on the cantilever 


sections of the frame and local bends on the side channels of the frame. 
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b 
O Fig. 8.29 Distribution of vertical displacements: 
a — frames with a vertical vibrator; b — frames with 
a vertical vibrator (the top sheet is conventionally not shown) 


The highest values of the equivalent stress reach 10.6 MPa (Fig. 8.30). They are observed 
in the places where vertical and horizontal channels are welded at the place where the vibrator 
is attached. The values of these stresses are less than the maximum standard stresses, which 
indicates that the strength of the welded seams is ensured. 
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O Fig. 8.30 Distribution of equivalent stresses: 
a — frames with a vertical vibrator; b — frames with 
a vertical vibrator (the top sheet is conventionally not shown) 


In the direction of the X axis, the distribution of displacements is circular in nature, with 
a maximum (0.042 mm) at the place of load application; the outermost parts of the frame are 


the least deformable. 
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The distribution of stresses in the frame elements is uneven. The concentration of stresses 
is observed in the places of welding of elements, but the stress raisers have small values of 
4.12 MPa in comparison with the ultimate strength of steel. 

The stress in the direction of the X axis is distributed circularly, similar to the distribution 
of displacements. Moreover, the highest stress values are observed in the two nearest cells to 
the place of load application. 


8.3.3 VIBRATION PLATFORM WITH HORIZONTAL VIBRATOR 


The structure of the vibration platform VPGP-6.3x3.6 (Fig. 8.32) includes a movable frame 
with overall dimensions in plan 6.28x3.66, made of channels and a steel sheet made of calm-melting 
steel st.3sp. In the windows of the moving frame, vibration plates with welded-on M30 nuts are rigidly 
attached. On the moving frame there are rigid transverse stops for fixing the form with the concrete 
mixture. The vibration exciter is bolted to the movable frame. With the help of a V-belt transmission, 
the vibration exciter shaft rotates from an electric motor with a power of 30 kW. The windows of 
the movable frame are closed with metal covers that prevent particles of the concrete mixture 
from entering the platform's vibration drive. The movable frame rests on 8 rubber-metal supports 
attached directly to the foundation of the vibration platform. The supports are attached to the moving 
frame with the help of protrusions on the supports entering the mounting holes in the moving frame. 
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O Fig. 8.31 Vibration platform VPGP-6.3x3.6 for the formation of ribbed slabs 
of coatings with the dimensions of the movable frame in the plan of 6.3x3.6 m: 
1 — movable frame, 2 — elastic support, 3 — vibration exciter of torsional 
vibrations, 4 — electric motor, 5 — foundation 
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The vibration platform model was loaded with the same uniformly distributed load, but the 
design and direction of vibration movement were changed. 


For this option, a deformed diagram, a displacement distribution diagram (Fig. 8.32) and 
a stress distribution diagram (Fig. 8.33) were also obtained. 
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O Fig. 8.32 Distribution of total displacements: 
a — frames with a horizontal vibrator; b — frames with 
a horizontal vibrator (the top sheet is conventionally not shown) 
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The most total displacements of the frame elements are 0.0782 mm (the area is marked in red 
in the diagram). Such deflections occur in the side channels of the frame on the cantilever sections 
free from supports (Fig. 8.33). In the center of symmetry of the frame, local deflections of up 
to 0.0608 mm are also observed. 
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O Fig. 8.33 Distribution of equivalent stresses: 
a — frames with a horizontal vibrator; b — frames with 
a horizontal vibrator (the top sheet is conventionally not shown) 
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Equivalent stresses have a maximum value of 16.6 MPa (Fig. 8.33). These values can be seen in 
the channels located above the frame supports and in the welds at the vibrator attachment points. 

Dynamic calculation of the form [10] (Fig. 8.34) is performed in the ANSYS/LS-Dyna finite ele- 
ment analysis program. The program simulates the compaction of a mixture in a metal form during 
the operation of a shock-vibration unit. For the numerical model, two types of finite elements are 
used that implement the elastic and plastic properties of the metal form and mixture (Fig. 8.35). 

The force resulting from the fall of the form from a height of 8 mm, taking into account the 
acceleration of gravity and the static and dynamic action of the loader on the metal form, is taken 
as the maximum load. 





O Fig. 8.34 Metal form 


ELEMENTS 


O Fig. 8.35 Numerical model of the form with wood concrete mixture 


169 


DYNAMIC PROCESSES IN TECHNOLOGICAL TECHNICAL SYSTEMS 


Fig. 8.37 shows the picture of stress distribution on a deformed structure diagram. The loader 
is not shown in the Fig. 8.37, since the distribution of stresses and strains in it is not important 
for the task at hand. As can be seen from the Fig. 8.37, the stress in the elements of the form 
is on average 10-60 kPa, and the subsidence of the concrete mixture during compaction is about 
38 mm, as shown in the work. 
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O Fig. 8.36 Stress distribution in structural members 


O Fig. 8.37 Stress concentration at the corners of the form 
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The greatest deformation and stress concentration (60-65 kPa) occurs in the corners of the 
form, where the welded seams of the structure are located (Fig. 8.37). 


8.4 INVESTIGATION OF THE STRESS-STRAIN STATE OF CRANES METAL STRUCTURES 


The study of static dynamic loads of the slewing ring of a truck crane was carried out in order 
to establish the position of the system and its elements with the highest stresses. The basic laws 
of the theory of elasticity and plasticity were used to construct a computational mathematical 
model of a loaded support-rotary contour. Equations of motion were drawn up taking into ac- 
count the behavior of structural elements and the machine as a whole under the simultaneous 
action of various loads. The design model was taken in the form of a structure consisting of rods 
rigidly connected to each other at the nodes. For each bar, this section is shown with geometric 
characteristics. As a calculation method, the finite element method was used with the definition 
of displacement (deformation), force (stress) at the mesh nodes of structural elements of the 
support-rotary contour. The distribution of stresses for each structural element was obtained 
depending on the type of load. 

In the theoretical studies [7—11], a scientific hypothesis was proposed according to which the 
development of reliable and most efficient machines for various technological purposes is based 
on the establishment and general purposeful use of the regularities of changes in the internal 
properties of the working bodies of machines and the technological load, which is a processing 
material. The implementation of the hypothesis is provided by the solution of a scientific idea [11]. 
Recently, the development and creation of the design of machines for various technological purpo- 
ses is carried out in the direction of searching for design solutions with a changing operating mode. 
In the process of operation of such machines, the maximum allowable use of the internal resources 
of the machine structure is assumed, which in turn necessitates ensuring strength and reliability 
with specified dynamic parameters. An integral assessment of the state of metal structures of the 
finished product can be carried out on the basis of dynamic tests. The idea of such tests is to find 
the actual dynamic characteristics of the structure (natural frequency, vibration amplitude, etc.) 
with the subsequent comparison of these characteristics with the characteristics obtained by the 
method of mathematical modeling of the given structure system. The mathematical model should be 
created in such a way as to adequately describe the real structure and its behavior under various 
loads, and also be able to reflect the modeling of various kinds of imperfections associated with 
manufacturing and direct operation as intended. 

Under dynamic loads on structural elements inherent in most machines for technological pur- 
poses, the main type of damage is the formation and development of cracks due to the accumula- 
tion of fatigue stresses. Durability in case of tedious destruction is determined by the long-term 
cyclicity of the action of loads on the supporting structures of the machine and its individual ele- 
ments, as a result of which a crack is formed, which gradually develops and leads to destruction. 
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The solution to the problem of determining the study of the stress-strain state is carried out 
by using an integrated approach that combines analytical calculations, mathematical modeling and 
experimental research directly. 

For the implementation of these tasks, the following sequence of research was envisaged: 

— analysis of the behavior of structural elements of the machine from the point of view of 
taking into account all types of loads that were carried out during its operation and determining the 
possible combination of loads acting on the elements; 

— development of a computational model of the research object, taking into account the ge- 
neral or individual most loaded nodes or structural elements: 

— drawing up equations of motion, taking into account the behavior of structural elements and 
the machine as a whole under the simultaneous action of various loads; 

— assessment of real loads for the development of an algorithm for determining the stress- 
strain state, identification of the most stressed structural elements; 

— development on a computer model of a matrix of control points for the limiting values of the 
integral characteristics of the state of the structure for subsequent use in field tests: 

— carrying out computer calculations to determine stresses and strains by applying a number 
of possible loads to the model; 

— development and implementation of experimental studies on samples reflecting the real ope- 
rating conditions of a part or machine unit; 

— adjusting the computer model until the comparison of the integral characteristics obtained by 
measurements at control points during the experiment and during simulation can differ from each 
other within the permissible error (the computer model obtained in this way will be adequate to the 
real design within the limits). points of adequacy — points of control of integral characteristics): 

— development of recommendations for improving the existing or creating a new design, cor- 
responding to the basic idea of the equal strength of stress distribution, regardless of the applied 
forces on the elements of the loaded product. The analysis of possible loads was carried out on the 
basis of the studies carried out using the example of truck cranes. This is due to the fact that these 
machines are a classic example of the structure of a system in which static and dynamic loads are 
organically combined, which require correct accounting for safe crane operation under real load 
conditions. The basic laws of the theory of elasticity and plasticity were used to construct a com- 
putational mathematical model of a loaded SC [11-14]. Equations of motion were drawn up taking 
into account the behavior of structural elements and the machine as a whole under the simultaneous 
action of various loads. In the case of a dynamic effect on its metal structure (Fig. 8.38) as a linear 
elastic body, the equations of motion in the adopted coordinate system (Fig. 8.39) are as follows: 


(8.21) 


where z,,. — stress tensors, (/, /=1, 2, 3); F, — components of volumetric forces; p — material 


density; U. — displacement vector. 
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O Fig. 8.39 Coordinate system of the acting forces 
on the support contour of the truck crane 


If to take into account the viscous properties, then let’s replace the operator from (21): 


a: a° 0 
-—esp-| — J+k—, 22 
are P (z+ dt oie 


where k — viscosity coefficient. 
Equation (22) in this form remains unchanged. 
Let's introduce designations for the speeds of wave propagation in the SC: 


9 1—-v E 


Ep 1 E (8.23) 
(1+v)(1-2v) p * 2(1+v) p | 


For the case when the force F(t) is ordinary and acts along the OX1 axis (at the origin of the 
coordinate system): 
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A pf e4)44 (1-2 ofc 4 )-[1-9%2). farte—aajal cea 
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xx, | 1 R) 1 7 

Ui OG Xiah = ie 2 mata Cara ata Cara Eee (8.25) 
1 
C; 


If the force acts along the OXe axis, then it is necessary to replace the indices: 


1 2; 
y (22) => 42-3; 
35> 1. 


If the force acts along the OX8 axis, then it is necessary to replace the indices: 


139 3; 
y(22) => 42-1, 
332. 


For a plane problem (x3=Q); for a one-dimensional problem: (x1 #0; x2=x3=0). 

The SC design scheme is taken in the form of a structure consisting of rods rigidly connected 
to each other at the nodes (Fig. 8.40) in a rectangular coordinate system X, Y, Z and is built along 
a line connecting the centers of gravity of the cross-sections of the rods. 

Each element is designated by small circles numbered with Arabic numerals. The serial num- 
ber of the bar (element) is shown in a circle. For each bar, this section is shown with geometric 
characteristics. As a calculation method, the finite element method was used with the definition 
of displacement (deformation), force (stress) at the nodes of the mesh of structural elements of 
the SC structure (Fig. 8.41). 

The correctness of the results of mathematical modeling largely depends on the mecha- 
nical characteristics of structural materials and the study of the operation of complex assem- 
blies of load-bearing structural elements. Experimental studies of the connection of parts to each 
other (Fig. 8.42, a) and nodes (Fig. 8.42, 6) were carried out by the test method for determin- 
ing the mechanical characteristics to check the results of analytical calculations. Tests on these 
stands allowed to propose the design of the elements used in real structures of the defense 
industry complex. 


174 


8 RESEARCH OF STRESS-STRAIN STATE OF ELEMENTS OF TECHNOLOGICAL TECHNICAL CONSTRUCTIONS 














O Fig. 8.40 Design diagram of the frame of a truck crane 


00323 


O Fig 8.41 Stress-strain state of the frame of a truck crane 


It is found that the maximum loads occur in the position of the crane unit perpendicular to the 
main axis of the vehicle on which the crane unit is fixed. At the same time, it is found that of all 
the nodes of the crane unit, the greatest stresses arise in the support contour (SC), which en- 
sures the rotation of the lifting mechanism by 360 degrees. 
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O Fig. 8.42 Testing of structural elements of metal structures: 
a — tensile welded joint; b — bracket test 


By changing the dynamic and static loads and analyzing the results obtained, similar to Fig. 8.41, 
the research results were introduced into the SC structure in such a way that all the elements had 
the same order of numerical values of stresses, which in the final design solution made it possible 
to reduce the SC mass by 30 %. 


8.9 DISCUSSION OF RESEARCH RESULTS 


The research results (Fig. 8.1) indicate that the maximum values of the equivalent stresses 
of the forming surface are 2.0 MPa. Higher values of stresses are concentrated in the area of 
application and distribution of the driving force. The frame design accepts a maximum stress of 
11.7 MPa (Fig. 8.2). The area of action of such stress values is limited by the contact zones of 
the frame with the supports and is of a local nature. Vibrograms of changes in the stresses of the 
forming surface are given for individual elements located near the application of a dynamic external 
force. This is a fundamentally new result. It consists in the fact that the detected transient process 
is intended to be taken into account when determining the parameters and locations of the vibra- 
tors. In such modes, the forms of natural vibrations of the system are realized with large vibration 
amplitudes and, accordingly, a lower frequency. And this opens up a real opportunity to reduce the 
energy consumption of the drives of the vibration machine. The stresses for structural elements 
located symmetrically on the opposite side from the applied force are investigated (Fig. 8.5). 
In contrast to the previous vibrograms, there is no in-phase change in stresses in relation to the 
forcing force, which confirms the inclusion of dissipative forces in the design model. An important 
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criterion for evaluating a design from the point of view of the efficiency of the compaction pro- 
cess is shear stress. Indeed, in the presence of such stresses in the medium, there is an inten- 
sive movement of particles and compaction. To analyze the shear stresses of the forming struc- 
ture, time intervals were selected corresponding to 1/4 of the period of oscillation of the driving 
force (Fig. 8.6). The zones with the maximum positive stress values, in which the force is applied, 
have been determined. A complex sign-like stress-strain state of the structure under the action 
of a spatial load has been confirmed. This study has a limitation in terms of the adopted law of the 
change of dissipative forces. 

A fundamentally new result has been obtained, which consists in the fact that the implemen- 
tation of a complex vibration mode as an effective method of accelerated compaction of concrete 
mixtures is ensured by the arrangement of unbalances at a certain angle on each separate vibration 
unit. AS a result, energy consumption is reduced by 30 %, and the process of forming a concrete 
product is reduced by 20 %. The presence of different forms of the polyphase spectrum is con- 
firmed by the vibration modes (Fig. 8.11-8.17) and the distribution of the vibration amplitudes 
of the form surface along the length of the structure for one period of vibration. The use of such 
effects is determined by the overall dimensions of the product in terms of plan and height, which 
affects not only the phase deployment of the unbalances along the central axis of the vibration unit, 
but also the magnitude of the static moment of the unbalances. 

The design diagram of an energy-saving vibration plant with a polyphase vibration spectrum has 
been substantiated. Modeling of the working process of a vibration plant is based on the use of the 
finite element method using the MSC.NASTRAN calculation complex (MSC.Software, Germany). 

As aresult of the research, the main modes of vibration have been determined by the numerical 
values of the frequencies of 1.32 Hz, 4.10 Hz, 16.61 Hz, 27.40 Hz. The calculated value of the 
oscillation frequency of the operating mode is 20 Hz. The most effective in the framework of the 
research carried out is the vibration form of 24.31 Hz. The implementation of this form of vibration 
provides an energy-saving mode of operation of the vibration unit. 

A polyphase spectrum of vibrations of the unit has been established and proposedat an exci- 
tation frequency of 20 Hz, according to which vibration amplitudes of 0.27...0.6 mm are realized. 
The investigated vibration mode of the computational model is a confirmation of the energy-saving 
mode, which is a prerequisite for calculating and creating a new class of vibration unit. 

Analysis of the stress-strain state of frames with horizontal and vertical vibrators shows that in 
the elements of the proposed frame structures during the operation of vibration platforms, stress- 
es arise significantly less than the design resistance of steel, which provides the condition for the 
strength of the structure. The deflections in the frame elements are much less than permissible, 
which indicates the sufficient planar and spatial stiffness of the created frames and metal forms. 

Frames with a vertically positioned vibrator work in tension and bending. The movements of the 
frame members are distributed circularly from the vibrator location. The greatest stresses arise 
at the welding points of the longitudinal and transverse channels and the top sheet in the vibrator 
attachment area. Their maximum values are 11.4 MPa. 
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For the distribution of displacements of frames with horizontally located vibrators, areas of 
both positive and negative deflections are characteristic. At the junction of the vibrator with the 
frame, a wave-like bend is observed, and in the cantilever sections, deflections of 0.0672 mm are 
observed. The maximum stress in the vertical direction is 16.6 MPa. 

In the manufacture of frames of such structures, special attention should be paid to the quality 
of welded seams in the area of the vibrator and in the attachment points of the supports in order 
to minimize the possibility of cracks in them. 


CONCLUSIONS TO SECTION 8 


1. In the course of the research, finite element models of vibration machines have been de- 
veloped to implement the transfer of energy in the vertical direction with a fixed frequency, with 
a polyphase direction of external forces, with a spatial direction of movement in the vertical and 
horizontal planes. 

2. The study of the nature of the distribution and the numerical values of stresses and strains 
has been done in the forming structure of the unit, depending on the angle of instantaneous action 
of the external force of the vibrators and vibration frequencies. 

3. The maximum values of the equivalent stresses of the forming surface are 2.0 MPa. Higher 
values of stresses are concentrated in the zone of application and distribution of force, the maxi- 
mum stresses are 11.7 MPa. 

4. The results obtained indicate the presence of vibration amplitudes of different shapes and 
numerical values over the area of a vibration unit with a polyphase vibration spectrum. 

9. otress concentrators in structural elements of machines for technological purposes have 
been determined, taking into account the working loads. 

6. The static and dynamic loads of the support contour of a truck crane have been investigated. 

7. The obtained research results are taken into account when making further design decisions 
when creating such machines. 
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